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ABSTRACT 


This paper describes the present program, recent progress, and 
future plans of the U. S. Naval Radiological Defense Laboratory (NRDL) 
in the field of nuclear-debris formation. 

Two fallout-formation prediction systems have been programmed 
for computer calculation. The predictions represent great advance- 
ments over systems that do not account for fractionation, but much 
room remains for improvement. The greatest: needs appear to be for 


the definition of particle-size distribution, the accounting for agglom- 


eration effects, and the development of a kinetic approach. Preliminary 
steps in the development of a kinetic approach are described. 

A by-product of the prediction system is the calculation of unfrac- 
tionated radiological properties. As an example of this work, decay 
rates, dose rates, and spectra are given for the most important case 
of *38y fission produced by a thermonuclear neutron spectrum. 

Some results of correlating radiochemical data from fractionated 
samples of air- -burst debris are discussed, but numerical results are 
reserved for a following paper* to facilitate comparison with correla- 
tions from other burst types. 

Preliminary laboratory investigations of fission- -product interac- 
tion with various substrate materials are described. These show large 
negative deviations from Raoult’s law, particularly when the interac- 
tion is between a basic fission product and acidic substrate or an acidic 
fission product and basic substrate. 


'*See paper by Glenn R, Crocker, Francis K. Kawahara, and Edward C, 
Freiling, this volume, 
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Finally, the directions that we expect future efforts to take are 
indicated. These include a sensitivity - analysis of our ‘prediction sys- 
tem, extended fractionation correlations, model refinements, model 
extensions to new conditions, the retrieval of information from old 
debris, and laboratory studies of fallout formation under irreversible 
conditions. 


INTRODUCTION 


We have been asked to describe in this paper the project on studies 
of nuclear-debris sees that: we are carrying out at NRDL for the 
Fallout Studies Branch, U. 8S. Atomic Energy Commission (AEC). This 
project is entitled “The athe Distribution and Characteristics of 
Nuclear Debris.” However, if we were to limit this paper strictly to 
this part of our nuclear-debris program, we would present a very in- 
complete picture of our efforts. We have therefore chosen to preface 
this paper with a few remarks :on. what we are trying to do and for 
whom we are doing it. Figure .1..is a schematic presentation of the 
Physical Chemistry Branch, NRDL, fallout studies program at the 
present. and in the recent past. Other, related work-at NRDL is shown 
only where necessary. The various boxes are-coded to identify spon- 
sorship. Next to the AEC, the Defense Atomic Support Agency (DASA) 
sponsors most of our work. The focal point, the ultimate goal, is pre- 
diction. capability. This is indicated by the box labeled “Model of For- ; 
mation and Properties of Nuclear Debris.” As input this box requires 
basic data from nuclear and: physical chemistry, .in addition to empir- 
ical relations such as scaling functions and fractionation correlations. 
The scaling. functions come from laboratory experiments and the 
fractionation correlations.from.the analysis of nuclear debris. In addi- 
tion to NRDL, data come from various sources. . 

‘This paper reports primarily, but not exclusively; on research 
completed or ‘nearly completed.in the four following subprojects of our 
work for the AEC: 


Prediction of fractionation effects 
‘Prediction of radiological properties 


Fractionation correlations 


Laboratory studies 


Also included in each section is an indication of the future direction we 
visualize these studies will take. : 
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Fig. 1—Nuclear-debris study program. 
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PREDICTION OF FRACTIONATION EFFECTS 


Clarification of Concepts 


If the reader will permit. a heterodox interpretation of the title of 
this subproject, it is appropriate to begin this discussion with the ef- 
fect of fractionation on thought and communication—the semantic 
effect, if you will. With a realization of the existence and significance 
of fractionation, many terms previously used to describe and discuss 
fallout have become undefined and meaningless. Thus phrases like 
“bomb fraction,” “fissions in the sample,” and “kilotons per square 
mile” are no longer “OK terms.” Moreover, as our understanding of 
fractionation has increased, deviation of fission-product abundances 
from those expected in instantaneous thermal-neutron fission of asbyy 
and of decay rates from the Way—Wigner t™'? rule are no longer con- 
sidered significant criteria or acceptable measures of fractionation 
effects. We have considered it an important part of our effort to sug- 
gest new definitions and, where possible, operational definitions for 
many of these terms. These are listed in the appendix to this paper. 
The essential features of some ofthese definitions are soundly grounded 
in the fundamental properties of the debris. Accidental features (e.g., 
the choice of standards for measuring fractionation) can be varied if 
sufficient justification appears for doing so. Other terms and defini- 
tions are merely matters of convenience, and these are appropriately 
indicated in the appendix. The adoption of these or similar definitions 
will be of great assistance, not only in promoting clear thinking and 
improving the signal-to-noise ratio when fractionation is discussed but 
also in planning analyses of those properties of nuclear debris af- 
fected by this phenomenon. Further background and justification for 


these definitions can be found in our set of reports with the series title , 


‘Fractionation.’’!~4 


First Phase of Predicting Fractionation Effects 


Early in the program we came to the conclusion that the prediction 
of fractionation effects could best be carried out in two concurrent 
phases.’ The first of these phases consisted of rapidly assembling a 
prediction system that would give answers to urgent questions in as 
reasonable a manner as time allowed. The second phase consisted of a 
more fundamental long-range approach. The effects of interest are 
gross radiological effects (exposure-dose rate, exposure dose, decay 
rate, gamma spectra); radiochemical-composition effects (radionuclide 
partition between local and worldwide fallout, radionuclide ratios in 
particles of different size, type, and location, and relative biological 
availability); effects related to mass deposit (specific activity, mass— 
dose ratio); and various combinations of these (e.g., *'I—dose ratio). 
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The first phase of this subprogram is virtually complete. We have 
available two models for prediction (a thermodynamic model anda 
radial-power-distribitidn model) and a means of modifying each semi- 
empirically. The information flow in the system is shown schemati- 
cally in Fig. 2. 


Thermodynamic Equilibrium Model The thermodynamic equilibrium model 
developed by C. F. Miller has already been adequately described in 
detail in several publications®~' and applied with some success to the 
case of reactor excursions by C. E. Miller, Jr.® Essentially it consists 
of distributing the available radionuclides in the nuclear cloud among: 
the particles present according to the predicted equilibrium distribu- 
tion. The model takes 1400°C as the temperature below which particles 
are impenetrable to condensing nuclides. Unfortunately, adequate ther- 
modynamic input data are not available to properly use and test the 
model in several critical cases (e.g., mass chains 89, 132, and 140 
and their adjacent chains), 


Radial-power-distribution Model The radial-power-distribution model? is an 
amplification of suggestions by R. D. Cadle, R. C. Tompkins, and P. W. 
Krey, namely, that refractorily behaving nuclides distribute themselves 
among particles according to the “available volume whereas volatilely 
behaving nuclides distribute themselves according to the available sur- 
face. By assuming that a collection of spherical particles exists and 
that all.mass chains are distributed according to some power of the 
radius, one finds that the model predicts logarithmic correlations of 
radionuclide ratios for monodisperse samples. This is a happy result 
because radionuclide ratios observed in fractionated nuclear debris can 
be correlated logarithmically, at least as well as they can linearly, 
and the correlation parameters then become useful for model predic- 
tions. If one further assumes a particle distribution such that the mass 
is distributed lognormally with particle diameter with modal diameter 
xy and variance o%, one finds that all mass chains are distributed 
among the particles with variance o* and modal diameters given by 
x, exp [(b,;-1)o’]. Here b, is a slope correlation parameter. This per- 
mits a great simplification in the calculations and makes hand calcula- 
tion feasible, whereas, in the case of the thermodynamic equilibrium 
model, a computer is required. Although the thermodynamic equilibrium 
model is a more fundamental approach, the assumptions, simplifica- 
tions, and approximations involved in its use, together with the state of 
the available input data, make the two approaches, at least for the 
present, of competitive reliability. The radial-power model has a 
great advantage in being presently applicable to air-, tower-, and 
surface-burst predictions. 

The semiempirical aspect of the radial-power-distribution model 
consists of the way it utilizes empirical correlation parameters (ob- 
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tained from the radiochemical analysis of fractionated fallout samples) 
to predict fractionation effects. When so used, with a lognormal 
particle-size distribution, mass-balance requirements are automati- 
cally fulfilled for each mass chain. The assumptions of spherical 
particles and radial power distribution guarantee logarithmic correla- 
tions for monodisperse samples, and mass-balance fulfillment can be 
expected for any overall particle-size distribution. 
The thermodynamic equilibrium model can be utilized in a semi- 
empirical manner by using it as originally described to calculate the 
f 2 
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Fig. 3—Semiempirical model prediction of the effect of fractionation 
on the 1.12-hr dose rate for a burst in which solidification occurs in 
6 Séc. ‘ 


degree of fractionation for each sample and then attaching to that sam- 
ple the effect empirically expected for that degree of fractionation. It 
has not yet been demonstrated, however, that mass-balance require- 
ments will then be fulfilled. The seriousness of this defect would be 
less for dose-rate estimates than for radionuclide partition and may 
not be serious in either case, compared to the magnitude of other un- 
certainties, 

Figure 3 shows the percentage of unfractionated exposure-dosé 
rate expected for different degrees of fractionation at 1.12 hr from a 
burst in which solidification occurred in 6 sec. According to Miller’s 
estimates’ this would correspond to a yield of 25 kt. The prediction 
indicates a maximum depletion of a factor of 5 in the dose rate from. 
local fallout. We have data from the 1962 test series in Nevada which 


indicate that depletion may be much greater in‘ actual situations. The* 


observations’ require further substantiation and verification, however. 
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land-surface burst. 


In the case of volatile-chain enrichment under extreme conditions (base 
surge from an underwater burst or gas venting from a subsurface 
burst), the ratio of dose rate to the unfractionated dose rate calculated 
from the mass-95 chain may greatly exceed a factor of 10. 

Another point of interest predicted by our model is that debris 
may adhere more closely to the Way—Wigner t7?? decay rule if it is 
fractionated than if it is unfractionated; this is supported by unpub- 
lished observations on ‘fallout from weapons tests. : 

» Finally, the model has been used to predict the partition of various 
mass ‘chains among local fallout (particles of 50 or greater in diame- 
ter), intermediate fallout (particles between 25 and 50 pin diameter), 
and worldwide fallout (particles less than 25u in diameter). The pre- 
dictions are shown in Fig. 4and are compared in Table 1 with partitions 
inferred from fallout analysis. *»!° The data from the high-yield coral- 
surface burst are of poor reliability and show only qualitative agree- 
ment, The data from the low-yield silicate-surface burst agree much 
better. The agreement is somewhat fortuitous, however, because the 


10 FREILING, CROCKER, AND ADAMS 


Table 1— COMPARISON OF RADIONUCLIDE PARTITIONS PREDICTED BY. THE SEMIEMPIRICAL 
METHOD WITH OTHERS INFERRED FROM ANALYSIS OF NUCLEAR DEBRIS | 


Fractions in 24-hr © 
cloud of a high- 
yield, coral- 


Fractions:in a low-yield, 


Predicted fractions silicate-surface burst 


Radionuclide Local Intermediate* Worldwide surface burst Local Intermediate} Worldwide 
BICs 0.1 0.1 0.8 0.36 + 0.36 0.00 0.22 0.78 
89Sy 0.1 0.1 0.8 0.02 0.24 0.74 
90Sy 0.15 0.1 0.75 0.11 + 0.11 0.07 0,24 0.69 
40Ba—40L_, = 0.25 0.1 0.65 0.20 0.26 0.54 
132Te 0.25 0.1 0.65 0.18 0.26 0.56 
Zr —5Nb 0.65 0.1 0.25 0.72 0.19 >» 0,09 
98Mo (coral) 0.65 0.1 0,25 0.02 + 0.02 0.72 0.19 0.09 


*Intermediate fraction taken as 25 to 50 yw in diameter, 
fIntermediate fraction taken as 18 to about 90 yw in diameter. ~ 


particle-size distribution in the latter case differed considerably from 
that used in the model. 

In summary, the present state of the art puts us in the position of 
the meteorologist who has to predict the weather even though he can’t. 
We make the most reasonable attempt our knowledge permits and cross 
our fingers. 


Second Phase of Predicting Fractionation Effects 


We will now discuss the second, long-term phase of predicting 
fractionation effects, the phase concerned with piecewise refinement of 
weak links in the calculational chain. Figure 5 shows the fallout- 
formation processes that have to be taken into account in a fundamental 
approach, Our thinking on this is still the same as when we presented 
this diagram at the last conference.’ Perhaps the weakest links are the 
definition of particle size and the absence of accounting for agglomera- 
tion effects, and next inimportance isthe transition from an equilibrium 
(thermodynamic) to a nonequilibrium (kinetic) approach. However, the 
‘definitive work being done by Russell on the first problem’ and the 
need to properly plan our high-temperature experimental work has led 


us to attack the last problem first. Our efforts in.this-direction form . 


the subject of this section of the report. One should not jump to the 
conclusion that the thermodynamic equilibrium treatment is inapplica- 
ble or that, even if it were, thermodynamic data are no longer required. 
As will be seen, each approach has its place in the overall development. 

The plan of developing a kinetic model has been, first, to decide 
on:a mechanism; second, to assemble useful, available theoretical and 
experimental results; and, third, to integrate these and{fill in the miss- 
ing steps to obtain the complete treatment. We are now in stage two. 
We are obviously not concerned so much at'this point with developing ~ 
a new theory of fallout formation as with the application of established 
theory. : 
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Fig: 5—Fallout-formation processes in air and surface bursts. 
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Mechanism With the help of Fig. 6, the mechanism can be described in 
the simple case from the viewpoint of a condensable molecule in the 
vapor phase. To be condensed, the molecule must first diffuse through 
the gas to the surface of a particle. Upon striking the surface the mol- 
ecule can either rebound into the vapor or cling to the surface. After 


Fig. 6—Mechanism for kinetic approach to radionuclide incorporation 
by fallout particles. 


clinging to the surface awhile, the molecule may reevaporate or dif- 
fuse into the bulk of the particle. If it reevaporates, it can easily be 
deflected back to the surface by the first collision with another vapor 
molecule, Inside the particle the molecule may diffuse back to the sur- 
face. In the actual case.these processes will be occurring while the 
temperature is dropping. Simultaneous condensation of macroscopic 
quantities of evaporated carrier can occur at early stages. Particles 
can collide and: stick. Various chemical reactions and changes of spe- 
cies can occur at different steps of the way. Radioactive transforma- 
tions will proceed concurrently. Evaporation of volatile daughters of 
condensed nuclides may be important, and radiation effects may be 
considerable. A rigorous treatment would obviously be highly complex. 
However, we can visualize certain simplifications under various condi- 
tions. 


Equilibrium Approximation One such simplification is the thermodynamic 
approach already described,*’® This will apply when the times for con- 
densation and diffusion are short in comparison to the rates of radio- 
active decay and temperature decrease. This approximation will be 
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favored by small particles and rapidly diffusing molecules. Einstein’s 
equation ae 
sear 


2 


E 


Ke 2D, 


provides a helpful rule of thumb in this regard. This equation relates 
the time t required for the mean-square displacement ¢? to the coef- 
ficient of interdiffusion D,.. To estimate the time required to approach 
a significant degree of equilibrium, we can calculate the time re- 
quired for — to equal one-tenth of a particle diameter. The degree 
achieved in this time will be appreciable as shown in the section en- 
titled “Particle Diffusion.” Thus, for D,, equal to 5 x 107? em?/sec, a 
100-y-diameter particle would be well on the way to equilibrium in 
about a second. Interdiffusion coefficients for fission products in 
fused silicates are scarce. Figure 7 shows some cases that have been 
measured, together with some for elements that reasonably approxi- 
mate fission products. Thus the curves for Rb* and Cs* in Na,O+CaO- 
4SiO, would probably be similar to, but below, that for Na*, with Cst 
being lowest. Similar relations would be expected among Sr**, Ba?t, 
and Ca’?, 

To determine the validity of the thermodynamic equilibrium treat- 
ment under different conditions, one must compare the diffusion times 
with the cooling rate. The cooling rate is easily calculated by com- 
bining Hillendahl’s equations,'® 


—0.34 
T(°K) = 7000 W(kt)-?:7 an 


fmax 
and 


trmax(Sec) = 0.037 W(kt)®-48 


to eliminate t,,,,, (the time of the final maximum) and differentiating 
the result. The cooling rate, in terms of either time or temperature, is 
then given by 


-dT _ 0-10 41.34 
ae = T76W t 


The last equation was used to prepare Fig. 8, which shows cooling rate 
as a function of temperature for a wide range of total yields. We can 
now see what would happen during the 1~sec equilibration time referred 
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Fig. 7—-Some experimental values of diffusion coefficients in molten 
silicates compared with the values necessary to achieve various root- 
mean-square displacements in given times.1-14 ' 


to previously under various conditions. At a temperature of 1400°C, a 
fireball from a 10-kt device would cool about 100°C in a second, and 
equilibration would never catch up to the falling temperature. For a 
- 10-Mt burst (which is more in the high-yield range for which the treat- 
ment was originally designed), the temperature would fall only about 
10°C, and the approximation.is much more realistic, ‘if still not.en- 
tirely satisfactory. It should not be surprising; then, if the accuracy of 
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‘COOLING RATE, °K/ SEC 
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Fig. 8—Cooling rates as a function of temperature and total weapon 
yield, 7 


predictions. made with the thermodynamic equilibrium model is sensi- 
tive to yield. * 

Examination of the previous figures clarifies the nature of another 
approximation used by the present system. The equilibrium treatment 
considers: fallout formation to occur in two distinct stages. According 
to this approach, the particles during the first stage form an assembly 
of liquid drops, and condensing radionuclides distribute themselves 
among the particles in proportion to their volumes. At 1400°C the 
particles solidify, and further condensation results in surface distribu- 
tions, 

The diffusion data in Fig. 7 show no breaks, either in the neigh- 
borhood of 1400°C or anywhere else, for any of the silicates studied. 


x 


*Unfortunately, Hillendahl’s equations are based on data which do not extend 
below 2000 degrees. Below 2000 degrees radiation becomes less important as a 
means of cooling than other processes, especially the engulfment of cold air. 
The cooling rates estimated here are therefore merely lower limits. Storebd*! 
has attempted to accountfor the engulfment of cold air in air bursts, His cooling 
equations lead to cooling rates of the order of 103 degrees per second at 1400°C. 
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Fig. 11—Relative disintegration rates as a function of time for indi- 
vidual fission products inthe case of **8U fission produced by a thermo- 
nuclear neutron spectrum. 
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Fig. 12—Gamima-photon spectra from fission products at various 
times after *8U fission produced by athermonuclear neutron spectrum. 


The relative contributions of the important gamma-emitting nu- 
clides to the total dose-rate have been plotted vs. time in Fig. 14. It is 
interesting to note that from a few days onward much of the dose rate 
at any given time can be accounted for in terms of two or three nu- 
clides and that volatilely behaving mass chains contribute prominently 
to the dose rate. The dose-rate contributions given here pertain, of 
course, to unfractionated debris. Since many of the hard emitters (e.g., 
cesium and iodine isotopes) are subject to fractionation, the relative 
dose-rate contributions will be greatly modified in the case of frac- 
tionated debris. 
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Fig. 13—Gross exposure-dose-rate decay for products of **8U fission 
broduced by a thermonuclear neutron spectrum. 
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Fig. 14—Relative contributions tothe exposure-dose rate as a function 
of time for individual fission products in the case of ?%8U fission pro- 
duced by a thermonuclear neutron spectrum. 
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APPENDIX 
TERMS AND DEFINITIONS 


RADIONUCLIDE FRACTIONATION Any alteration of radionuclide com- 
position occurring between the time of detonation and the time of 
radiochemical analysis which causes the debris sample to be non- 
representative of the detonation products taken as a whole. 

FRACTIONATED Two substances, either radionuclides or inert ma- 
terial, are fractionated with respect to each other in a sample of 
debris if they are not present in their representative ratio. The 
term is meaningless when applied to a single substance. 

*POTENTIAL FRACTIONATION The existence of different composi- 
tions in various portions of a sample of debris when these portions 
are subject to separation by subsequent natural processes, e.g., 
the existence of particles of different size with different composi- 
tions in the same portion of a nuclear cloud or of a slurry in which 
different nuclides are distributed Heteen liquid and solid phases 
in different proportions. 

*NATURAL FRACTIONATION Fractionation “anadieea by the pro- 
cesses occurring subsequent to detonation. 

*PRIMARY FRACTIONATION IN AIR, TOWER, ANDSURFACE BURSTS 
Actual fractionation produced by meteorological, gravitational, and 


centrifugal separation of potentially fractionated particles in a- 


cloud. 

*SECONDARY FRACTIONATION IN AIR, TOWER, AND SURFACE 
BURSTS Fractionation produced by interaction of the debris with 
radioactively inert environment (e.g., leaching with water or pref- 

- erential small particle adsorption on surfaces). 

*ARTIFACTITIOUS FRACTIONATION Fractionation resulting from 
human intervention (e.g., biased sampling, incomplete decontami- 
nation of collector surfaces, and faulty analytical techniques). 

FRACTION OF CONSITITUENT OR EXTENSIVE PROPERTY IN A 
SAMPLE The ratio of the quantity of said constituent or property 
to the total amount produced by the device and measured at the 
same time whennecessary (replaces the “bomb-fraction” concept). 


* Asterisks indicate convenient, but not essential, terms, 
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EQUIVALENT FISSIONS The number of equivalent fissions present in 
a sample as determined by analysis for any constituent is the prod- 
uct of the total number of fissions produced by the device and the 
fraction of the total amount of constituent present (replaces the 
“fissions in the sample” concept). 

FRACTIONATION COEFFICIENT OF TWO CONSTITUENTS IN A SAM- 
PLE The ratio of their fractions or equivalent fissions. . 

*THE FRACTIONATION RATIO The ratio of fractions or equivalent 
fissions of Sr to ®Zr. 

*DEGREE OF FRACTIONATION The Briggsian logarithm of the frac- 
tionation ratio. 

*TERACALORIE 10'* calories (the kiloton equivalent of energy). 

*TERACALORIE EQUIVALENT OF ANUCLIDE The ratio of quantity 
of nuclide produced in a detonation to the energy release in ter-. 
acalories by fission. 

TERACALORIE EQUIVALENT PER SQUARE MILE A unit of con- 
tamination surface density (replaces the “kilotons per square mile” 
concept). 
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SMALL BOY SHOT FALLOUT 
RESEARCH PROGRAM 


CARL F. MILLER and JAMES D, SARTOR 
Stanford Research Institute, Menlo Park, California 


ABSTRACT 

In this paper the general outline of a field test program is presented 
along* with the application of the outline to the design of the experi- 
mental program for the Small Boy shot. Also, a preliminary analysis 
and correlation of the data obtained on two of the projects is presented 
in a discussion and an evaluation of the intensity —activity ratio and the 
intensity—area integral per unit fission yield (i.e., the ratio of roent- 
gens per hour at 1 hr per kiloton per square mile) ey the fallout from 
the Small Boy: shot. 


INTRODUCTION 


Over the past decade many significant advances indevelopment and 
field testing of nuclear explosives have been made. In the field testing 
of nuclear weapons, explosion products, called fallout, received in- 
creasing attention as the testing proceeded. One of the major reasons 
for the increased attention to fallout was that, during this same period 
of time, better understanding and recognition of the radiological 
hazards to biological species from exposures to nuclear radiation 
were accomplished. In addition, studies of the application of nuclear 
explosives in possible future warfare showed that exposure of human 
beings and other living species to the nuclear radiations associated 
with fallout could result in many fatalities. 

Field experimental programs designed to obtain information on 
fallout became increasingly complex as time progressed. The evolve- 
ment of experimental procedures in the field tended to reflect, at the 
time, both the understanding of the radiological hazards and the 
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technical parameters that influenced the character and magnitude of 
the hazards. 

The fallout study program for the Small Boy shot, as part of the 
last series of aboveground field tests of nuclear weapons, was co- 
operatively organized and sponsored by the Office of Civil Defense 
(OCD), the Department of Defense (DOD), and the Atomic Energy 
Commission (AEC), The projects sponsored by OCD were under the 
operational control of the Defense Atomic Support Agency (DASA); 
those sponsored by AEC were under the operational control of their 
Civil Effects Test Organization (CETO). 


WEAPON-TEST-FALLOUT STUDY PROGRAM 


A general outline of the information and data requirements and 
the program objectives was prepared prior to the establishment of the 
experimental projects for the Small Boy shot field test. This outline 
was reviewed by various representatives of OCD, DASA, and AEC. 

The outline, reflecting operational data requirements by the three 
agencies, included an overall design of field experimental programs 
for improving the data base on fallout, The data needs were divided 
into the following four subprograms or study areas: 


1. Fallout formation 

2. Distribution of fallout 

3. Definition of radiological exposure environments 

4. Alteration of exposure environments by radiological counter- 
measures 


For each of the four study areas, a general list of the parameters 
for direct or indirect experimental evaluation was prepared. These 
are given in Tables 1 to 4. The first three study areas were organized 
to provide input data for fallout-model development and radiological- 
hazard evaluations needed by all three agencies. The fourth study area 
was designed somewhat more specifically for providing information for 
use in the evaluation of radiological-safety programs for AEC and in 
the design and evaluation of radiological-defense systems for OCD, : 

The main overall objective of the weapon-test-fallout study was 
to make experimental determinations and to obtain and report experi- 
mental data needed for describing and defining the parameters in each 
of the four study areas. The objective for the field-test program was 
limited to the taking and reporting of experimental data because each 
agency sponsored other research tasks for the correlation and evalua-~ 
tion of the observed data and because inclusion of data analysis in 
field-test reports previously had usually resulted in extended delays in 
report publication. 


46 Table 1— FALLOUT-FORMATION PARAMETERS 


Initial condition variables 
Composition of the device 
Composition and physical properties of the environmental materials at shot 
point 


Geometric parameters 
Dimensions and location of the visible fireball and cloud as a function of time 
after detonation 
Dimensions of the apparent crater and crater lip 
Spatial distribution of radioactivity in the fireball and cloud volume as a func- 
tion of time after detonation 
Spatial geometry of toroidal circulation of the rising fireball and cloud 


Energy balance parameters 

Air overpressures and ground-shock energies 

Radiant energy emitted as a function of time after detonation 

Temperature of fireball gas as a function of time after detonation 

Material composition of fireball and influx, mixing, and ejection of gases and 
particles as a function of particle size and time after detonation 

Energetics (dynamics) of particle circulations in and about the rising fireball 
and cloud 


Related resultant parameters 

Chemical composition of fallout particles 

Radiochemical composition of particles as a function of particle diameter, in- 
cluding specific activity of the particles 

Radioactivity and mass-distribution functions 

Solubility of radioelements as a function of particle size 

Shape, structure, and density of fallout particles and internal distributions of 
crystalline and glass phases and of the radioelements within the fallout 
particles 7 


Table 2— FALLOUT-DISTRIBUTION PARAMETERS 


Initial-condition variables 
Spatial distribution of radioactivity and mass of particles at the time of cloud 
stabilization and at the times when particles start free-gravity fall 


Particle-fall-rate parameters 
Density distribution of fallout particles as a function of particle diameter, 
shape, and number or mass 
Atmospheric conditions-(air density, viscosity, temperature, pressure, rela- 
tive humidity, etc.) as a function of altitude over the fallout area 
Particle shape and drag factors 
Agglomeration effects (including incorporation of solid particles in raindrops) 


Meteorological parameters 
Wind speeds and direction as a function of altitude and time over the whole 
fallout area 
Vertical air motions 
Effect of terrain features and diurnal heating and cooling on the air-flow 
patterns 
Deposition dynamics 
Time of arrival of fallout particles at given locations on the ground or in the 
air as a function of particle size 
Duration of fallout deposition at given locations in the fallout area 


(Table 2 continues on page 47) 
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Rate of accumulation of particles in terms of radioactivity and mass at given 
locations in the fallout area 
Effect of rainfall on deposition rates and air concentration of particles 


Area distribution of the fallout particles 

Standard intensity contours (usual fallout pattern expressed in roentgens per 
hour at 1 hr at 3 ft above an open contaminated area) 

Intensity— activity ratios (in roentgens per hour at 1 hr per kiloton per square 
mile at various locations in the fallout area 

Intensity— area integral and fraction of the device in the fallout area 

Fraction of the total radioactivity and mass of particles deposited on the 
earth’s surface as a function of time after detonation 

Surface-density contours (mass of particles deposited per unit area) 

Gross fractionation numbers for the fallout deposited at given locations in the 
fallout area as a function of particle diameter 


Table 3— RADIOLOGICAL-~EXPOSURE- ENVIRONMENT PARAMETERS © 


Initial contamination of surfaces 

Surface wind speeds and directions and relative humidity 

Effect of large- and small-scale surface roughness features of terrain on the 
airflow over the surface 

Shape, size, general orientation, and character of the surface 

Impaction and retention efficiency of particles and surfaces as a function of 
the above parameters, particle diameter, and surface density of the deposit 

Initial contamination factors for plants, animals, and humans 

Inhalation of fallout particles during fallout deposition by animals and humans 
as a function of particle diameter, wind speed, and related factors 


Natural redistribution of deposited fallout particles and radiation sources 

Transfer of soluble radionuclides from the fallout particles to surfaces by 
dew and rain as a function of time, particle size. and amount of dew and/or 
rain 

Effect of erosion by wind and rain of the deposited fallout particles from and 
onto surfaces on the ionization rate and particle surface density as a func- 
tion of particle diameter, type of surface, wind speed, amount of rainfall, 
time period over which erosion occurs, and general terrain features 

Movement of soluble radionuclides and particles in and over soils due to rain- 
fall as a function of particle diameter, gross deposit surface density, soil 
type, number of rain showers, and amount of rainfall 

Movement of soluble radionuclides and particles in streams and rivers as a 
function of particle size, deposit surface density, water flow rates, and 
other stream or river characteristics 

Variation of retention efficiencies and contamination factors on meteorological 

_ parameters (wind, rainfall, relative humidity. etc.) 

Inhalation of redistributed (wind-blown) particles 

Movement of particles by people and vehicular traffic 


External gamma-radiation-exposure parameters 

Variation of air ionization rates with time after detonation at given locations 
in the fallout region (reference geometry: 3 ft above an open contaminated 
area) ; 

Gross ionization-rate decay of the radioactivity carried by particles as a 
function of particle diameter 

Attenuation of gamma rays by terrain-roughness features 

Variation of the air ionization rate with altitude over given locations in the 
fallout area 


(Table 3 continues on page 48) 
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Response factors for radiation detecting instruments: 

Area radiation-contribution factors for discontinuous source geometries and 
variable source intensities 

Skyshine-contribution factors as a function of surface roughness and source 
geometry 


Contact radiation exposures 
Beta dose and dose rates in air over contaminated surfaces as a function of 
distance from the surface 
Absorbed beta dose to the tissue of contaminated plants, animals, and humans 
Dependence of absorbed contact beta doses on meteorological factors, ex-~ 
posure times, and environmental factors (all biological species, including 
insects) 


Internal radiation-exposure parameters 

Biological availability of individual radionuclides as a functionof particle size 

Foliar absorption of radionuclides by edible plants or plant parts as a function 
of particle diameter 

Uptake rate of radionuclides by edible plants grown on tilled and untilled soils 
as a function of particle diameter or location, plant growth rates, soil type, 
and climatic variables : 

Ingestion and assimilation rates of radionuclides by tissues of animals and 
radionuclide concentrations in animal tissues and food products from ani- 
mals fed with contaminated plant parts and water obtained from different 
locations in the fallout area 

Absorbed doses to plants, animals, and humans from ingestion of radio- 
nuclides (all intake paths) 


Table 4—- ALTERATION OF EXPOSURE ENVIRONMENTS BY © 
RADIOLOGICAL COUNTERMEASURES 


Attenuation of gamma-ray intensities 

Shelter protection factors as a function of wall materials and thickness, 
radiation-source geometry, and compartment size, shape, and location in 
building 

Attenuation by barriers as a function of barrier height and thickness and other 
geometric parameters, such as the distance from barrier or the dimensions 
of area enclosed by the barrier ; 

Attenuation by burial of the fallout: particles by covering them with soil or by 
plowing (i.e., mixing the fallout with soil to various depths) land areas 


Decontamination 

Effectiveness of the removal of fallout particles from contaminated surfaces 
as a function of the method and its procedural parameters, particle diam- 
eter of the fallout, surface density of the fallout, weathering, type of sur- 
face, and effort applied 

Disposal of fallout particles removed by decontamination 

Exposure dose to decontamination crews during decontamination Gperatious as 
a function of the type of area that is cleaned, the operation schedule, and the 
methods of decontamination 

Effectiveness of water purification methods in the removal of soluble radio- 
nuclides and fallout particles from exposed water sources 

Effectiveness of food processing and treatment methods in the removal of 
soluble radionuclides from contaminated. food products (i.e., milk ‘decon- 
tamination, removal of particles from grains in threshing and milling, 
desorption of radionuclides from. vegetables into water, etc.) 


(Table 4 continues on page 49) 
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Exposure control methods 

Reduction in exposure doses by scheduling of exposures for operations 

Reduction in exposure doses by adjustment of living routines as a function of 
time after detonation and radiation intensity 

Reduction in absorbed doses to animals and humans by storage and allocation 
of contaminated foodstuffs 

Ingress of fallout particles into shelter spaces through openings and ventila- 
tion systems and effect on the exposure dose or on the shelter protection 
factor . 

Verification of shielding provided by specially designed shelter~entryway 
configurations 

Proof tests of operational radiological-defense systems and overall effective- 
ness evaluations 


The individual project titles, as parts of the field-test fallout 
program for the Small Boy shot, are listed in Table 5. The projects, 
with a few exceptions, were very successful in meeting their individual 
objectives. This general success was due, first, to the detailed de- 
velopment of the project experimental designs on the part of the 
project leaders and the efforts of all the project personnel in carrying 
out the work and, second, to the favorable winds that carried the fallout 
over the established station array. It should be noted that projects 
involving the design and testing of sampling equipment for possible use 
in future programs were included in the program. 


Table 5— PROJECT PARTICIPATION IN SMALL BOY SHOT 
FALLOUT PROGRAM 


. Thermal measurements and fireball photography 

. Meteorological measurements 

. Radiological survey of the fallout area 

Aerial survey of the fallout area 

Fallout collection and gross sample analysis (onsite and offsite) 
Field ionization-rate measurements 

Long-range fallout collection and radiochemical analysis 
Physiochemical and radiochemical analysis of fallout samples 

. Terrain shielding measurements and field spectra (onsite and offsite) 
10. Contamination of plants 

11. Ingestion of radionuclides in fallout by native animals 

12. Assimilation of radionuclides in fallout by humans 

13: Contamination of milk (unscheduled) 

14. Proof-test of a prototype rocket-mounted collector 

15. Test of a prototype fallout sampler 

16. Comparison of aerial-survey instruments 


woaonooakhwnde 


THEORY OF THE INTENSITY-ACTIVITY RATIO 
AND THE INTENSITY-AREA INTEGRAL 
Known relations between the observed radiation rates and the 


radioactivity carried by fallout particles spread over real (open) 
terrain are required to evaluate radioactive-material balances and 
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the relative worldwide disposition of all the radionuclides produced in 
a nuclear detonation. The ratio of the roentgens per hour at 1 hr to 
kilotons per square mile, where the kilotons are taken as being a 
measure of the radioactivity carried by the fallout particles, is used 
in the integration of fallout-pattern contours to estimate the fraction of 
the device accounted for in the pattern. Other more detailed uses of 
the ratio include the effects of terrain roughness and radionuclide 
fractionation on gamma-radiation intensities, Some mathematical fall- 
out models utilize this ratio in estimating standard intensities (i.e., 
the roentgens per hour at 1 hr values at 3 ft above a uniformly coh: 
taminated plane) in computing fallout patterns. 

It is generally known that the intensity—activity ratio .can be 
defined .for the intensity as observed at a given location or, as an 
average, for a whole fallout area, where it is sometimes called the 
intensity—area integral per unit fission yield. It is convenient in data 
analyses to express the surface density of the radioactivity in fissions 
per square foot rather than in kilotons per square mile. Both repre- 
sentations of the surface density of radioactivity are independent of 
time after detonation. Since the relation between the number of fissions 
and the energy released in fission:is about the same for most common 
fissile materials, namely, (1.45 + 0.03) x 1073 fissions per kiloton of 
fission yield, the two PepR eR eutanigns of the surface density of radio- 
activity are related by . 


Ap = 5.20% 10° A, ; (1) 


where A; is in fissions per square foot and Ay is in kilotons per 
square mile. a 

The intensity — activity ratio for a given location in a fallout region 
is defined by : 


Ky =3- , — (2) 


where I, is an observed or a measured value of the (standard) intensity 
in roentgens per hour at 1 hr at 3 ft above an extended open flat area 
uniformly covered with fallout particles carrying the appropriate 
amount of radionuclides to result in the surface density, A;, in fissions 
per square foot or the equivalent surface density, A,,; in kilotons per 
square mile, as given by Eq. 1. The standard intensity, I,, is usually 
determined. from a measurement of the roentgens per hour at the 
location after all the fallout has been deposited and a decay correction: 
of the observed intensity Be the ‘standard time of 1 hr after detonation 
has been made. - 
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The value of K, has been computed for uniform depositions of un- 
fractionated fission products on ideal smooth planes. However, the 
calculated values of K, from observed values of I, over real terrain 
for a given value of A, or A, are lower than those computed for the 
ideal plane because of attenuation of the gamma-ray intensity by both 
small- and large-scale roughness of real terrain, nonidealized re- 
sponses of radiac instruments used in measuring the intensity, and 
decreased contributions from the more volatile fission-product radio- 
nuclides whose relative concentrations in the fallout are decreased 
(i.e., the radioactive composition is altered or fractionated). The 
values of I,, relative to those for the fission products, are increased 
by the production of induced radionuclides (usually through neutron 
capture) in both weapon components and nearby environmental mate- 
rials. The representation of K,, including consideration of the four 
factors discussed in this paragraph, is given by 


Ky = Da(r,k,, + r,k,) (3) 


where D= instrument response factor for the assumed detector— 
radiation source geometry 
q = terrain attenuation factor 
kK, = computed value ‘of I,/A; for 3 ft above an ideal plane uni- 
formly contaminated with unfractionated fission products 
r= gross fractionation number and: is equal to the ratio of I, 
for the fractionated mixture of radionuclides to the I, for the 
unfractionated mixture . . 
k;= computed value of the I,/A; contributions from induced 
radionuclides for the same detector—source geometry as 
_ for ky se 
r, = gross fractionation number for the induced radionuclides 
and is related to k; in the same way that r,, is related to k,, 


If it is accepted that the radionuclide composition varies with the 
size of the fallout particles, the major factors in Eq. 3 that are de- 
pendent on particle diameter are ry, and r;. In more precise treat- 
ments, D and q-also depend on. the radionuclide composition because 
-both of these parameters are functions of gamma-ray energy, which, in 
turn, depends on the relative abundance of the radionuclides in the 
fallout. : : : 

The value of the intensity—area integral for a fallout pattern usu- 
ally is determined by integration of the standard intensity contours 
from the highest intensity to a selected low-valued contour. The frac- 
tions.of the radionuclides deposited outside the lowest contour are not 
included in the pattern summations, Also, the quantities that are buried 


a a eh a 
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in the crater and crater lip would not be included in the pattern sum- 
mations. The intensity—area integral is defined mathematically by 


Jp= fi, da | (4) 
R 

= f KA, da (5) 
R 


where da is the incremental area, dx dy, and R is the region within the 
lowest I, contour included in the integral. If K) is designated as the 
weighted average value of the K,’s and C as the ratio of intensity area 
for the region R to the integral for the region enclosed by. the contour 
of I, equal to zero, then Jp, can be represented by 


J, = CK)BW . (6) 


where W is the total weapon yield and B is the fraction of the total 
yield due to fission. If the Ky value for the unfractionated mistune of 
radionuclides produced in the detonation is designated as K} on = 1, 
r;= 1), then the fraction of the device deposited in the region, R, is 
given by 


igi ae (7) 
D Ko 
or 
P= — (8) 
Kp 


If Ky is written in the form of Eq. 3, then Eq. 8 becomes 


a C[F,) ky, + Fiki | 
eee ee) 


(9) 


In Eq. 9 the two bracketed terms refer to the intensities over an ideal 
smooth plane on the assumption that the product Dq for Xe is equal to 
Dq for Kj}. — 

The values of CK, and Fp for a fallout pattern can tba estimated 
from the intensity—area integral and the values of Kens k,, D, and q or 
q. The value of C can be estimated if K, is evaluated separately, and Ky 
can be evaluated if the: variation of Ky with particle size and the frac- 
tion of the total radioactivity carried by particles of different sizes are 
- known. With these two types of data, Ky is calculated from 
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Ky = 20 £,Ko; (10) 
j 


where f, is the fraction of the total activity carried by particles with a 
range in diameter designated. by j and Kg, is the value of Ky for the 
same particle-size groups. However, for a determination of the varia- 
tion of Ky with particle size, the variation of r,, and yr, with particle 
size must be known, and approximate values of q or Dq must be 
determined. 

The average value of Vip (and of r,) for the intensity summation 
can be calculated in the same way as K, so that 


\ 


Fp = Un, (11) 


where r, is the Te, value for the jz particle group. 


TREATMENT OF THE DATA 


Preliminary analyses and correlations of the radiochemical data 
and gross-sample-activity, sieved-sample-activity, decay-rate, and 
intensity measurements in the field on the fallout from the Small Boy 
shot have been made, leading to the evaluation of preceding param- 
eters. Certain simplifying assumptions are made for portions of the 
analyses to facilitate the treatment of the data and the presentation of 
the computational results. 

First, it is assumed that, in Eq. 3, the values of r;k, are small 
compared with the values of r,,, and Ky is represented only by 


Ky = Dar,,K . (12) 


The values of Ky) were determined according to Eq. 2, in which the 
I, value of each fallout-collecting station was derived from correla- 
tions between data from gamma intensity time recorders (GITR’s) and 
_activity measurements on the collected samples.'” The values of A; 
were derived from radiochemical analyses of the collected fallout 
samples, and the fission yield equivalent for each sample was based on 
the amount of *Zr found in each sample.’ , 

“The values of r;, were determined by taking the ratio of the ion 
current per fission of a gross fallout sample, as measured in a cali- 
brated 47 ion chamber‘ to the computed ion current per fission for an 
unfractionated fission-product mixture.’ In the field project,’ the 
radioactive content of all collected gross samples and all sieved 
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fractions of gross samples was measured either with the ion chamber 
or with a calibrated gamma Scintillation counter. The activity of many 
samples was measured with both instruments. The measured ion cur- 
rents in milliamperes and the photon pulse rates in counts per minute 
for each sample were decay~-corrected to 100 hr after detonation. 

If the measured ion current per *Zr fission at 100 hr from the 
ion chamber is designated as i(100) and the calculated ion current per 
fission for the unfractionated fission product at 100 hr after fission is 
designated as i,(100), then the gross fractionation number for the 
measured sample at 100 hr is 


i(100) 


r,, (100) = ;,(100) 


(13) 


The definition of the r,, by ratios of the ion currents from the Pa 
chamber implies the assumption that i/iy is equal to I Js where a is 
the roentgens per hour for an unfractionated fission-product mire: 
This assumption is fairly valid for most measurements on fission- 
product mixtures.’ For correlation purposes, a fractionation number, 
Tio is defined for the calibrated gamma-scintillation-counter mea- 
surements as 


i’(100) 
son) 


rf, (100) = (14) 
where i’(100) is the observed. counting rate per *7r fission at 100 hr 
after detonation for a sample of fallout and ig(100) is the calculated 
counting rate per fission for the unfractionated fission-product mixture. 
The r,, of Eq. 12 is for the i ratios at 1 hr after detonation rather 
than at 100 hr. For evaluation of r,, the i(100) values were decay- 
corrected to 1 hr after detonation (or after fission), The assumed 
decay-correction factors for i and i’ are represented by 


i(1) = 250 i(100) (15) 


and 


\ 


i(1)=200i(100) ¢ eae (16) 


The. only available computed values of i) and ig for the 47 ion 
chamber ang the calibrated gamma scintillation counter, respectively, 
are. those’'® for the products from thermal-neutron fission of 7°U, For 
this mixture i)(100) is 2.48 x 10° ma/fission and ij(100) is 2.10 x 107° 
counts/min per fission. The two respective calculated decay-correction 
factors from 100 to 1 hr after fission are represented by _ 
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ig(1) = 435 i,(100) (17) 
and 
ig(1) = 246 i5(100) (18) 
Thus , \ 
_ i(100) 
ry, = 0.575 (100) (19) 
and 
7a i(100) 
Yip = 0.813 ;,(100) (20) 


Actually the decay factors for i and i’ should vary with r,, or with 
the particle size of the fallout. Thus, for samples in which the degree 
of fractionation (u products) is small, the decay factor from 100 
to 1 hr should approach 435. In this preliminary analysis, the previ- 
ously assumed decay factors (Eqs. 15 to 18) were used for all samples; 
neither observed nor specially computed decay factors for all fallout 
samples were available for application in Eq. 19. The assumption that 
all samples had the same decay factor could result in an error ofa 
factor of 2 in the r,, value for some of the fallout samples. 

' The variation of i(100), r;)(100), and rj with particle size was in- 


‘vestigated by using ion-chamber and radiochemical-analysis data’ on 


sieved fractions of collected.fallout samples. Two sets of the ion- 
chamber measurements were reported. In addition to the field mea- 
surements, which were corrected to absolute ion-current values, 
an uncorrected second set of measurements was reported by the 
radiochemical-analysis project.° These data were corrected to 100 hr 
after fission by using (t/100)'-" for computing the decay factors. The 
decay-corrected data were then’ correlated with the field measure- 
ments to obtain a correction factor for conversion of thé data to 
absolute units. This was done by computing the ratio of the two mea- 
surements (corrected to 100 hr) on each sample. The median value of 
this correction factor was calculated to be 0.788, which is close to the 
value of 0.78 obtained in the past with a radium standard.’ After appli- 
cation of this second correction factor to the second set of measure- 
ments, the geometric mean value of i(100) was calculated from the two 
sets of the ion-chamber measurements and the reported® fission-yield 
content of each sample. The i(100) values thus obtained are given in 
Table 6 and are’ plotted in Fig. 1 as a function of the midrange particle 
diameter, d, for each sieve fraction. Subsieve analyses indicated the 
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Table 6—ION CURRENT PER FISSION FOR SIEVED FRACTIONS OF 
FALLOUT SAMPLES 


[i(100) (1072° ma at 100 hr per fission)] 


Midrange particle diameter, p 


Station 22 66 132 264 530 1060 2120 >4000 
100 3.88 3.30 0.950 1.02 0.676 0.844 1.32 
101 3.74 1.48 1.14 0.855 0.680 0.791 1.48 
201 4,50 2.07 1,02 0.668 0.967 0.751 0.846 15,.9* 
203 0.271* 1.70 1.72 0.678 0.997 0.662 1.60 
305 3.96 1.78 1.10 1.04 0.983 0.784 0.741 1.08 
403 3.90 2.34 0.782 1.02 0.693 1.04 4.62* 3.24* 
505 4,42 2.39 0.669 1.00 0.768 16.8* 8.44* - 4,28* 
507 0.821* 1,37 1.01 0.656 0.726 19.7* 
693 3.52 1et2 0.610 0.462 0.682 35.9* 0.880 
707 5.49 1.86 1.06 1.00 0.615 10.3* 
813 3.56 0.942 1.25 0.810 3,.90* 1.53 1.52 
Mean 4.06 1.82 0.989 0,848 0.759 0.871 1.05 1.26 
Smoothed 

curve 4.06 1.82 1.10 0.78 0.76 0.87 1.05 1.26 
¥ ,, (100) 1.64 0.734 0.444 0.815 0.306 0.351 0.424 0.508 
Y fy 0.942 0.422 0.255 0.181 0.176 0.202 0.244 0.292 


* Values not used in calculating mean values. 
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Fig. 1—Variation of i(100 with particle diameter. 
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presence of particles with diameters as small as -1to3vuw in all 
samples; therefore the midrange diameter of the 0- to 44-p fraction 
was taken as 22 yu. 

' The mean values of i(100) were computed for each particle-size 
fraction on the assumption that the average radionuclide composition 
carried by particles of a given size is the same everywhere and is 
independent of where the particles land. After correlation the two sets 

' of ion-chamber data were generally in very good agreement; therefore 
most of the spread in the i(100) values for each particle size is due 
mainly to spread in the radiochemical results. The spread in the values 
is largest for the smallest and largest particle diameters where the 
analytical sample sizes were Smallest. In the midrange of the curve, 
the minimum-to-maximum spread in i(100) is a factor of about 3. 

The major significance of the data is in the variation of i(100) 
and therefore of r,(100) and r,;, with particle size. From the average 
particle diameter of 22 1 up to an average particle diameter of about 
400 pt, i(100) decreases by a factor of about 5.5. At 400 » diameter, a 
minimum in the curve exists. A similar variation with particle size in 
the gross solubility of the radionuclides is observed.? 

The Ky, values for the collecting-station locations in the fallout 
area from the Small Boy shot were calculated with adjusted and 
averaged A, values; correlations of the variation in i’(100) and i(100) 
with the median particle size of the fallout particles at each station 
were used to make the adjustments in A,;. Initial correlations of i’(100) 
indicated a consistent difference in value depending on the time at 
which the samples were allocated for chemical analysis; the median 
value of the correlation factor for this difference was calculated to be 
1.56. Since evidence for locating the cause of this difference is 
_presently not available and since the calibrated gamma-scintillation- 
counter measurements were a highly unlikely source of error, the 
reported A, values® for the samples with the higher values of i’(100) 
were increased by a factor of 1.25, and the A; values for the remainder 

~ of the samples were decreased by a factor of 1.25. If the same type of 
correlation, including comparisons with the observed I, values, are 
applied to a third set of samples collected at offsite locations, the 
reported A,r values’ are found to be more or less consistently high by a 
factor of 10. ; 

The data for the fallout samples leading to the calculation of Ky for 
each location are summarized in Table 7. Included are the activity 
median particle size, ds); the ion-chamber measurements on the gross 
samples, 1,, I,, and I; the values of i,(100), i,(100), and i’(100), as 

calculated from the initially adjusted values of A; (given as Aj); the 
smoothed-curve values of i(100) and i’(100), as taken from the curves 
in Figs. 2 and 3; the recomputed values of A;, where A, is I,/i(100), 
A; is I,/i(100), and A, is 1’/i’(100); the average value of A,; I,; and Ky. 
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Table 7—-SUMMARY OF CALCULATIONS* OF A; AND Ky FOR THE SMALL BOY SHOT FALLOUT 


Station— 1’ (160), i,(100), —i,(1.00), i’(100), i (100), ¥ (100), 
sample 1,(100), I,(100),  counts/min 107% ma 107% ma 107? counts/min 107% ma 1079 counts/min 
No, deo ma/ft? ma/ft? per sqft fission fission ~ fission fission fission 
100-1 640 212x108 =220x10-8 ~=—-2,22x 105 0.763 0.791 0.802 0.836 0.785 
100-2 640 220x10-8 220x1078 = 2.37108 1.048 1.048- 1,130 0.836 0.785 
101-1 1100 174x107? ~—-1,30 x 108 1.193 0,890 0.880 0.893 
101-6 1100 174x107 ~— 1,29 x 108 0.650 . 0.482 0.880 6.893 
103-9 131x 1071! 0.0570f (1.0) (0.93) 
200-9 130. -122x1079 «1143x1078 ~=—- 1.16 x 104 1.100 1.288 1,046 1.13, 0.985 
201-7 750 831x10-§ 7.95 « 108 0.950 0.908 0.842 - 0.835 
201-10 750  113x107' 831x10-8 = 7,92« 108 0.876 0.645 0,614 0,842 0.835 
202-9 827 x 10-12 0.484 (0.843) (0.84) 
203-16 780 246x1078 258x1078 0.779 0.816 0.844 0.840 
203-21 780 266x10-8 258x1078 2,65x105 0.873 0.846 0.870 0.844 0.840 
204-9 “ 300x107! 0.848 (0.85) (0.85) 
207-9 45110712 0.780 (0.90) (0.91) 
209-9 488 x 10712 9,175F (0.95) " (0,92) 
300-9 259x 10710 2.53103 1,006 ; 0.984 (0.85) (0.74) 
303-9 530 746x1078 673x1078§ 8,34x1i0® 10.10 9.12T 11.28f 0.830 0.758 
305-2 660 460«1078 = 3,34 105 0.788 0.572 0,837 0.800 
305-4 660 . 452x108 4600x1078 0.820 0.834 0.837 0.800 
305-9 660 456x107° 460x109 5.31x10 0.818 0.825 0,954 0,837 0.800 
309-16 167x107" 4,75+ i (0.85) (0.85) 
311-9 3841x1071? 0.573 (0.90) (0.91) 
400-9 469 x 10710 * 3.41x 103 0.635 0.462 (0.84) 0.74 
401-9 380 473x10~§ 4456x1078 = 5.48 108 0.663 0.639 0.768 0.852 0.744 
403-2 340 629x 1078 = 3.45 x 105 0.955 0.524 0.856 0.750 
403-3 340 629x 10-8 = 5,43 x 108 1.182 1.022 0,856. 0.750 
403-4 340 -697x1078 6 29x 1078 0.817 0.737 0.856 0.750 
405-9 310 192x108 §=214x1078 ~—1,23 x 105 0.725 0.808 0.465 0.868 0.760 
407-9 310 158x107! 159x107! = 1.20x 103 0.963 0.969 0.731 0.868 0.760 
413-9 596 x 107!2 0.390T (0.95) (0.84) 
501-9 250 406x107! 893x107! 4,70x 103 0.708 1.556 0.820 0.896 0.792 
502-9 379 x 1071? 0.907 (0.88) - (0.78) 
503-9 310 = 1821078 =182x1078 = 1,59. x 105 0.847 0.847 0.740 0,868 0.760 
505-1 260 - 9334x1078 — 2.99 108 0.994 0.890 0.890 0.780 
505-6 260 345x10-§ 3341078 3.42 x 105 0.876 0.849 0.869 0.890 0.780 
507-4 180 = 890x107 1241078 = 1.53 « 108 0.940 1.310 1.616 * 0.990 0.870 
507-5 180 =-128x 10-8 =124x1078 ~=—-1.53x 108 1,202 1.164 1.437 0.990 0.870 
509-9 286x 10-9 650x107 = 2,98x 104 1.173 2.665f 1,222 (1.0) (0.88) 
513-9 296x107? 0.304} (1.1) (0.95) 
601-9 104x107? 1881078 == 1.07x 104 0.971 1,758 1.000 (1.1) (0.95) 
603-1 150 9169x1078 =173x10-8 = 2.04 x 108 0.478t 0.489 0.577 1,06 0.930 
603-4 150 173«1078 = 2.01 x 108 0.112 0.130 1.06 0.930 
605-9 160) 9-177 1078 144K 10-8 = 1,64x105 =. 1,106 0.900 1,025 1.03 0,905 
607-9 128x 1078 1,16 x 104 0.800 0.725 (1.0) (0.88) 
700-9 145 x 107! 0.900 (1,4) (1,2) 
701-9 176 107! 0,682 (1,3) (1,1) 
703-9 1711079 1.75 x 104 0.768 0.785 (1,2) (1.0) - 
704-9 1311071! 1,928 (1.1) (0.95) 
707-3 145 683x10-® =767x 1079 6,23 104 0.905 1.016 0,826 1.07 0,942 
107-6 145 767x10-§ = 7,44 104 1.056 1,023 1.07 0.942 
813-10 90 914x107 813x107! 9,61 103 1.057 0.910 , 1.112 1.42 1,17 
814-4 88 850x107!9 gs6x10- 7.88% 103 0.265¢ 0.276t 0.246 1.43 1,18 
815-10 37 133x107 8.86 x 103 3.03 2.02° 2.73 1.94 
816-16 40 961x107! 1008x1079 8,65 « 103 2.29 2.57 2.06 2.57 1.85 
18-9 (77) 9.28x103 . 1,232 1.58 1.27 
18-10 (79) 9.60x 103 1.207 1.55 1,25 
18-11 (79) 1.11x 104 1,115 1.55 1.25 
18-12 (79) 1.43 x 104 1.560 1,55 1.25 
27-1 (64) 2.91x 103 38,75t 1.82 1.39 
27-2 (64) 2.43 x 103 1,273 1.82 1.39 
27-3 64,5 2.44103 1,220 1.80 1.38 
27-4 (64) 7.79 x 10? 0.815 1,82 1.39 
35-17 (43) 4.41 10? 1,522 2.43 1.77 
35-18 (43) 7.13 « 10? : 1.873 2.43 1.77 
35-19 (42) 1.05 x 103 2.55 2,48 1.79 
35-20 38 1.10x 103 ' 1.280 2.70 1,90 
56-11 (23) 7.69 x 10? 3.37 3.90 2.67 
56-12 (21) 7.88 x 10? 3.97 4.20 . 2.87 
56-14 (20) 3.47x 10? 1.753 . 4,386. . 2.97 
*Values in parentheses are estimated, tValues not used in the analyses. 


Notes: (1) 1, is the ion-chamber measurements by the radiochemical project corrected to 100 hr and to absolute units 
by a factor of 0.788. ; 


Table 7—(Continued) 


Ay An Ay Ag Ag Kos Ko, 
fission fission fission fission fission Ig, r/hr 107% r/hr at Lhe r/hr atl hr 
sq it sq ft sq ft sq ft sq ft at 1 hr fission/sq it kt/sq mile 

2.781034 2.54 x 1014 2.63x10%  . 2,83x10" 2.70x10'4 26.0 0.963 301 
2.101014 2.631014 2.63x 10! 3.02x 1014 2.60x 1044 26.0 1.00 520 
1.46 x 104 1.98 1015 1,46 1045 1.63x 1045 206 1.26 657 
2.68 1015 2.98x 105 1.4419 2.03 x 1045 206 1.01 528 
2.30 x 101? 1.31 x10! 1.3110! 0.015 1.14 595 
1211x2018 1,08 x 1013 1.26 10% 1,18x 10% 1.16 1013 1.69 1.46 758 
8.76 x 1014 9.88x 10! 9,52x 10! 9.39 x 1014 98.4 1.05 545 
1,29 108 1,34 x 1015 9.88 10" 9.48x 10! 1.14= 10% 98.4 0.862 449 
171x101 9,81x 1010 : 9.81 10% 0.010 1.02 530 
3.16 x 104 2,92x 10 3.06 1044 3.05 x 104 30.5 1.00 520 
3.05 x 1014 3.15 x 104 3.06 x 1014 3.15 x 1oM 3.10x10!4 30.5 0.985 512 
3.54x 1010 3.53 «10 3.5410! 0.0035 0.990 515 
5.79 x 101° 5.0110! 5.01 x 10° 0.0052 1.04 540 
2.79101! 5,13 x 1010 5.13x 10! 0.0057 1.11 578 
2.57x 1012 3.05 x 10!? 3.42 x 10%? 3.01x 10! 0.30 0.998 519 
7.39 x 1048 8.9910! 8.11x 10! 1.10x 1045 8.55x 104 79.7 0.932 485 
584x104 5.50x 10! 4,18x 10! 5.17x 10! 54.4 1.05 547 
5.52x 104 5.4.x 10! 5.50 10" 5.46 x 10! 54.4 0.995 518 
5.58x 1014 5.45 x10! 5.501014 6.64 10!4 5.79 1014 54.4 0.940 489 
3.521010 1,96x 10"! 1.96 x 101 0,020 1.02 531 
5.95 x 1010 3.79 x 1010 3.79 1019 0.0040 1.06 549 
7.39 x 1012 5.58x 10" 4.61x 10! 5.09 x 10" 0.54 1.06 551 
7.14x 10" 5.55x 1044 5.35 x 1014 7.37x 10% 6.35 10! 53.9 0.848 441 
6.59 x 10!4 7.35 1014 4.60 10 6.18 1014 74,5 1.20 627 
§.32x 10" 7.35 x 10" 7.24% 10% 6.64x 10'4 74.5 1.12 584 
8.54 x 10! 8.14 x 1014 7.35 x10" 8.01x 10% 74.5 0.930 483 
2.65 x 10'4 2.22x10!4 2.47= 1044 1.62 10! 2.24x 10" 25.4 1.13 590 
1.64 x 101? 1,82x 19! 1.83x 10! 1,58x 10%? 1,72x 10! 0.18_ 1.05 545 
1,53 x 10!! 6.18% 10! : 6.18x 10° 0.0070 1.13 589 
5.74101? 4.5310! 9.96 101? 5.93 10!? 5.40 101 0.57 1.06 549 
4.18x 101° 4.31% 1010 4.24 1010 0.0044 1,04 540 
2.15x 1014 2.10 104 2.10 1044 2.09 104 2.11« 10! 21.5 1.02 530 
3.36x 1014 . 3.76 x 1044 3.83x 1014 3.65 x 10!4 39.6 1.08 565 
3.94x 1044 3.88x 10! 3.76 x 104 4.38 10! 3.86 x 10! 39.6 1.03 533 
9.48 1043 *8.99x 1013 1.25x 1014 - 1.76 1014 1.22 10! 14.7 1.20 626 
1.065 x 1014 1,29x 10 1,25x 10! 1.76x 10" 1.34x 10" 14.7 1,10 570 
2.44x 108 2,86 x 1013 6.50x 104 3.39x 1013 2.90x 10! 3.4 1,17 610 
9.74 x 1010 2.69 x 1010 2.69 10!0 0.0034 1.26 657 
1,07 1018 9.46 x10! 1.71x 1018 1.13 x 104% 1.05x 1019 1.4 1.33 693 
3.54x 1014 1.60x 104 1.63 10!4 2.19 10" 1,81x 10" 20.4 1.13 586 
1,55 x 1015 1.63 1014 2.16% 10%4 1.90x 1044 20.4 1.07 559 
1.60x 10% 1.72x 10" 1.40x 1044 1.81 10!4 1.63 x 10% 17.0 1.04 542 
1,60x 108 1,28x 10% 1,32 1048 1.41108 1.5 1.06 553 
1.61% 10! 1.0410!" 1.321091 0.017 1.29 670 
2.58x 10! 1.35 x 10! 1.96 x 1044 0.020 1.02 530 
2.23 x 1013 142x108 1.75 x 1013 1,82x 108 2.0 1.10 571 
6.80x 1010 1.19x of 1.19x 101) 0.015 1.26 655 
7.55 x 1013 6.38 1018 7.17 1018 6.61™ 10!3 6.93x 10!3 9.08 1.31 682 
7,.28x 1018 7.17x 1053 7.90x 10% 7:45 x 1013 9,08 1,22 634 
8.65x 101? 6.44x 1012 5.72 10!2 8.21~ 10! 7.25x10!2 1.58 2.18 1133 
3.21x 1018 5.95 x10" 6.20 1012 6.68 x 1012 6.28x 10!2 1,32 2.10 1094 
4,39 x 10!? 4.87x 1012 4.57% 10"? 4.61 10% 4.18 2.56 1330 
4,20x 10!2 3.74x 10% 4.20x 10? 4.68™ 10!2 4.20x 1012 1.17 2.79 1450 
7,54x 10! 7.56 x 1012 7.55 «10! 4,17 1.55 805 
7.95 x 10" 7.42™ 10! 742x102 1,18 1.59 830 
9.96 x 10!2 8.88 10! 8.88x 1012 1,32 1.49 115 
9.17« 103? , 1,14 1018 1.14x 10" 1.58 1,39 725 
7.52x 1010 2.09 1012 2.09 x 10!2 0.28 1.34 697 
1.91x 10? 1.75 x 10! 1.75 x 10! 0.47 2.68 1395 
2.00 x 10!2 1.77™ 10% 1.77x 10!2 0.27 1.52 793 
9.68x10" 5.60% 101! 5.60 10!! 0.15 2.68 1390 
2.90x 10! 2.49» 19013 2.49x 1011 0.12 4.82 2510 
3.81 101! 4.03 10" 4.03x 10! 0.18 4.47 2320 
4,12x10!! 5.87 101! 5.87x 10!" 0.27 4.60 2390 
8.60x 10! 5,82x 104 5.82x 191 0.26 4.47 2320 
2.28x 1041 2.88x 10!! 2.88x 103! 0.074 2.57 1335 
1,98x 10! 2.75104 2.75101! 0.16 5,72 3020 
1.98 10!! 1.17x 10! 1.5810! 0.087 5.50 2860 


Notes: (2) I, is the ion-chamber measurements by the field project (average for each station). (3) I’ is the calibrated 
gamma-scintillation-counter measurements by the field project. (4) i and i are values read from the curves in Figs. 2 
and 3, respectively. 
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The computed values of Ky are plotted as a function of ds) in 
Fig, 4, It may be notéd that the curves of i(100), i(100), and Ky given 
as a function of dsj) in Figs. 2, 3, and 4 are very similar to the curve 
for i(100) plotted as a function of d in Fig. 1. ; 
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Fig. 4—Variation of Ko with dsq for fallout-collecting stations. 


The high values of K, for the stations in the 35- and 56~numbered 
series (Table 7) suggest that the I, values are Oyen ee males by cor- 
rection of the monitoring data! ic 1 hr with the t~'” function. one 
adjustment of the two apparently different sets of reported A, values’ 
to -a midvalue rather than correction of one of the two sets of data to 
be consistent with the other may have resulted in Ky values that are 
either low or high by 25%, Further investigation of the reported data? 
together with theoretical ie wae | suggest what corrections are 
appropriate, 

The Dq terms of Eq. 12 can be evaluated using an appropriate 
value of k, > The selection of a value here is limited to the k,, value 
for the figsion products from the thermal-neutron fission of aaeTy. 
because the i(100) and -r,- values, as deduced from the ion-chamber 
measurements, are based on the response of the ion-chamber to the 
*85y products. The calculated value’ of k,, is 3950 r/hr at 1 hr per 
kiloton per square mile. The Dq of Eq. 12 is given by 


a Ky . . 
Dd = 35950 rs at) 
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The values of r,,(100), r/(100), ry, rf, Ko, and Dq for each 
station are summarized in Table 8. The value of K, in Table 8 for 
stations having more than one estimated value of K, in Table 7 is a 
median value. The terrain about the stations generally increased in 
roughness! proceeding from station 100 to station 707; however, this 
general increase in roughness is not reflected in a corresponding 
general decrease in Dq values for these stations, At the stations of 
the 800 series, the terrain was fairly smooth and flat at stations 813 
and 814 and somewhat rougher at stations 815 and 816. The terrain 
features at the remainder of the stations were varied. At the stations 
in series 100 through 800, the GITR’s. were calibrated to have a 
response factor, D, near unity; thus for these stations the computed 
value of Dq is practically equal to q. For the other stations the re- 
sponse factor for the instruments and monitor has not been reported, 
but it is assumed to be about 0.7. Therefore, where the computed 
value of Dq is in excess of 0.7, either the I, value is too large or the 
A, value is too small (i.e., the r;, value is too small). 


Table 8—SUMMARY OF CALCULATED rvp, rfp, Ke AND Dq VALUES* 
FOR THE SMALL BOY SHOT FALLOUT-COLLECTING STATIONS 


ike 
Station and ; r/hr at 1 hr 
sample No, T,,(100) Tf, (100) , es To kt/sq mile Dq 

100 0.338 0.369 0.194 0.304 511 0.667 
101 0.355 0.425 0.204 0.346 589 0.731 
103 (0.403) (0.443) (0,232) (0.361) 595 0.650 
200 0.456 0.469 © 0.262 0.382 758 0.724 
201 0.340 0.398 0.195 0.324 494 0.642 
202 (0.340) (0.400) (0.196). (0.326) 530 0.685 
203 0.341 0.400 0.196 0.326 515 0.666 
204 (0.343) (0.405) (0.197) ~ (0.330) 515° 0.663 
207 (0.363) (0.433) (0,209) (0.353) 540 0.678 
209 (0.383) (0.438) (0.220) = (0.357) 578 0.665 
300 (0.343) (0.352) (0.197) (0.287) 519 0.691 
303 0.335 0.361 0,192 0.294 485 0.640 
305 0.338 0.381 0.195 0.310 518 0.673 
309 (0.343) (0.405) (0.197) (0.330) 531 0.684 
31r (0.363) (0.433) (0.209) (0.353) ~ 549 ‘0.665 
400 (0.339) (0.352) (0.195) (0.287) 551 0.717 
401 0.344 0.354 0.198 0.289 441 0.564 
403 0.346 0,357 0.199 0.291 561 0.715 
405 0.350 0.362 0.202 0.295 590 0.741 
407 0.350 0.362 0.202 0.295 545 ~ * 0,684 
409 (0.383) (0.400) (0.220) (0.326) 589 ‘ 0.678 
501 0.362 0.378 0.208 0.308 549 0.669 
502 (0.355) (0.372). - (0.204) (0.303) 540 0.671 
503 0.350 0.362 0.202 0.295 530 0.665 


’ (Table 8 continues on page 63) 
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Table S— (Continued) 


Kg, 
Station and ' r/hr at 1 hr 
sample No. r,,(100) Vfp (100) lfp Tip kt/sq mile Dq 
505 0.359 0.372 0.206 0.303 549 0.674 
507 : 0.399 0.415 0.230 0.338 598 0.659 
509 (0.403) (0.419) (0.232) (0.342) 610 0.666 
513 (0.444) (0.452) (0.255) (0.369) 657 0.654 
601 (0.444) (0.452) (0,255) (0.369) 693 0.688 
603 0.428 0.443 0.246 0.361 572 0.590 
605 0.415 0.431 0.239 0.351 542 0.575 
607 (0.403) (0.419) (0.232) (0.342) §53 0.605 
700 (0.565) (0.571) (0.325) (0.466) 670 0.523 
701 (0.524) (0.524) (0.302) (0.427) 530 . 0.445 
703 (0.484) (0.476) (0.278) (0.388) 571 0.522 
704 © (0.444) (0.462) (0.255) (0.369) 655 0.644 
707 0.431 0.449 0.248 0.366 658 0.673 
813 0.573 0.557 0.330 0.454 1,133 0.870 
814 0.576 0.562 0.332 0.458 1,094 0.836 
815 1.100 0.924 0.633 0.753 1,330 0.533 
816 1.037 0.880 0.596 0.718 1,450 0.617 
18-9 0.638 0.605 0.367 0.482 805 0.556 
18—10 0.625 0.595 0.360 0.485 830 0.585 
18-11 0.625 0.595 0.360 0.484 775 0.546 
18-12 0.625 0.595 0.360 0.485 725 0.511 
27-1 0.735 0.662 0.422 0.540 697 0.418 
27-2 0.735 0.662 0.422 0.540 1,395 0.838 
27-3° 0.726 0.657 0.418 0.535 793 0.482 
27-4 0.735 0.662 0.422 . 0.540 1,390 0.835 
35-17 0.980 0.843 0.564 0.687 . 2,510 1.128 
35-18 0.980 0.843 0.564 0.687 2,320 1.043 
35-19 1.000 0.852 0.575 0.695 2,390 ..1.053 
35-20 ~ 1.090 0.905 0.626 0.737 2,320 0.939 
56-11 1.573 1.272 0.905 1.037 1,335 0.374 
56-12 1.695 1.367 0.975 1.113 3,020 0.785 
56—14 1.760 1.414 1.012 1.152 2,860 0.717 


*Values in parentheses are derived from estimated ds) values for the station. 


' The mean value of Dg for all the stations in the 100 to 700 series, 
for which ds is given, is 0.686 + 0.025. With this value of Dq, Eq. 12 
becomes 


Ky = 2710 ry, (22) 


Thus K3 for all the indicated stations is 2710 (i.e., Ky with lg = 1). 

The f, values for the Small Boy shot faliout pattern were computed 
by constructing fallout patterns for a given range in particle size and 
evaluating the intensity—area integral for each pattern. The fraction of 
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I, at each station contributed by each particle-size group was calcu- 
lated from activity measurements on sieved fractions of the fallout 
samples. The constructed fallout patterns for several of the particle- 
size ranges and for the total fallout pattern are shown in Figs. 5 to 8. 
The activity—size distribution derived from this analysis is shown in 
Fig. 9. The median particle diameter for the distribution is 210 yu. The 
derived distribution is not lognormal, 


The {; values from which the distribution ‘was determined are 
summarized in Table 9; also given in the table are the f,r,; and f,K), 
values leading to the estimates of Tf and Ky for the whole fallout 
pattern. The sum of the intensity—area integrals, Jp, of the fallout 
patterns for the different particle groups is 580 r/hr at 1 hr per 
square mile; the value of J,’for the pattern of the gross I, contours 
is 640 r/hr:at 1 hr per square mile, about 11% larger. The yield of the 
Small Boy shot has not been reported, except that the yield was small,’ 
With the use of the intensity—area integral for the I, contours and the 
yield, the value of J,/BW for the pattern thus could be estimated. With 
a K, value of 1070 r/hr at 1 hr per kiloton per square mile, the value 
of C could be estimated; and the fraction of device, F,, accounted for 
within the constructed fallout pattern of Fig. 8 could similarly be 
estimated. , 


The fallout patterns constructed for the various particle-size 
fractions do not include contributions from the regions of high I, 
values near ground zero nor from the large area of low-level fallout 
that was deposited in northern Utah, These contributions are not in- 
cluded in the analysis because fallout samples for determining the 
distribution of radioactivity among the fallout particles deposited in 
these regions were not available for analysis. For the total fallout 
pattern, including the contributions to I, from these areas, a Jp value 
of 1460 r/hr per square mile is obtained; this value is a factor of 
about 2.3 larger than that calculated for the fallout pattern constructed 
from the I, values at the stations for which particle-size data are 
available. 


The fraction of the device, Fp, accounted for within the larger 
fallout pattern, would therefore be 2.3 times that for the smaller 
pattern mentioned earlier. -The ratio of the contributions to the 
intensity—area integral for the high-intensity region to the low- 
intensity region, as determined from integrals, is about 3 to 1. If 


an average r, value of 0.29 is assumed for the fallout in the high-- 


intensity region and an average r, value of 1.5 for the low-intensity 
region, a value of K, can be estimated for the larger pattern from the 
data of Table 9. Allocation of appropriate fractions of the intensity — 
area integral to the three areas gives a weighted average value for Ky 
of. 1400 r/hr at 1 hr per kiloton per square mile; the corresponding 
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Fig. 6—Isointensity contours for 200- to 350-u particle diameters in 
milliroentgens per hour at 1 hr. 
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Fig. 8—Total fallout isointensity contours 
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Table 9—-SUMMARY OF THE f; VALUES FOR THE SMALL BOY SHOT 
FALLOUT AND THE CALCULATIONS FOR a AND Ky 
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value of C for the larger fallout pattern could then be calculated if the 
yield were known. , 

A commonly used value of K§ in Eq. 7 for estimating the fraction 
of a device in the fallout pattern, especially when experimental values 
of Dg are not directly available, is 2000 r/hr at 1 hr per kiloton per 
square mile. If this value is used with the Jp value, 1460 r/hr at 1 hr 
per square mile, for the Small Boy fallout pattern, the estimated value 
of F, could be compared to that obtained through use of the radio- 
chemical analysis of the fallout material and the methods outlined in 
this report. 


COMMENTS AND CONCLUSIONS 


The computed l tp values appear to be somewhat low, especially for 
the median particle diameters between about 200 and 800 u. Further 
examination of the radiochemical data with respect to the fractionation 
of individual nuclides and to the fission yield of ane as a measure of 
the number of fissions, could provide further information on the r, 
values, The major factors in determining the r,, values, if the data 
from ion-chamber measurements are accepted as being by far the 
most accurate, are the fission content of the samples and the decay 
factors from 100 to 1 hr after detonation for the ion-chamber data. 
The absolute value of iy for the unfractionated mixture of radionuclides 
from thermal-neutron fission of *U rather than for fission by fission 
neutrons would result in a relatively small error in the r,, estimates. 
It should be expected that the decay factors for the samples would 
approach those for i) as the r,, (100) values approach unity, indicating a 
relatively unfractionated mixture of radionuclides. 

The Dq values derived from the data are consistent with other 
previously derived values of q for the Nevada Test Site terrain.’ These 
values generally would not be influenced by possible future changes in 
the r,| values since any change in ty would result in an equivalent 
change in the values of the Kp. 

The example analysis of some of the Small Boy shot data applica- 
ble to evaluation of the intensity activity ratio shows that the ratio is 
not a constant for a given fallout pattern but varies over the pattern 
depending on the particle sizes in the deposited fallout. However, when 
gross activity—size distribution data are evaluated, an average value 
of the ratio can be derived for the fallout pattern. On the other hand, 
the average value of the intensity—activity ratio is not required for 
estimating the fraction of the device within a fallout pattern. A value 
for this ratio for the fallout deposited uniformly over an ideal plane, 
an average terrain attenuation factor, and an instrument response 
factor, however, are needed for estimating the fraction of device in the 
pattern. 


| 
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Further detailed analyses of the available data are needed to 


verify or adjust the: Values of several of the parameters, as described 
in this preliminary analysis of the reported data, for more accurate 
evaluations of the intensity—activity ratios and the fraction of device 
accounted for within the fallout pattern. 
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RADIOCHEMICAL-DATA CORRELATIONS 
ON DEBRIS FROM SILICATE BURSTS 


GLENN R. CROCKER, FRANCIS K. KAWAHARA, and EDWARD C. FREILING 
U. 8. Naval Radiological Defense Laboratory, San Francisco, California 


ABSTRACT 


Local-fallout samples collected in 1962 in the Johnie Boy, Small Boy, 
and Sedan shots were analyzed radiochemically for sr, *"sr, *ty, zr, 
99Mo, Ru, Ru, 1, Pe, 136Cg, log MBA 'iCe, 144c@ 23ND and 
239pu. These results, as reported by the project officers, have been 
correlated by plotting the log of the ratio of equivalent fissions of each 
nuclide, i, to equivalent fissions of Zr, Yr; .95, against the log of the 
similar ratio for “Sr and *Zr, roo 95. The data were fitted to straight 
lines by linear regression; and slopes, intercepts, coefficients of cor- 
relation, and confidence limits were determined. The slope of such 
a line for a nuclide i is an indication of the degree of fractionation of 
the nuclide relative to the fractionation of *Sr. For all cases observed 
in these Nevada shots, the same group of nuclides (Sr, *°Sr, *ty, !Ru, 
“106Ru, ty, 82pe Cs MBA Cs, MICe, and ?*Np) fractionated from 
*5Zr. The nuclides *°Mo, '“*Ce, and 2 Pu did not fractionate appreciably 
from Zr. For the Sedan shot the slope of the '*"Cs plot is near 1.0, 
but, for the remaining fractionating nuclides, the slopes lie within a 
narrow intermediate range. For the Johnie Boy and Small Boy shots, 
for which the results are quite similar, these slopes show much wider 
. variation. The results of the correlation have been compared with 
those from similar correlations for a coral-surface burst and some 
air bursts. Aside from the fact that °*Mo fractionated from Zr in the 
air bursts and that *8Np did not fractionate from °°Zr in the coral- 
surface burst, the results indicate that differences in the fractionation 
behavior of the nuclides are of degree rather than of kind. Correlation 
of radiochemical results with particle size indicates increasing degree 
of fractionation with increasing particle size. 
72 


o 


CORRELATIONS ON DEBRIS FROM SILICATE BURSTS 73 


INTRODUCTION 


In the summer of 1962, the U. S. Naval Radiological Defense Labo- 
ratory (NRDL) collected fallout from three shots at the Nevada Test 
Site (NTS)—Johnie Boy, Small Boy, and Sedan. Radiochemical studies 
were made of the debris from these shots and were reported.!-? Dur- 
ing the past several months, we have made log—log correlations of 
these radiochemical results and have attempted some interpretation 
of them. 

Before the results are discussed, it seems appropriate to review 
the shot conditions since important differences existed among the three 
events. These differences undoubtedly had effects on the fractionation 
behavior of the fallout. 


SHOT CONDITIONS 


Small Boy was a low-yield shot fired from atop a 10-ft-high 
wooden tower above alluvial soil in Area 5 at NTS. NRDL collected 
many fallout* samples of debris at 43 stations within 8.7 miles of 
ground zero and took other samples on a line 15 miles from ground 
zero. Further samples were supplied by the University of California 
at Los Angeles sampling network at distances of 18, 27, 36, 52, and 
70 miles. Four cloud samples were supplied by Los Alamos Scientific 
Laboratory. The discussion in this report is mainly restricted to 
samples from within 8.7 miles of ground zero and from the cloud 


‘samples. The number of samples used from each station varied, and 


many samples were subdivided by sieving into seven fractions that 


-were analyzed individually. A total of about 187 samples is dealt with 


here, all of which were analyzed for ®Sr, *°sr, Sty, and zr. In addi- 
tion, about one-third of them were analyzed for **Mo, '“Ru, ‘Ru, 


'8os, ics, Ba, Mice, “4ce, 2 Np, and 78%Pu. Some of this last 


group of samples were also analyzed for ‘I and '*Te. The numbers 
quoted do not include a fairly large number of radiochemical analyses 
on samples used for solubility studies. 

NRDL’s participation in the Johnie Boy and Sedan-.shots was on a 


rather limited basis. Johnie Boy was a low-yield burst 23 in. below the 


surface of basaltic material in Area 18 at NTS. Forty-four fallout 
samples from the area out to about 1Y, miles from ground zero and 
two cloud samples were studied radiochemically. All of these were 
analyzed for ®sr, *°Sr, °'y, and *Zr, and about one-third of them were 
analyzed for the long list of nuclides previously given for the Small 
Boy samples. 


*As distinguished from cloud and air-filter samples. 
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Sedan was the well-known 100-kt cratering shot for the Plowshare 


program. The device was buried 635 ft below the surface in alluvial 
material in Area 10. NRDL collected samples at stations within 1 to 
3.6 miles of ground zero. Analyses were performed on about 50 sam- 
ples for the short list of nuclides previously given, and about one-third 
of these were analyzed for the long list. 


TREATMENT OF THE DATA 


The radiochemical analytical work was contracted out to three 
commercial ‘laboratories since NRDL does not have facilities for 
handling a large volume of samples for routine analysis. The con- 
tractors were selected on the basis of qualification tests, which were 
also intended to serve for interlaboratory calibration purposes. It re- 
quired nearly a year for all the contractors to complete and report 
the analyses. The results were reported to NRDL as disintegrations 
per minute or equivalent ?%U thermal-neutron fissions at shot time. 
NRDL then converted the values to equivalent fissions of the device 
at shot time, using mass-chain yield values supplied by the weapons 
laboratories, At the same time the calibration factors derived from 
the qualification-test analyses were applied. 

As a reference nuclide Zr was chosen for fractionation studies. 
Some such choice seems unavoidable if fractionation is to be dis- 
cussed. Zirconium-95 belongs to a class of nuclides (other members 
are »°Mo, !4Ce, and !4"Nd) that are generally present in debris in equal 
quantities when expressed as equivalent fissions. Strontium-89, on the 


other hand, shows wide variation when compared with any of these. 


The ratio of equivalent fissions of ™Sr to equivalent fissions of Zr 
is, in fact, a fairly good indicator of the degree to which a sample is 
fractionated. The ratio will be 1 for an unfractionated sample, less 


than 1 for a sample depleted in "Sr, and greater than 1 for a sample | 
-enriched in ®Sr. 


The behavior of other fractionating nuclides can be correlated 
with that of “Sr by log—log plotting. The log of the ratio, r,; 95? of 
equivalent fissions of nuclide-i to equivalent fissions of ®Zr is plotted 
against the log of the similar ratio for ™Sr and ®Zr. The data so 
treated can be fitted more or less satisfactorily to a straight line. 
The slope of the line is then an indication of the extent of fractionation 
of nuclide i. 


RESULTS AND DISCUSSION 


For the Johnie Boy and Small Boy shots, the local-fallout collec- 
tions were sufficient to define the geographical extent of the close- -in 


wo 
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fallout fairly well. The Small Boy field was the cigar-shaped down- 


wind area typically associated with such shots. The Johnie Boy field 
was very, perhaps atypically, narrow with a very hot line down the 
center, which. was visible on the ground as a darkened streak. The 
NRDL collection array at Sedan was not sufficiently widespread to 
define the limits of the close-in fallout field. 

With each of these shots, one can associate a sort of average, or 
typical, value (or narrow range of values) for the ratio T5995 observed 
in gross samples. It is to be understood that this range of rg, 4, is much 
extended when cloud samples, sieve-fraction samples, and samples 
from the peripheral stations are considered. Nonetheless, the average 
value is useful for describing the overall degree of fractionation of the 
shot. For Johnie Boy a weighted average of this ratio for the hot-line 
stations is around 0.03, indicating very severe fractionation. For 
fmall Boy the values for most stations are in the range 0.1 to 0.2, 
indicating more moderate fractionation. All Sedan samples were 
sieved, but the T5995 values for gross samples were reconstructed by 
properly weighting the values for the sieve fractions and were found 
to range from about 0.5 to 3.8. Valuesin this range would be anomalous 
for local fallout from a true surface burst, but they seem to be charac- 
teristic of cratering (buried) shots, venting underground explosions, 
and venting underwater bursts. 

All of the Johnie Boy, Small Boy, and Sedan r, g, values have been 
log—log plotted against rgg 95. Figure 1 indicates how some of the data 
accommodate themselves to this treatment. These are the data on ‘Ba 
for the Johnie Boy shot, and they illustrate a particularly satisfactory 
fit. The slope, the y intercept at log x = 0 (x = 1), and the coefficierit of 
correlation* have been determined by linear regression to be 0.61, 
2.4, and 0.98, respectively. Not all.of the data, by any means, fit a log— 
log plot so neatly. Figure 2 shows a much less satisfactory and more 
typical example. These data are part of the results for 50Sr in the Small 
Boy shot. The coefficient of correlation here is only 0.684. 

The data were originally fitted by lumping together data points 
from all three contracting laboratories. In the process of investigating 
a poor fit of ‘Ce data for Small Boy, we noticed that points from 


Laboratory A could be fitted fairly well, whereas those from Labo- 


ratory B were scattered so badly as to suggest that Laboratory B was 
experiencing difficulty with this analysis. Figure 3 shows these data. 
Investigation of the ‘Ru and !8Ru data showed a similar situation, but 
it was reversed with respect to the two laboratories— Laboratory B 
points correlated notably better than did Laboratory A points. 


*The square of the coefficient of correlation reflects that part of the varia- 
tion in one set of measurements which can be explained by their dependence on 
the other. 
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Fig. 1—Fractionation plot for Ba in Johnie Boy. 
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Fig. 2—Fractionation plot for *Sr in Small Boy. 
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Fig. 3—-Fractionation plot for 4#Ce in Small Boy. 


All the data were reinvestigated by plotting and fitting the data 
points of each laboratory separately. Differences in the degree of 
correlation appeared which were sufficient to cast suspicion on part. 
of the data for !Te, cs and '"Cs. Calibration differences among 
the laboratories also appeared to be obscuring the: correlation in 
some instances. Figure 4 illustrates this effect. When the data from 
the three laboratories are treated separately, the three lines shown 
are obtained with a reasonable degree of confidence. Note that their 
slopes do not differ greatly. However, when the data are lumped to- 
gether, the correlation is much less clear. 

In spite of the difficulties involved in correlating the data, some 
facts emerge fairly clearly, as can be seen in Table 1. First, the 
same group of nuclides previously mentioned (®Sr, *°sr, *ty, ‘Ru, 
Ru, 341, Te, cs, cs, Ba, ICe, and.?%*Np) fractionated from 
*°Zr in all three shots, whereas *Mo, ‘4Ce, and ?%°Pu did not. Second, 
ics appears to fractionate about as severely as Sr. Beyond these 
basic points, the Sedan data display some differences from the data 
obtained in the other two shots. This is not surprising in view of the 
radical differences in the Sedan shot conditions. For Sedan the frac- 
tionating nuclides other than '*’Cs all showed a more or less inter- 
mediate degree of fractionation—the slopes of the log—log plots were 
between about 0.3 and 0.6—-whereas, for Johnie Boy and Small Boy, 
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Fig. 4—Correlation lines for *Y in Small Boy from three different 
laboratories. 
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*Certain reservations relative to these values are expressed in the text. 
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the range was considerably wider. The slope values obtained from the 
Johnie Boy and the Small Boy data are in quite good agreement with 
each other. The Sedan data are also given for purposes of comparison 
along with slopes from previous correlations on a high-yield coral- 
surface burst’ and on some air bursts. A few interesting points with 
regard to the table may be mentioned. The nonfractionating behavior 
of 38ND in the coral-surface burst is at variance with its behavior in 
the air bursts and the silicate bursts. On the other hand, °8Mo frac- 
tionates in the air bursts but not in the other shots. There is an in- 
dication of possible slight fractionation of '44Ce in Sedan and in the air 
bursts. In the silicate bursts '°Cg fractionates less than '"Cs. The 
fractionation of the latter is explainable on the basis of its rare-gas 
precursor, '*’Xe. Cesium-136, on the other hand, is shielded; i.e., it 
has no precursors. A possible explanation of its fractionation behavior 
lies in the basicity of the alkali-metal oxide or hydroxide relative to 
the silicate soil medium. 

A reservation should be noted for the slope values reported for the 
103-, 106-, and 132-mass chains for Johnie Boy and Small Boy. There 
is some indication in each of these cases that the data points belong to 
two different families and should not be fitted to the same straight line. 
The available data are unfortunately inadequate to resolve the question. 

Aside from the Mo and *Np anomalies just mentioned, insofar 
as the table permits one to judge, the differences in fractionation in 
the various kinds of bursts are of degree rather than kind. 

Analysis on many samples separated into sieve fractions displays 
the trend of fractionation with particle size in a rough way. Figure 5 
shows an example of these data from the Small Boy shot and also 
illustrates the variability of the ratios observed. The ratio of *Sr to 
Zr has been plotted against the sieve-opening diameter. There is, 
of course, no really accurate way of plotting these data. For example, 
any point on the right side of the 1580-y line corresponds to a sample 
that passed a No. 7 sieve (2800-p openings) and was caught on a No. 12 
sieve (1580-y. openings); presumably the sample may contain particles 
up to 2800. The scattering of the points is thought to be due to the 
combined effects of heterogeneity and very small sample size. At least 
two kinds of active particles, smooth spheres and irregular particles, 
were discernible in these samples. Since many of the sieve fractions 
consisted of only a few particles, their compositions may not have been 
truly representative of their particle-size classes. If an average 
weighted according to sample masses is taken, one obtains the points 
of Fig. 5 through which a line has- been drawn. The figure is quite 
representative of the behavior of fractionating nuclides in both Small 
Boy and Sedan. The sieve data on the Johnie Boy samples were too 
few to draw any conclusions relative to particle size. 
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Fig. 5—Effect of particle size on fractionation ratio. 
SUMMARY 


Correlations indicate that the fractionation behavior of *Sr, 30Sr, 
Sty, ics, 14Ba and '4ce relative to Zr in.silicate-surface bursts 
does not differ essentially from that observed in a coral-surface burst 
or in air bursts. Differences are noted for Mo and ?°Np. For Ru, 
M6Ru, 1, pe, Bos, “Ice, and Pu, some differences in degree of 
fractionation are indicated, and further study of these nuclides is 
recommended. 
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Richmond, California 


ABSTRACT 


Particles, isolated from high-altitude nuclear explosions, were sepa- 

~ rated and their physical parameters determined. The particles varied 
in color from colorless to black. They were generally spherical, and 
some had satellites that were a result of collisions of small particles 
with larger particles before complete solidification had occurred. Par- 
ticles with diameters greater than 2 » were isolated. 

The beta radioactivity of the particles was measured at 25 days 
and was found to be related to the diameter expressed in the following 
way: A = aD’, where b showed values of 3.06 + 0.08, 3.22 + 0.16, and 
3.01 + 0.30 for the three shots studied. Since the standard deviations 
were due to variations in activities of the particles and not in experi- 
mental errors, it was concluded that the activity was directly propor- 
tional to the particle volume. The standard deviations were found to be 
unrelated to particle size and color, but the value of a was found to be 
a function of the weapon yield. The value of b, however, was not related 
to weapon yield. 

Individual particles were allowed to decay to determine the decay 
rate of the fission-product mixtures. On the average, the particles 
followed the relation A = kt~!"!°. Values of the decay slopes varied from 
0.7 to 1.4. The distributions of the decay slopes of the individual par- 
ticles from the average was found to be Gaussian. The deviations were 
not a function of particle size or color; they are explained as being a 
result of the widely different radionuclide content of the particles. 
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PHYSICAL CHARACTERISTICS 


- The fireball of a nuclear explosion contains material from two 
sources, the weapon itself and the explosion environment. In studies of 
the physical phenomena occurring during cloud cooling and of the pre- 
cipitation of the vaporized ‘materials and radioactivities in the form of 
a particulate. aerosol prior to the Operation Dominic I test series, the 
processes have been hard to resolve because of incomplete vaporiza- 
tion of the environmental materials and the incorporation of new ma- 
terials into the fireball as it rose. During the Operation Dominic I test 
series, however, it was possible to study high-altitude nuclear explo- 
sions in which the material incorporated into the fireball was only that 
of the weapon itself plus its ballistic casing. This relatively small 
amount of material was completely vaporized by the heat of the nuclear 
explosion. 

As the fireball expands and cools, the vaporized materials con- 
dense in the form of small droplets that solidify, forming the small 
particles which we are examining in the investigation reported here. 
Figure 1 shows photomicrographs of particles isolated from high-yield 
air bursts from the Operation Dominic I series. These particles show 
the spherical-shape characteristic of solidified droplets. When viewed 
under the optical microscope, these particles show great differences 
in color, varying from colorless to black, This may be due to variation 
in matrix-element composition or to surface effects. 

The radioactive part of the nuclear debris consists of fission 
‘products, radionuclides produced by neutron reactions with the matrix 
materials, and the unreacted fissile material of the nuclear device. 
The majority of the radionuclides are incorporated into the particles, 
with the fission products contributing the dominant part of the radio- 

‘activities. The.composition of the fission-product mixtures varies with 
fission mode and time after fission. Table 1 shows the composition of 
an unfractionated *“y fission-product mixture 25 days after the time 
of fission. 

Figure 2 shows the activities of the isolated particles plotted as a 
function of particle diameter. The analysis of particle samples, 2. in 
diameter or larger, shows that the beta activity increases as the cube 
of the particle diameter. Three shots were analyzed in this study, and 
a least-squares fit of the logarithm of the- activity vs, the logarithm 
of the particle diameter of the form log A = b log D + a gave values of 
b of 3.06 + 0.08, 3.22 + 0.16, and 3.01 4.0.30. The standard deviations 
observed are due to variations in the activities of the various particles 
and not to experimental errors in the measurement of these activities. 

The cubic relation of beta activity to size is used to describe one 
of the important physical properties of particles, the activity per unit 
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Fig. i1—Pholomicrographs of particles. 
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Table 1— THE RADIONUCLIDE CONTENT OF AN UNFRACTIONATED 
2355 FISSION-PRODUCT MIXTURE AT 25 DAYS 


Activity, Activity, 
Radionuclide % Radionuclide % 
89Sy 5.4 131] 1.7 
Sp — Ny 0.05 133X¢ 4.75 
aly 6.4 BICs 0.05 
3577 6.5 140Rq 8.1 
35Nb 2.8 MOT a 9.2 
8\Mo—8Te 0.3 - 141Ce@ 9.1 
10Ry 3.8 144C@ 1.5 
106Ru 0.08 143 py 25.6 
129mTe 1.4 MINd 12.5 
132Te —132] 0.5 M7 Dm 0.2 


particle volume. Figure 3 shows the activity per unit particle volume 
as a function of particle diameter for two shots. As is to be expected, | 
the activities per unit particle volume tend to fall at a common average 
for each of the three shots studied, depending on the fission yield of 
the nuclear device. For each case, however, there are particles that 
vary considerably from the norm. This variation appears not to be 
related to particle color, and ona fractional basis it appears to be as 
large in large particles as it is in small particles. 

The beta-energy spectrum, of course, varies as a function of time 
as the composition of the fission-product mixture changes. Figure 4 
shows aluminum absorption measurements on a particle sample at dif- 
ferent times. Close to zero time the absorption curve approaches a 
straight line on semilog plots, but as time passes the absorption curve 
shows a shape of a low-energy component plus a high-energy component. 
When the sample has decayed 200 days, the major high-energy com- 
ponent is due to '4Pr, the daughter of '#4Ce. 

In the discussions so far, the similarities between particles have 
been noted. Striking dissimilarities become apparent, however, after de- 
cay. Figure 5 shows the beta decay of 7*°U and *°°Du fission-product 
mixtures aS measured with a methane-end-window proportional 
counter, Although it is readily apparent that the rate of decay of these 
fission-product mixtures varies with time, a relation of A=kt7’isa 
close approximation of the decay from 3 to 4000 days. 

‘The decay of isolated particles was followed and the data during 
the period 3 to 200 days were chosen to determine the decay rate by a 
least-squares curve fit. Figure 6 shows plots of the number of particles 
vs, the decay slope of those particles, As analytical aids the standard 
deviation of the decay slope, the standard error of the estimate, the 
variance in each variable, and the difference between each experi- 
mental value and the unfractionated fission-product mixtures, derived 
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Fig. 3—Radioactivity per unit particle volume for two shots as a func- 
tion of particle size. 
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Fig. 4—Aluminum absorption measurements as a function of time. 
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Fig. 5—Beta-activity decay of unfractionated 2354 and Pu fission- 
product mixtures. 
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Fig. 6—The number of particles observed with various decay slopes. 
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from the straight-line approximation, were calculated. For over 90% of 
the decay curves, the standard error in the exponent was less than 
+3%. The excellence of. the curve fit was further demonstrated by the 
magnitude of the correlation coefficient, which exceeded 0.99 for over 
90% of the decay curves. 

Values of the decay constant were 0.7 to 1.4, with an average value 
for all shots of 1.10. These variations were found among particles of 
approximately the same size from each shot. There was no correla~ 
tion between the decay slope and the particle size or the color. A sta- 
tistical analysis of the variance indicates that a negligible variation 
can be associated with experimental errors or statistical fluctuations, 
and it is concluded that the wide variations in the beta-decay rates 
result from widely varying fractionation patterns among the particles. 
These distributions appear to be Gaussian and independent of weapon 
sield. 

The physical characteristics of single particles can yield informa- 
tion concerning the radiochemical content. Although single-particle 
beta and gamma decay curves, the beta-energy spectrum (as well as the 
gamma-energy spectrum), and the activity per unit particle volume can 
give valuable information from the health-hazard point of view, the 
radionuclide concentrations as a function of particle size are of prime 
importance for an understanding of the fractionation process. In any 
laboratory program dealing with single-particle analysis, it is desira- 
ble to provide simple test procedures to select those particles which 
will yield particles exhibiting the interesting radionuclide composi- 
tions. A correlation of the radionuclide compositions with the physical 
parameters of the particles will provide a field test-particle selection 


procedure. The radiochemical composition of these particles was also 


investigated, and the program concerned will be reported in the next 
paper.* 


*P, Benson, C. E, Gleit, and L. Leventhal, this volume. 
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ABSTRACT 


The radiochemical fractionation characteristics from two Operation 
Dominic high-yield air bursts were studied. Particles measuring from 
2.5 to 20.0 u were isolated from filter-paper matrices by autoradio- 
graphic registration and confirmed by optical microscopy. After physi- 
cal measurements radiochemical and/or gamma-spectroscopy studies 
were performed on individual and aggregates of particles. Of particular 
interest was the distribution of volatile and refractory nuclides for each 
particle size. The nuclides radiochemically analyzed were **Mo, ®**Sr, 
gr BZp 1MmPe 12%pe BiCg MB, Mice Moe Ming %'y and Np, 
Gamma-spectra analysis was used to-determine the Zr and ‘!°Ba con- 
centrations in 334 particles. A computer technique for resolving these 
nuclides and °°Mo, 1035py 131], 132p@ Bigg Moe, Mog MING and Np 
was developed. 

The particles were found to be extremely fractionated. A logarith- 
mic fractionation correlation plot for °Ba was determined which com- 
pared within the limits of error with data plotted by Freiling for high- 
yield water bursts. The logarithms of Zr and “°Ba activities were 
plotted as functions of the particle diameters for the two shots. The 
slopes for the regression plot for =Zr were 3.41 + 0.18 and 3.09 + 0.37. 
These values imply a cubic relation. The slopes for the regression plot 
for “°Ba were 2.63 + 0.24 and 2,32 + 0.50, somewhat between a square 
and a cubic relation. These results point toward condensation of the re- 
fractory nuclides initially with subsequent precipitation of the volatiles 
on the particles at later times. 
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These experiments have proved that important radiochemical frac- 
tionation information can be obtained from single particles. Both the 
limitations and the advantages of this technique have been delineated. 


INTRODUCTION 


The 1961 Dominic test series presented the first opportunity to 
study radionuclide fractionation in single particles from high-yield air 
bursts. In these air bursts the condensation matrices were dependent 
only on the weapon composition. Masses of contaminating materials 
such as tower, water, and salt, which would present indeterminable fac- 
tors, were not present. The radionuclide distributions in selected par- 
ticles obtained from several tests in this series were analyzed radio- 
chemically and by gamma spectrometry. Of particular interest was the 
distribution of volatile and refractory nuclides for each particle size. 
The nuclides radiochemically analyzed were *Mo, Sr, "Sr, Zr, 
amt, Te, Blog “Ba Mice, Me, “Ing, 23%y, and “Np. It soon be- 
came apparent that the radiochemical method was too time consuming, 
and as an alternative a nondestructive gamma-spectrometric technique, 
which also yielded the required data on radionuclide behavior, was em- 
ployed. A computer technique was developed for resolving the spectra, 
and 334 of the particles were examined by this method. Analysis of the 
radionuclide concentrations in individual particles so far has been con- 
fined to the determination of correlations between the data and the par- 
ticle diameters. 


PROCEDURE 


The radioactive particles were located by autoradiography and iso- 
lated from the filter~paper matrix. After particle size and other physi- 
cal parameters were measured, the particles were subjected to sequen- 
tial radiochemical analysis (see Fig. 1). Particles ranged in size from 
4 to 20 jzand contained between 10° and 10!° fissions based on ??Mo 
analysis. Aggregates of particles representing groups 4to 5 p» and 4.5 to 
9u were similarly analyzed to determine radionuclides that were in too 


low abundance to be measured in single particles. A complete study was 


made on integral papers. The nuclides were assayed by standard count- 
ing techniques, and the results were converted to R* values for fission 
products or atoms for induced radionuclides. ° 


*An R value is defined as the ratio of the ratio (r) of calculated fissions in 
isotope (i) to fissions as calculated from %Mo to a similar ratio (r) determined 
from thermal-neutron fission of 235U; i.e., (2;.99)X/(rj,g9) th. 
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Fig, 1—Sequential separation of particle radionuclides. 


Gamma spectra of individual particles were obtained with a 3- by 
3-in, Nal detector and a 400-channel pulse-height analyzer. The parti- 
cles were analyzed as early as two days after the shot and thereafter 
once every two particle half-lives. ; 

In general, six spectra per particle were obtained over a two-year 
period. Particles studied ranged in size from 2.5 to 25 and contained 
10° to 101! zr fissions. The gamma spectra were normalized prior to 
gamma-spectra stripping by a computer program written in FORTRAN 


for the IBM-7094. This program has a capability for the determination. 


of *Zr, °Mo, "cs, “ce, '4ce, °Ba, and “Np, as well as for others, 
and is based on the method proposed by Carnahan.! The isotopes men- 
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tioned represent the predominant peaks in the spectra obtained. Figure 
2 indicates a typical gamma-spectrometer scan observed over a two- 
year period. The background radiation is the lowest line. The peaks of 
interest are singled out on the figure. 


RADIOCHEMICAL STUDIES 


Calculations of the counting-rate yield at zero time for a particle 
representing 10° fissions indicated that measurable activities would be 
obtained for the nuclides of interest if nofractionation is assumed (see 
Table 1). However, fractionation effects in single particles caused wide 


Table 1— THEORETICAL RADIOCHEMICAL SENSITIVITY 
FOR A PARTICLE OF 108 FISSIONS 


(Basis: Counting Rate at Zero Time for 100% Chemical Yield) 


Nuclide Counts/min Remarks 
89Sr 23 R=1 

Sr 0.1 R=1 

S2r 10 R=1 

iA 1000 110 abs, R= 1 
399Mo 562 R=1 

‘lag 0.5 R=1 

18Cq 2 R=1 

103Ru 2 Zero abs, R= 1 
106Ru 0.2 Zero abs, R= 1 
18Cs 0.05 R=1 

BICs 0.09 R=1 

40Ba 94 R=1 

141C@ 24 Zero abs, R=1 
44C@ * 7.4 Zero abs, R=1 
‘Sty 25 R= 

6ry 0.06 74 abs, R=1 
MINd 40 R= 

"Be 0.004 N/f* = 100 x 1075 
4Na 44 N/f = 100 x 1075 
®4Cu 200 N/f = 1000 x 107% 
8iCu 0.06 N/f=1x 1075 
18TW ; 0.23 N/f=1x 107 
203 pb 0.015 N/f=1x 10% 
31Th 900 : N/f = 100 x 1075 
23m 690 N/f = 0.3 

29Np 2880 N/f = 0.3 

2401 2880 N/f = 0.1 

238py 0.00005 N/f = 10 x 1075 
239 Pu / Jt 0.0009 N/f= 0.5 


*N'= atoms, f = fissions. 
+J = [atoms 89Pu/atoms 24°Pu/(atoms 23?Pu/atoms *49Pu) + 3.7]. 
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variations from predicted values. Particles ranging in size from 5 to 
20 . were analyzed from each of two shots. Two sets of selected aggre- 
gate particles between 4 to 9 were also analyzed.” 

.A logarithmic fractionation correlation plot of ry49,39 VS. T95,39* in- 
corporating filter-paper aggregate particles and individual particle data 
from both shots fitted linearly (see Fig. 3). Freiling® has plotted ryyo ¢9 
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_ Fig. 3—Logavithmic fractionation correlation for Ba. 


data for megaton-range coral-surface bursts and three types of water 
bursts. Within the limits of error, the slope and intercept of the re- 
gression line of Fig. 3 are comparable with those of Freiling’s plot. 

‘The single particles. are severely fractionated with respect to the 
total sample. It is interesting to note that for these particles the frac- 
tionation plot had to be extended another decade to include the particle 
r-value ratios. 


*1;,g9 is defined as the ratio of calculated isotope (i) fissions to calculated 
89Sr fissions, f,/fg;, where f equals fissions. 
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Table 2— PRINCIPAL GAMMA RAYS OF THE MAJOR FISSION 
PRODUCTS FOR A PARTICLE OF 108 FISSIONS 


(Basis: Gammas Per Minute at Zero Time) 


Radionuclide Energy, Mev Dis/min Gammas/min 
9577 — 95 Nb 0.717 47.0 47.0 
0.745 8.6 ; 8.5 
0.235 0.17 
9Mo—8mTe 0.780 1050.0 210.0 
0.740 146.0 
0.180 10.5 
0.041 105.0 
: 0.141 930.0 930.0 
3 Ry—103mMRh 0.498 38.0 36.0 
0.040 38.0 
1317 0.364 170.0 138.0 
132T¢ — 132] 0.231 670.0 670.0 
1.40 76.0 
0.96 130.0 
0.777 490.0 
0.673 670.0 
0.528 170.0 
1379, —BImpa 0.29 
0.661 0.270 
140 Bq — 1407 2 0.54 230.0 70.0 
0.30 23.0 
0.16 23.0 
0.03 230.0 
1.60 220.0 
0.815 67.0 
0.49 91.0 
'41Ce 0.145 . 89,0 60.0 
143Ce 0.660 2100.0 325.0 
0.356 260.0 
0.289 1285.0 
0.126 260.0 
14Ce— 144pr 0.134 10.0 : 10.0 
0.081 10.0 
0.034 10.0 
1.480 0.20 
0.695 0.41 
MINg 0.532 100.0 25.0 
0,318 1.50 
0.092 60.0 
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Similar graphs were prepared for 406 and *°"y. The plot of the 
particle data for M44o9 also had to be extended a decade higher. The 
uranium activity per particle was too low to get accurate data, anda 
linear fractionation plot could not be developed from these data. 

Analyses were made for Sr, '%cs, '*Te, ce, and *Np, but they 
were not detected in the samples. Satisfactory values were obtained for 
Pompe, “4ICe, and “Nd. 


GAMMA-RAY SPECTRA 


The individual radionuclides that are susceptible to gamma-spectra 
analysis in particles were determined. Limits of detection were ob- 
tained for the radionuclides of interest by assuming 10° thermal-neutron 
fissions. This was done by correction of the estimated gamma rays per 
minute for the efficiency of the 3- by 3-in. Nal crystal at a particular 
energy. Further’ selection was made by determining if interfering 
gamma energies were present that produced unresolvable spectra in the 
region of interest. Of the heavy elements, 235Np was the only nuclide 
high enough in yield for analysis. However, because of its short half- 
life, no attempt was made to resolve this nuclide. The induced activities 
were too low for spectral analysis. Among the fission products listed in ” 
Table 2, **Mo also was too low for spectral analysis. The radionuclides 
‘Ru, 1, and ‘Cs were not detected, but Zr and “Ba were re- 
solved. Also '*Te, '!ce, '4Ce, and '4’Nd can be resolved. 

The Ba and “Zr fissions for two high-yield shots have been cal- 
culated! from gamma spectra obtained from 21 particles from shot A 
and 45 particles from shot B. Peak resolution was confirmed by decay 
of the respective peaks. Three of the particles were sacrificed after 
gamma stripping and were radiochemically analyzed for '°Ba and "Zr. 
Agreement between, physical and radiochemical measurements was + 5%. 

The logarithm of '°Ba fissions was plotted as a function of the 
logarithm of the particle diameter (Figs. 4 and 5), and a linear re- 
gression calculation was performed. The Slopes of the regression line 
obtained for the shots A and Bdata were 2.63 + 0.24 and 2.32 + 0.50, re- 
spectively. The hypotheses that the populate slopes are square or cubic 
(equal 2.0 or 3.0) were tested. Because of the scatter of the points about 
the regression line, both hypotheses could notbe rejected except for the 
hypothesis that the shot A slope was two, which was rejected at the 95% 
confidence interval. The intermediate value may be explained as due to 
the decaying of volatile precursors before particle condensation. The 
slopes obtained for the regression plots for 7, (Figs. 6 and 7) for 
shots A and B were 3.41 + 0.18 and 3.09 + 0.37. For this nuclide a cubic 
relation appears to hold. An explanation for the large value for 57 r for 
shot A is not apparent. Spectra stripping of additional particles may re- 
duce this value. 
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Fig. 4—Barium-140 activity as a function of particle 
diameter for shot A. 
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Fig. 5—Barium-140 activity as a function of particle 
diameter for shot B. é 
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Fig. 6—-Zirconium-95 activity as a function of particle 
diameter for shot A. 
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Fig. 7—Zirconinum-95 activity as a function of particle 
diameter for shot B,. 
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Although the surface relation between 404 activity and particle 2 
size cannot be demonstrated on the basis of these data, the variance in 
the activity per unit area is quite obvious. One would expect this in nu- 
clear-explosion clouds in which the particles form rather quickly and 
the volatile constituents precipitate on the particle surfaces much later. 
This hypothesis appears to be confirmed if one considers the *Zr as 
representative of a refractory nuclide. The cubic relation then:-would 
indicate that the Zr is condensed initially andis distributed as a func- 
tion of the volume of the particle. In effect, regularity is observed 
through approximately a hundredfold change in volume. 

Another interesting observation that stems from a comparison of 
the shots A and B plots is that the intercept of the ‘4°Ba regression line 
of the nuclides of a particular particle size is a function of the yield to 
mass ratio. 


CONCLUSIONS 


The feasibility of isolating and analyzing individual particles by 
radiochemical and gamma-spectrometric techniques has been demon- 
strated. : 

Particles under 10 u and as low as 10° fissions have yielded accu- 
rate values for the important volatile and refractory fission products. 
Heavy-element analysis is borderline for *°"U and **Np and is de- 
pendent on analysis at early times. Unless an induced activity is pres- 
ent in abnormally large amounts, the analyses of these radionuclides 
in single particles are not feasible. 

Gamma spectrometry is a satisfactory technique for the resolution 
of nuclides in large numbers of individual particles and may be used to 
obtain °Ba and *Zr values to + 5%. 

Logarithm fractionation plots of particles from high-yield air tests 
have been found to correlate well with rj4o,39 and Yrs g9 values from a’ » 
coral-surface burst and three types of water bursts. 

Barium-140 fissions have been found to correlate with between a 
square and a cubic relation to particle diameter. Zirconium-95 fissions 
have demonstrated a cubic relation to particle diameter. 
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PREDICTION OF FALLOUT 
FROM SUBSURFACE NUCLEAR DETONATIONS 


JOSEPH B. KNOX 
University of California, Lawrence Radiation Laboratory, 
Livermore, California 


ABSTRACT | 


A numerical simulation model has been developed for the prediction of 
fallout from subsurface nuclear detonations that produce craters 
through spall and the action of the cavity gas. The physical processes 
modeled are atmospheric transport, lateral eddy diffusion, and gravi- 
tational sedimentation of radioactive particulates. This cratering fall- 
out model is normalized to the observed external gamma-dose-rate 
fields of the Sedan (100 kt) and the Danny Boy (0.43 kt) cratering shots 
conducted at the Nevada Test Site, Calculations of the fallout patterns 
for additional shots, used for testing the prediction capability of the 
cratering fallout model, indicate that the model gives estimates of the 
external gamma dose rate at H+ 1 hr with a maximum error of a fac- 
tor of 2 to 3 in the gamma dose rate vs. distance along the hot line of 
the pattern, 


INTRODUCTION 


During the past few years, a small but continuing effort has been 
expended in developing a model for predicting fallout from subsurface 
nuclear detonations., In this report the experience in predicting fallout 
from surface bursts that is transferable to the construction of a fallout 
model for subsurface detonation is summarized, the development of the 
fallout model for subsurface nuclear detonations (the cratering fallout 
model) is described, the predictive capability of the model is illustrated 
by means of independent test cases, and some of the problems asso- 

331 


332 KNOX 


ciated with the prediction of fallout from. row-charge subsurface nu- 
clear detonations are discussed. This paper is limited to research or 
development in which the investigator has been personally involved 
within the Plowshare Program of the Lawrence Radiation Laboratory. 


CRATERING FALLOUT MODEL 


Basis of Model Construction 


For prediction of the area affected by radioactive fallout from a 
subsurface detonation of a nuclear explosive and of the gross external 
gamma dose rate in the surface-fallout pattern, knowledge of the fol- 
lowing factors are needed: 

1. The heights of the base, the top, and the radius of each radio- 
active cloud (i.e., the main cloud and the base surge) formed by the 
detonation at the time the clouds cease to rise in the atmosphere, This 
time is defined as the time of cloud stabilization. The cloud heights are 
prescribed in terms of height above surface zero. 

2. The total yield of the explosive, W_; the fission yield of the ex- 
plosive, W,; the depth of burial of the explosive, z; the fraction of the 
fission-product gamma emitters expected to appear in the fallout pat- 
tern beyond the estimated radius of direct ejecta, F.; and the equivalent 
fission-yield gamma needed to simulate the gamma dose from induced 
activities, é 

3. The activity—particle size distribution in both the main cloud 
and the base surge and the fraction of F, in each. 

4, The terminal fall velocity of the fallout particles (in still air) 
as a function of particle size and height in the atmosphere. 

5. The time and space prediction of the horizontal wind at the level 
of each cloud top along with the specification of the wind-shear tangen- 
tial and normal to the wind for the layer through which particles fall. 
This knowledge is required because in a fallout calculation involving 
two .clouds (main cloud and base surge), the calculation of the fallout 
pattern for each cloud is done separately. The total surface-fallout 
pattern of the shot is found by summing the patterns from the base 
surge and the main cloud. 

6. The effect of horizontal eddy diffusion on the growth of the 
horizontal radius of the disks of radioactive particles as the disks 
fall earthward. The initial debris cloud is subdivided into disks of 
debris in the model as a function of initial height in the cloud and of 
particle size. 


The preceding information (items 1 through 6) provides the sim- 
plest, but still adequate, basis for constructing a cratering fallout 
model. In addition to this information, it is assumed that the fission- 
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product radioactivity is unfractionated and that 1 kt of unfractionated 


fission products spread uniformly over 1 sq mile corresponds to an 
(H+ 1)-hr dose rate of 3380 r/hr at a height of 3 ft above an infinite 
plane.! This normalization constant of 3380 r/hr is corrected for ter- 
rain shielding by a factor of 0.8. a 

Certain atmospheric processes or effects have been neglected in 
developing the cratering fallout model, just as these same effects were 
neglected in developing the land surface-burst fallout model:” 

1. The effect of synoptic-scale vertical motions on the vertical 
displacement of the falling disk of particles. 

2. The effect of the mean divergence of the horizontal wind on the 
radius of the disk of falling particles during descent. 

3. The effect of vertical eddy diffusion. 

4. The time from detonation to cloud stabilization. : 

5. The development of subsynoptic-scale wind systems. 

6. Changes in time and space of the shear tangential and normal 
to the horizontal wind at cloud-top level. 

7, The effect of water products of condensation on the size and fall 
rate of radioactive particles, 


Fraction of Gamma Activity in the Close-in Pattern 


For calculation of the close-in fallout pattern from a subsurface 
detonation, the fraction of the gamma activity produced by the fission 
yield and appearing in the close-in fallout pattern, F,, must be known. 
Thus far, experimental data from the Sedan, the Teapot ESS, the Jangle 
U, the Neptune, the Jangle S, and other surface-burst shots have pro- 
vided a basis for estimating F.. The measured fallout patterns from 
these shots have been integrated’ from the radius of direct throwout to 
the limit of the measured pattern. The results are shown in Table 1. 

These F, calculations were performed by assuming a normalization 
constant of 3380 r/hr per kiloton of fission yield per square mile at 
H+1 hr and a terrain shielding factor of 0.8. The graphical presenta- 
tion of F, as a function of z/W _ is shown in Fig. 1. In construction of 
the experimental F. curve, the following asymptotes were used: (1) an 


Table 1— EXPERIMENTAL DATA 


Shot W, kt z, ft FE. Medium 
Sedan 100 635 ~0.10 Alluvium 
Teapot ESS 1.2 67 0.46 Alluvium 
Jangle U 1.2 17 0.64 Alluvium 
Neptune 0.115 100 0.005 Tuff. 
Jangle S 1.2 0 0.50 Alluvium 
Danny Boy 0.43 109 0.04 Basalt 
Blanca 19 835 0.0005 Tuff 


334 KNOX 


10° 


° 
JANGLE U 


JANGLE S : : 
~ASYMPTOTE FOR 
FL AS z>0 


107! 


DANNY BOY q 


PRELIMINARY EMPIRICAL - 
Fe 107? F. CURVE FOR 1 KT 


NEPTUNE 9 


ASYMPTOTE FOR NO 
DYNAMIC VENTING 


BLANCA 


-lov4 
107! 10° 10! 107 10° 


2/w"/3, ery Ti“ 


Fig. 1—Fraction of gamma activity appearing in the close-in fallout 
pattern. : 


asyinptote of F, = 0.50 for z = 0, suggested by the Jangle S data in Table 
1 and supported by previous work,’ and (2) an asymptote of no dynamic 
venting for z/w* = 330 ft/kt”, supported by experimental evidence’ and 
by studies of containment physics.°® ; 

It should be noted that the value of F, (estimated by pattern inte- 
gration) for the Sedan shot departs considerably from the curve fitted 
to. the whole sample of F, data. This departure could arise because of: 
an F, yield dependency that is inadequately known at this time. In the 
absence of knowledge of such a dependency, F, estimates made from 
the curve in Fig, 1 for high-yield cratering events (of the order of 100 
kt and above) should be considered as uncertain by a factor of 2. It 
should be further stated that F, estimates derived from Fig. 1 are 
made with the tacit assumption that the nuclear explosive is fully 
tamped and that no additives for fission-product gamma-radiation 
suppression have been placed around the explosive. 

It should also be mentioned that the maximum value of F, (0.75 at 
z/W * © 30) occurs at very nearly the same value of z/W” as the maxi- 
mum base-surge radius (crosswind) in a neutral atmosphere, reported 
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in Ref. 7. The author believes that this coincidence of maximums at 
z/W* = 30 is physically consistent. 
Cloud Geometries 


The geometrical definition of the top, the base, and the radius of 
both the main cloud and the base surge at the time of cloud stabiliza- 
tion is shown in Fig. 2. In addition to the definitions given in Fig, 2, the 
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Hinr MAIN- CLOUD HEIGHT 
Rp, BASE -SURGE RADIUS 

Rr MAIN- CLOUD RADIUS 


BASE SURGE | 


SURFACE ZERO i Rp SFI 


Fig. 2—Definition of cloud dimensions and symbols. 


| 


height of the base of the main cloud is defined as being equal to the 
height of the base~surge‘top, Hy, in the model. 

Evidence suggests that the geometry of these two clouds at the 
time of stabilization is a function of the total explosive yield, the ma- 
terial in which the detonation occurs, the depth of burial of the explo- 
sive, and the meteorological conditions existing during the development 
of the clouds.’ At present, the cloud-geometry parameters. (R,, H,, Rin, 
and H,,) must be evaluated experimentally as functions of total yield 
and depth of burial, Reasonable samples of experimental data exist for 
alluvium and basalt materials. Examples from one of the most useful 
summaries of cloud-geometry data for alluvium® are shown in Figs. 3b 
to 3g. (Figure 3a is a computational aid to the acquisition of input to 
Figs. 3b to 3g. In these figures, zdenotes depth of burial and D,z denotes 
the depth of apparent crater.) This summary utilizes all the known 
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Fig. 3a—Computational aid to the acquisition of input data to the cloud- 
dimension graphs of Figs. 3b to 3g. 
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Fig. 3b—Base-surge radius as a function of the total yield and the pa- 
rameter 2/Dg. 
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Fig. 3¢ —Base-surge radius as a function of the total yield and the pa- 
rameter z/Da. 
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Fig. 3d—Base-surge height as a function of the total yield and the pa- 
rameter z/Da. 
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Fig. 3e—Base-surge height as a function of the total yield and the pa- 
rameter 2/Da. 
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Fig. 3f—Main-cloud radius as a function of the total yield and the pa- 
rameter 2/Da. 
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Fig. 3g—Main-cloud height as a function of the total yield and the pa- 
rameter z/Da. 


‘cloud-geometry data from high-explosive and nuclear-explosive detona-~ 
tions conducted by the U. S. Atomic Energy Commission in alluvium 
and basalt. Implicit in the summary is the assumption that a high ex- 
plosive and a nuclear explosive detonated in the same material at iden- 
tical depths of burial and under similar meteorological conditions pro- 
duce the same cloud geometries, 


Activity-Particle Size Distributions 


In a typical subsurface nuclear detonation in alluvium, two clouds 
are formed. The main cloud is composed mostly of vented cavity gas 
and particulates (originating from either condensation or injection and 
entrainment of soil). The base surge is composed of ejecta and sus- 
pended fine particulates. For a nuclear cratering shot in alluvium, it 
is assumed that 80% of the F, gamma activity is in the main cloud and 
20% is in the base surge. The 0.8 F, main-cloud activity is assumed to 
be subdivided between two lognormal activity—particle size distribu- 
tions, The first activity particle size distribution contains the activity 
0.8 w,(1)F, and is characterized by the mean ln f,,(1) and the standard 
deviation o,,(1). The second activity—particle size distribution con- 
tains the activity 0.2 w,,(2)F, and is characterized by the mean In F,,(2) 
and the standard deviation o,,(2). The activity of the first distribution 
is assumed to be homogeneously mixed through the whole main cloud, 
whereas that of the second distribution is assumed to be homoge- 


340 KNOX 

neously mixed only in the lower fifth of the cloud. A similar prescrip- 
tion of activity—particle size is used for the base surge. Figure 4 
summarizes the parameters governing the activity—particle size dis- 
tributions in the two clouds at the time of cloud stabilization. These 
parameters governing the activity—particle size distributions in the 
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Fig. 4—Schematic drawing of an idealized cloud from a subsurface 
detonation showing spatial relations of the activity—particle size dis- 
tribution assumed in the model. 


cratering fallout model have been determined by mathematical experi- 
mentation with the model in the re-creation of the observed fallout 
patterns for the Sedan and the Danny Boy shots. Results of these cali- 
bration calculations will be discussed in a later section. 


Terminal Fali Velocity of Fallout Particles 


The vertical fall velocities of the fallout particles are modeled as 
the terminal fall velocities of smooth spheres of density 2.5 g/cm? in 
an International Civil Aeronautical Organization standard atmosphere 
as computed by McDonald® for both the Stokes’-law region and the 
aerodynamic region (wherein the Reynolds number exceeds 1). If some 
fallout particles are a cluster of small spheres attached to a large 
central particle, these complex particles are assumed to fall with the 
speed of the equivalent smooth spherical particle of the same mass. 
Experimental evidence has been obtained by Rapp and Sartor’® to sup- 
port this assumption. 
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Specification of the Horizontal Wind Field for the Model 


The horizontal wind field that transports the debris-disk centroids 
during their fall to the earth’s surface may be specified in two ways in 
the cratering fallout model: 

"1. Idealized Wind Hodograph. If H denotes height above surface 
zero, then the horizontal wind, v,, at height H for a simple wind hodo- 
graph (see Fig. 5) is 


A(p) 
A(p,,) 


V;,(H) = Vn (H,,) 


~ S(H,, — H)n 


where v,(H,,) = steady-state wind at cloud-top level H,, (or it can be 
specified as a function of time and space in either wind 
component form or by means of a stream function) 

A(p) = wind-shear component tangential to the horizontal wind 
at cloud-top ‘level and is evaluated from shot-time 
winds in a diagnostic fallout calculation or from pre- 
shot wind information for a predictive fallout calcula- 
tion (it is held constant in time for the period of fallout 
deposition) 

A(p.,) = value of A(p) at cloud-top level (normally it is set equal 
to 1) 

S = wind-shear component normal to the horizontal wind at 
cloud-top level (it is evaluated from shot-time winds or 
preshot wind information, depending on the purpose of 
the fallout calculation, and is held constant in time) 

n = unit vector normal to y,(H,,) in a right-handed system 

p = atmospheric pressure corresponding to H 


2° Arbitrary Hodograph. The horizontal wind can be specified in 
wind-component form for as fine a vertical interval as desired or for 
which wind information exists. 


Debris-disk Radius as a Function of Time 


For estimation of the radius of a debris disk expanding by hori- 
zontal eddy diffusion during its fall to earth, it is proposed that the 
disk radius as a function of time R,(t) be represented by 


Re(t) = (Re,o + 2Dt’)* 
where t’ is the distance traveled by the disk centroid divided by the 
mean horizontal wind speed in the layer through which the disk has 
Settled, R.» is the debris-disk radius at time of stabilization, and D 
is the horizontal eddy-diffusion coefficient.? The diffusion coefficient 
D is estimated as the Richardson’s diffusion coefficient 0.2 x }*, where 
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Fig. 5—Simple wind hodograph assumed in the 
model for wind specification. ({dealized wind 
hodograph.) 


Ris the standard deviation of the position of the particles (in the disk) 
from the disk centroid. Since } is poorly known in nuclear-debris 
clouds, it is usually set equal to Req or 2R.,) if accelerated relative 
diffusion is to be approximated in the fallout model. 


Physical Processes Simulated in the Model 


The principal physical or meteorological processes simulated by 
the cratering fallout model are (1) the transport of the debris disks by 
the mean wind in the layer through which the disks are falling, (2) the 
relative advection of the debris disks by the horizontal wind field con- 
taining both a speed and a directional shear, and (3) the lateral eddy 
diffusion, which expands the disks falling earthward. 

The first two processes are simulated by calculating the transport 
of the disk centroids by the ambient horizontal wind field during the 
disk’s earthward fall until the disk centroid is on the ground surface 
by using either the idealized wind hodograph or the arbitrary hodo- 
graph for the horizontal wind specification. The result of this calcula- 
tion is the position vector from surface zero to the predicted ground 
position of each disk centroid and the time of arrival at the ground 
surface of each disk. The lateral eddy-diffusion process is modeled by 
the expression for R,(t) which estimates the debris-disk radius at the- 
time of arrival of the disk centroid on the ground surface. 


where I is equal to 0, 1, 2,.. 
ticle modeled in the cloud; and r, is the radius of the smallest particle 
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The debris disks tracked earthward are defined as follows: In each 
cloud (main cloud and base surge), 11 disks of particle size r, are de- 
fined at each isobaric surface Py such that 


= 


, 10; r, is the radius of the largest par- 


modeled in the cloud. The term p, is equal to p, — [P/10(p, — pm)], 
where P is equal to 0, 1, 2,..., 10; p, is the pressure at the base of 
the cloud at the time of cloud stabilization; and p,, is the pressure at 
the top of the cloud at the time of cloud stabilization. 

The (H + 1)-hr external gamma dose rate (for a height 3 ft above 
an infinite plane) is calculated by using the method of Batten, Iglehart, 
and Rapp!! modified to account for the effect of normal shear and 
lateral eddy-diffusive disk growth. 

The following quantities are calculated in the model and are output 
in the indicated modes:* 
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(H + 1)-hr dose rate as a function of distance 

: along the! ‘hot line feach cloud and total cae 
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DIAGNOSTIC CALCULATION FOR SEDAN = 


For calibration of the cratering fallout model on the Sedan shot, 
the observed shot-time winds, the observed cloud geometry of the main 
cloud and the base surge, the estimated F, = 10%, and the appropriate 
fission yield were input to the model. A first guess of the 12 activity— 
particle size parameters, to be discussed later, was also input to the 


* About 0.25 min of IBM 7094 computer time is required for the calculations 
for a two-cloud nuclear-cratering fallout problem, and 2.4 min of Livermore 
Advanced Research Computer time is required for the cathode-ray-tube dis- 
plays indicated above. 
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model, The main cloud was assumed to contain 80% of F.; the base 
surge, 20%. The parameters governing the activity—particle size dis- 
tributions have been determined from the observed Sedan gross gamma 
fallout pattern by mathematical experimentation. The values of these 
parameters are as follows: 


w,, (1) = w,(1) = 0.9 
W,, (2) = w,(2) = 0.1 
In F,;(1) = In ¥,(1) = 2.9 


ln r,,, (2) = In r,(2) = 5.0 
 Om(1) =0,(1) = 0.69 
O (2) = 0,(2) = 0.59 
Figure 6 shows the calculated and the observed gamma dose rate* at 


H+ 1 hr from fission products vs. distance along the hot line for the 
Sedan shot. 


DIAGNOSTIC CALCULATION FOR DANNY BOY 


The Danny Boy shot was a 0.42-kt nuclear cratering detonation 
emplaced at a depth of 109 ft in dry basalt. The observed gamma fall- 
out pattern for Danny Boy has been published.’ By mathematical ex- 
perimentation with the cratering fallout model, the activity—particle 
size distribution parameters can be adjusted to duplicate the observed 
Danny Boy fallout pattern. In this shot no visible main cloud was ob- 
served.’ Thus, in the diagnostic calculation with the cratering fallout, 
model, it is assumed that 100% of the F, gamma activity is in the base 
surge and that there is a preliminary value of 0.05 for F.. The activity — 
particle size distribution parameters determined in this calculation are 
given below: 


w,(1) = 0.9 
~ ow, (2) = 0.1 | 
In F,(1) = 3.0 
In F,(2) = 5.7 
o, (1) = 0.69 
o,(2) = 0.59 


_ *The observed (H + 1)-hr gamma dose rate for fission products derived 
from the observed total gamma dose rate at H + 24 hr, assuming that 52% of the 
(H + 24)-hr gamma dose rate was from tungsten.!? 
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Fig. 6—The calculated and observed gamma dose rate at H+ lhrasa 
function of distance along the hot line of the Sedan pattern (diagnostic 
calculation). 


Figure 7 shows the calculated (H + 1)-hr dose rate vs. distance curve 
from the model, and the observed (H + 1)-hr dose rate vs. distance. As 
was previously mentioned, the value of F, used in the cratering-fallout- 
model calculation was 0.05. If a value of 0.04 for F, (as reported in 
Ref. 4) had been used, the agreement between the calculated (H + 1)-hr 
dose rate and the observed (H + 1)-hr dose ratevs. distance would have 
been better than shown in Fig. 7. 

Recently effort has been expended to obtain machine capability of 
plotting the calculated (H + 1)-hr dose-rate patterns for an arbitrarily 
Selected interval of dose rate. Figure 8 shows the machine-plotted dose- 
rate pattern for the Danny Boy diagnostic calculation. The computer- 
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pattern breadth at 25 miles downwind is about 7 miles; the observed- 
pattern breadth at 25 miles downwind is 5.5 miles. It should be noted 
that the closure of the isodose-rate line of 1074 r/hr at H+ 1 hr in the 
plotted pattern is artificial and is the result of the logic used for the 
computer plotting rather than of the logic used for the cratering fallout 
model, All isodose-rate lines will apparently be closed at the downwind 
edge of the pattern if computed dose-rate informationis insufficient for 
their appropriate extension downwind. The observed (H + 1)-hr gamma- 
dose-rate pattern for Danny Boy is shown in Fig. 9 for comparison with 
the computer-plotted pattern. 


INDEPENDENT TESTS OF THE MODEL 


_ An independent test calculation of the cratering fallout model (cal- 
ibrated on Sedan) was performed by using data from the Teapot ESS 
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Fig. 8—The machine-plotied (H+ 1)-hr gamina-dose-rate pattern for 
the Danny Boy shot (diagnostic calculation), roentgens per hour. 


shot. The observed winds at shot time, the observed cloud geometry, 
the published fission yield of 1.2 kt, and a value of 0,85 (an early pre- 
liminary value) for F, were input ito the model. Figure 10 shows the cal- 
culated and the observed (H+ 1)-hr gamma dose rate as a function of 
distance along the hot line of the pattern. A comparison of these two 


_ dose rate vs. distance curves shows that the largest error between 


calculation and observation is of the order of a factor of 2.5. An exam- 
ination of the radiosonde observation near shot time at the Nevada Test 
Site indicates that the vertical temperature distribution of the layer 
through which the particles were falling was slightly superadiabatic. 
Under such conditions, it is possible that the vertical eddy diffusion on 
the day of the Teapot ESS shot was larger than average. If enhanced 
vertical eddy diffusion were operative, the dose rate would have been 
slightly less near ground zero and would have been enhanced downwind. 

Other model confirmation tests have been performed. These include 
calculations with the surface-burst version of the model on the Apple II 
and the Zucchini shots? and on several other atmospheric shots for 
which the data are still classified. The results of these additional tests 
were satisfactory. 
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Fig. 9—Preliminary gamma-isodose-rate contours observed at H+ 1 
hr.4 Readings in roentgens per hour minus H+ 1hr. Data taken from 
Nuclear Defense Laboratory ground surveys (close in); Edgerton, Ger- 
meshausen & Grier,Inc., aerial surveys (intermediate range); ae U.S, 
Geological Survey aerial survey (long range). 
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Fig. 10—The calculated and the observed gamma dose rates at H+ 1 
hy as a function of distance along the hot line of the pattern (predictive 
calculation) from the Teapot ESS shot. 


PREDICTION OF FALLOUT FROM ROW-CHARGE SHOTS 


In principle, the fallout pattern from a row of subsurface nuclear 
detonations may be estimated with the cratering fallout model, pro- 
vided the model input parameters are adequately known for row-charge 
events. In the current state of knowledge, the cloud-geometry input 
parameters and the value of F, appear to be the most difficult to specify. 
Study of time-lapse photography of the main-cloud and the base-surge 
evolution from past high-explosive row-~charge tests give significant 
experimental information on the cloud-geometry parameters. For ex- 
ample, empirical methods of predicting crosswind radius and height of 
the base surges originating from high-explosive row-charge shots have 
been studied.’ In this study it was shown that for five equal-weight 
charges, equally spaced and emplaced at the same depth of burial, the 
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resulting base surge has approximately the same radius as the base 
surge from a large single shot of the same total yield and emplaced at 
the same scaled depth. The base-surge height for a five-charge row 
shot is reasonably well approximated by scaling the height of the base 
surge for a single-charge event by the 0.2 power of the total yield of 
the row~charge shot. 


A first approximation to the geometries of the individual main 
cloud may be obtained by treating each main cloud independently and. 
estimating the top and the radius of each cloud from the work of Day.® 
Examination of the high-explosive row-charge documentary photographs 
of the Rowboat, the Dugout, and the Pre-Buggy shots indicates that such 
an approximation is reasonable. This approximation can, of course, be 
in error if main-cloud interactions occur. There is a need to evaluate 
the uncertainty in fallout prediction for multiple-charge shots in cases 
where cloud interactions lead to the injection of radioactivity at levels 
higher in the atmosphere than predicted. 


Concerning the F, for row-charge shots, preliminary experimental 
results for high-explosive single- and row-charge shots have been re- 
ported previously.'® In this study, results of the measured vented frac- 
tion of “°La tracer from single- and row-charge high-explosive shots 
were given. The experimental evidence suggests that the vented frac- 
tion from row-charge shots may be about ‘twice that from single- 
charge shots. One of the most pressing needs for fallout prediction 
from row-charge shots is the establishment, through either experi- 
ment or theory, of the dependency of the F, on yield, depth of burial, or 
charge spacing. 

The development of computer aids for the row-charge fallout- 
prediction problem, however, can proceed independently of the solution 
of the two previously’ cited problems: (1) the specification of initial 
cloud geometry and (2) the specification of F,. Therefore the capability 
of plotting fallout patterns for multicloud and/or multidetonation events 
has been developed. Figures 11 and 12 show, respectively, the (H + 1)- 
hr dose-rate patterns computed for 10 Danny Boy detonations on an 
east to west line with charge centers separated by 33.5 m and for 10 
Danny Boy detonations on a north to south line with charge centers 
separated by 33.5 m. For the purpose of these calculations, each deto- 
nation is assumed to vent 5% (e.g., F, = 0.05), each cloud is assumed 
to be the same as that for Danny Boy, and the input wind for each cloud 
fallout problem is assumed to be the same as the shot-time wind for 
Danny Boy. These two row~charge fallout patterns, although detectably 
different as determined from the printout, appear to be very similar. 
Calculational results suggest that the fallout pattern from a small- 
scale row~-charge shot is not sensitive to the orientation of the wind to 
the alignment of the charges. 
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Fig. 11—The (H+ 1)-hr gamma-dose-rate pattern computed for 10 
Danny Boy detonations on an east to west line with charge centers sep- 
arated by 33.5m. Values are in roentgens per hour. F 
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Fig: 12—The (H+ I)-hr gamma-dose-rate pattern computed for 10 
Danny Boy detonations on a north to southline with charge centers sep- 
arated by 33.6 m. Values are in roentgens per hour. 
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CONCLUSIONS 


In conclusion, it may be stated that the cratering fallout model de- 
veloped gives reasonable results for the prediction of the area of the 
pattern, the geometry of the dose-rate contours, andthe (H + 1)-hr dose 
rate vs. distance along the hot line. The independent testing of the 
model should be extended to include more than the Teapot ESS case 
described and the several others cited. 

Meteorological improvements in the cratering fallout model could 
well include the following items: 


1. The prediction of synoptic-scale Changes in the normal- and 
tangential-shear components. 

2, A better understanding of the growth of the debris disks through 
horizontal eddy diffusion during their earthward fall. 

3. The prediction of diurnal change of wind at low levels and close 
to ground zero. 

4, The inclusion of the effects of topography on the evolution of the 
horizontal wind field. 

5. The effect of terrain-induced éiedilations on fallout deposition. 

6. A better solution to the cloud-rise.problem for cratering deto- 
nations. 


Before embarking on the development of these meteorological im- 
provements in fallout predictions, however, one must consider potential 
improvements that may come from other areas. Areas of promise are 
cratering physics and filtration theory for the vented-fraction problem, 
special emplacement for control of vented fraction, and improvement 
of nuclear explosives, 
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DISTRIBUTION OF RADIOACTIVITY WITH HEIGHT 
IN NUCLEAR CLOUDS 


GILBERT J, FERBER 
U. S. Weather Bureau, Washington, D. C. 


ABSTRACT 


During Operation Dominic I at Christmas Island in 1962, aircraft 
sampling of nuclear clouds was done soon after cloud stabilization 
to investigate the amount of radioactive debris that stabilizes in the 
troposphere and its distribution with height. The detonations studied 
were all air bursts over water. Some data for surface bursts obtained 
during Operation Redwing in 1956 are used for comparison, Results 
indicate that for air bursts less than 1% of the total radioactivity is 
present in the stem of the nuclear cloud. It is estimated that about 
‘one-third of the total debris from the Christmas Island clouds initially 
stabilized in the troposphere. 

Project Stemwinder has shown that in-cloud dose-rate monitoring 
by aircraft is a relatively simple and economical way to obtain informa- 
tion on the distribution of radioactive debris in a nuclear cloud, Used 
in conjunction with limited radiochemical analysis of samples, this 
type of monitoring could produce a reliable inventory of the debris in 
a nuclear cloud. 


INTRODUCTION 


The objective of Project Stemwinder was to probe and sample nu- 
clear clouds as soon as possible after cloud stabilization to investi- 
gate the amount of radioactive debris that stabilizes in the troposphere 
and its distribution with height. Sampling was done with RB-57 air- 
craft of the 1211th Test Squadron under the scientific direction of the 
Atmospheric Radioactivity Research Branch, U. S. Weather Bureau. 
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The detonations investigated were all air bursts over water during 
Operation Dominic I at Christmas Island in 1962. Some data for sur- 
face detonations obtained by aircraft sampling during Operation Red- 
wing in 1956 are used to compare withthe Project Stemwinder data, 

The project was conceived as an attempt to utilize available sam- 
pling aircraft to narrow the area of uncertainty involved in two re- 
lated problems. The first problem was concerned with the operational 
need for prediction of the possible local hazards due to rainout of 
radioactive debris from a portion of a nuclear cloud that might pass 
over Christmas Island shortly after an air burst. Since the tops of 
rain clouds in the Christmas Island area were generally below 20,000 
ft. and often below 10,000 ft, the amount and distribution of activity in 
the stem of the mushroom cloud was of primary concern. The second 
problem was concerned with the partitioning of nuclear debris between 
the stratosphere and the troposphere as a function of nuclear yield, 
tropopause height, burst height, and, possibly, other factors. Such 
partitioning has been an important consideration in estimating the 
long-range fallout from nuclear tests since nuclear debris has a mean 
residence time of several weeks in the troposphere (intermediate fall- 
out) aS opposed to many months or years in the stratosphere (worldwide 
fallout), depending on the latitude and altitude of injection. The frac- 
tion of the debris which remains in the troposphere may be particularly 
important in considering the possible hazards from relatively short- 
lived nuclides such as ‘7 since the stratospheric portion usually de- 
eays to insignificant amounts before it can return to the surface of the 
earth. 

It is emphasized that the preceding remarks apply only to the very 
small particles that contribute to the intermediate and the worldwide . 
fallout. In the case of surface detonations, much of the radioactivity is 
associated with relatively large particles that comprise the local fall- 
out. These large particles are not affected by the tropopause and 
appear in the local fallout regardless of whether they are initially in- 
jected into the troposphere or into the stratosphere. : 


CLOUD HEIGHTS OF AIR BURSTS IN A TROPICAL ATMOSPHERE 


Operation Dominic I shot data (including those for yield, burst 
height, and cloud-top, base, and tropopause heights) are given in the 
Project Stemwinder final report.! Since there was no scientific pro- 
gram to document cloud heights, a “best guess” was made for each 
cloud by evaluating estimates made by observers on the ground and in 
the sampling aircraft and by using the dose rates reported at the 
various sampling altitudes to verify, where possible, the visual ob- 
servations. Variations in the burst heights did not appear to have any 
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consistent effect on the ‘cloud heights. Evidently, the effect of the burst 
height was masked by the influence of meteorological factors and/or 
the errors in the cloud-height estimates. 

Selected data? were added from other Pacific test series to aid in 
drawing the mean curve and the curves for the estimated range of 
cloud heights (see Fig. 1). Almost all the detonations in previous U.S. 
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Fig. 1—Cloud -top heights and estimated range for air bursts (180Y°-4 
< burst height < 0.15 Hy) in a tropical atmosphere. 


tests in the Pacific had been surface bursts, and the documentation of 
nuclear cloud dimensions had been generally poor. The curves shown in 
Fig. 1 are intended to be valid only for air bursts in a tropical atmo-- 
sphere and for burst heights less than about 15% of the expected cloud- 
top heights. For this purpose, an air burst may be defined as a detona- 
tion at an altitude equal to or greater than 180y°*, where Y is the total 
- yield in-kilotons. It is emphasized that there are no reliable cloud-top 
data for yields greater than about 5 Mt, and the extrapolation of the 
curves beyond this point represents little more than an educated guess. 
Indeed, over the entire range of yields shown in Fig. 1, the dashed 
curves indicate only the expected range of cloud heights for the stated 
conditions and should not be interpreted as representing absolute 
limits, ; 
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STEM-CLOUD PENETRATIONS 


An RB-57 aircraft was available for stem-penetration missions 
immediately following seven of the Dominic Idetonations. The navigator 
was provided with a dose rate meter with a range from 0.01 to 2000 
mr/hr and was instructed to record the dose rate as the aircraft 
penetrated the stem cloud at specified altitudes. The dose rates mea- 
sured in the cockpit were then used to estimate the amount of activity 
in the cloud, : 

The relation between cloud concentration and dose rate in a uni- 
form infinite cloud’ is given by 


Cc Dp 8405 
~ po (3.7 x 104)(1.6 x 10 “JE 


(1) 


where C = cloud concentration, pc/em*® 
D = dose-rate, r/sec 
Py = Standard densi of air at sea level = 1.293 x 107 g/cm? 
o = density of air at sampling altitude, g/cm? 
E = average gamma energy, Mev 
84 = energy absorbed per roentgen, ergs per gram of air 
1.6 x 10°® = ergs/Mev 
3.7 x 104 = dis/sec per microcurie 


en 


If the unit of dose rate is converted to roentgens per hour and the con- 
centration to megacuries per cubic mile, 


Soe 
C= 217 (2) 


For stem penetrations within an hour after the burst, E was as- 
sumed to be 1 Mev. For the sampling missions between 2 and 5 hr after 
the burst, a value of 0.86 Mev was used.’ Figure 2 gives the value of 
p/P) as a function of altitude for a typical tropical atmosphere.” With 
the use of appropriate values for E and p/pg in Eq, 2, the dose rates 
recorded during stem penetration were converted to cloud concentra- 
tions. An estimate of the stem diameter was then used to estimate the 
total volume of cloud in a 1000-ft-thick layer, The total amount of ac- 
tivity in the layer and the fraction of the bomb represented by that 
activity was then determined by multiplying the concentration by the 
volume. The results are shown in Fig. 3 as a plot of the fraction of the 
bomb present in a 1000-ft-thick layer of the stem cloud vs. height 
(indicated as percent of the total stem height). The three highest 
Dominic I data points are derived from the extended sampling missions 
described later in this paper. The Redwing data used for the figure 
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Fig. 2—Ratio of air density (p) 
to standard sea-levelairy density 
(ppt= 1.293 x 107% g/cm’) as a 
function of altitude. 
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Fig. 3——Variation with height of the percent of the total activity resid- 
ing in a 1000-ft-thick horizontal layer of the cloud. 
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also are discussed later. The curve is intended: to represent a con- 
servative estimate (for safety considerations) of the stem activity as a 
function of height for air bursts. 

The rather large scatter in the data may be attributed to several 
factors. It appears that stem visibility may be the most important of 
these. Most of the higher activity readings occurred during penetra- 
tions when the stem cloud was visible to the pilot. The relatively low 
readings were obtained when the cloud was not visible or when it was 
not known whether the cloud was visible. It is quite possible that the 
aircraft did not actually penetrate the stem on these occasions. The 
dose rates measured inside the aircraft may have been due to “shine” 
from the stem cloud or to activity from diffuse material outside of the 
stem core. In those cases where the stem was not visible and where 
several passes were made at the same altitude, only the highest read- 
ing has been used. 

Unfortunately, for the lower 80% of the stem, virtually all the data 
for the larger detonations are questionable becauSe of the stem-visibil- 
ity problem. Therefore it is impossible to say whether or not the low 
stem activity found for these shots indicates a real decrease in the 
fraction of activity in the lower part of the stem with increasing nu- 
clear yield. . 

The following factors also contribute to the uncertainty in the 
results: 


1. Stem-volume estimates. For the determination of the total ac- 
tivity present in a 1000-ft-thick layer, the stem diameter at the pene- 
tration altitude was estimated. The values used were based on visual 
estimates made by observers on the ground and in the sampling air- 
craft or, where necessary, on estimates for other detonations in the 
same yield range. The estimated diameter could be in error by as 
much as a factor of 2 in some cases. 


2. Stem-height estimates. The stem was considered to extend 
from sea-level to the base of the cloud regardless of the burst height. 
The cloud bases used were based on visual observations from the 
ground and from the sampling aircraft and verified, where possible, 
by radiation readings taken in the sampling aircraft. The uncertainty 
in the height of the cloud-base (stem height) is about. 10%. 


3. Representativeness of dose-rate readings. The measured dose 
rates are assumed to represent those in a uniform, infinite cloud. The 
assumption appears to be reasonably valid for those penetrations 
where the stem was visible. The aircraft required 20 sec or more to 
traverse the cloud at a speed of about 7 miles/min while the mean free 
path of gamma radiation in air is on the order of a few hundred feet. 
The navigator reported that the dose rate usually rose sharply on en- 
tering the cloud, remained fairly steady (within a factor of 2) during 
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penetration, and then dropped sharply. (It would be advantageous to 
use automatic time—intensity recorders in future operations.) 

The effect of aircraft shielding on the dose rate in the cockpit is 
also uncertain. Tests made at the ground, using a point source outside 
the aircraft, indicated that there was no appreciable shielding effect 
on gamma radiation by the aircraft. Equation 2 assumes that the re- 
ceptor is completely surrounded by a uniform radiation field. Actually, 
of course, the receptor was surrounded by a “blank space” equivalent 
to the volume of the aircraft. No attempt has been made to correct for 
this. However, the effect should be small, probably less than a factor 
of 2, since the mean free path of the gamma radiation is large com- 
pared to the dimensions of the aircraft. Experimental determination of 
the correction factor should be planned in connection with any future 
operation of this type. 


AIRCRAFT SAMPLING IN THE VICINITY OF THE CLOUD BASE 


Aircraft equipped with Los Alamos Scientific Laboratory air- 
filter tanks were used for sampling after five Dominic I detonations. 
Approximately 1-hr sampling missions were flown at altitudes from 
35,000 to 48,000 ft at 2 to 5 hr after detonation. The two sampling 
tanks were-opened simultaneously when contact with the cloud was 
made and remained open for the entire sampling period. As the sam- 
pling patterns were flown, readings in the cockpit were made at 1-min 
intervals with a hand-held AN/PDR-27J3 Radiacmeter capable of mea- 
suring activity in a range from 0,01 to 500 mr/hr. Sampling missions 
were successful after four of the five detonations. The radiochemical 
analyses of the samples are reported elsewhere. ° 

The dose-rate readings obtained during three extended sampling 
missions were sufficient for estimations of the distribution and the 
amount of activity in the clouds at the sampling altitudes. Similar 
analyses were done for all three missions, The results are included in 
Fig. 3. 5: , 

One of the clouds was Sampled at an altitude of 45,000 ft at approxi- 
mately 3 to 4 hr after detonation. The base of the cloud was reported 
to be at about 45,000 ft. Both the shot-time wind data and the position 
of the cloud indicate east-southeast cloud travel at about 15 knots. For 
correction for the movement of the cloud during the sampling period, 
the reported aircraft positions were adjusted to the sampling midtime 
of 3, hr after detonation. The corrected radiation field and the actual 
unadjusted sampling track are shown in Fig, 4. If a decay exponent of 
-1.2 is assumed, integration of the pattern yields 520 r/hr-cu mi 
at 1 hr in a 1000-ft-thick layer. From Eq. 2 this is equivalent to 270 
Mc, or 4.9 x 107 of the total fission products produced by the detona- 
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Fig. 4—Sampling track and radiation pattern in a nuclear cloud at. 
45,000 ft. Values given are in millivoentgens per hour at H+ 373;hr. 


tion. The cloud covered an area of 5200 square miles at the 45,000- 
ft altitude. co 

The sampling track appears to have covered the cloud very well. 
However, the pilot reported shine from higher portions of the cloud 
during the last few minutes of sampling when the aircraft was outside 
the visible cloud. The readings due to shine were about 100 mr/hr. It 


- is possible that a significant fraction of the in-cloud readings may also 


have been due to shine from the upper portions of the cloud; therefore 
the estimate of the amount of activity at 45,000 ft may be high. 


COMPARISON WITH RADIOCHEMICAL RESULTS 


The radiochemical analyses of the samples obtained on the three 
extended sampling missions provide a check on the method of calculat- 
ing the amount of debris present from the dose-rate readings in the 


cloud. The total number of fissions collected in each sample has been 


determined,’ based on the number of atoms of *Mo present corrected 
for the fission yield of ‘Mo for thermal fission of **U, If the total 
volume of air passing through the sampling tank and the average dose 
rate along the sampling path are known, the fissions per sample can be 
estimated by the same method that was used to estimate the fraction of 


ay 
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the bomb in a 1000-ft-thick layer. The volume sampled is determined 
from the altitude, the air temperature, the aircraft speed, the sampling 
time, and the sampling-tank and filter-paper characteristics.® The aver- 
age dose rate is determined from the readings taken in the cockpit at 
1-min intervals during ‘the sampling period. From Eq. 2 the gamma 
megacuries in the sample can be calculated. If it is assumed that 1 kt 
of fission (1.4 x 10% fissions) is equivalent to 550 gamma megacuries 
at 1 hr, the following conversion factor can be used: 1 Mc (H+ 1)= 
2.64 x 107° fissions. Table 1 gives the pertinent data and the fissions 


Table 1— COMPARISON OF CALCULATED AND ANALYZED 
FISSIONS PER SAMPLE 


Fissions 
Sample Average Calculated per sample 
* Altitude, volume, dose rate, fissions (radiochemical 
Mission ft cu ft. mr/hr atH+1 per sample analyses) 
A 43,000 1.06 x 108 190 1.9 x 1014 3.3 x 10!4 
B 48,000 1.11 x 108 270 2.3 x 1014 4.9 x 1014 
Cc 45,000 1.10 x 108 560 5.6 x 1014 5.8 x lo! 


per sample as calculated from dose-rate readings and as determined 
from radiochemical analyses of the samples. 

The agreement between the calculated values and the results of 
the sample analyses is remarkably good, considering the uncertainties 
due to the possibility of shine from other portions of the cloud, air- 
craft shielding, and aircraft contamination. The calculated values for 
the samples from missions A and B are low by about a factor of 2, 
possibly because of the effect of the blank space previously mentioned. 
The calculated value for the sample from mission C is in almost per- 
fect agreement with the result of the sample analysis. Mission C is the 
one depicted in Fig. 4 and for which there was reason to suspect a shine 
contribution to the dose rates which may have compensated for the 
blank-space effect. 

' Additional experimental data are needed to evaluate all the factors 
involved, but the results indicate that the method employed on these 
missions is a practical and promising way to obtain the distribution of 
activity in a nuclear cloud. 


REDWING IN-CLOUD DOSE-RATE DATA 


The doses and dose rates at various altitudes in several nuclear 
clouds (all but one from surface bursts) were investigated by aircraft 
penetrations’? during Operation Redwing in 1956. Some of these pene- 
trations were complete traverses through the cloud. Since the altitude, 


N 
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the mean speed of the aircraft, the time in the cloud, and the average 
dose rate are reported, these data can be utilized in the same manner 
as the Dominic I stem penetrations to compute the device fraction con- 
tained in the cloud at the penetration altitudes. The computed device - 
fractions are plotted in Fig. 3 for comparison with the Dominic I data. 

Several interesting features may be noted. It appears that the 
activity in the upper half of the stem is greater for surface bursts 
than for air bursts and that the difference increases with altitude. The 
largest gradient of activity with altitude appears at about 70 to 80% of 
the stem height; this implies that for surface bursts the “radiological 
base” lies below the visual cloud base. However, this inference may 
not be warranted since the high activities encountered below the base 
may be due to the descent of fallout particles. 

The values computed for the lower portion of the mushroom indi- 
cate about 1 to 2% of the total fission products per 1000 ft. Since the 
mushroom portion of the clouds investigated averaged about 30,000 ft 
in vertical extent, the average activity in the mushroom must have been 
about 3% per 1000 ft. Thus we have some basis for believing that this 
admittedly crude method can give at least the right order of magnitude 
for the activity at a given altitude, even when the average dose rate re- 
corded during a single pass through the cloud is used. 

We note that the one Redwing data point for an air burst gives 
about five times the activity indicated by the curve estimated from the 
Dominic I data. This greater activity might be attributed to the fact that 
the detonation took place at a lower scaled height than any of the 
Dominic I air bursts. The burst height was somewhat below the mini- 
mum altitude for a true air burst according to our definition (burst 
height > 180Y°"*); therefore the activity distribution might be expected 
to be intermediate between those for air bursts and surface bursts. 
Although the close-in fallout measured after this detonation was very 
light, it was considerably more than that found after any Dominic I 
shot, where shipboard dose rates never exceeded 0.1 mr/hr. However, 
since only one surface vessel was available for fallout measurements 
during the Dominic I tests, the very limited number of measurements 
obtained does not permit the drawing of firm conclusions. 


RESULTS AND CONCLUSIONS 


Activity in the Stem Cloud for Air Bursts 


Although the Dominic I stem-penetration data leave a good deal to 
be desired for defining the distribution of activity in the stem, the curve 
in. Fig. 3 represents a best estimate based on our interpretation of 
these data. A major uncertainty lies in the assumption that the dis- 
tribution in the stem does not vary with yield. As indicated previously, 
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Fig. 5—Cumulative activity as a function of height in the nuclear cloud 
for airy bursts. 


this curve may represent an overestimation of the activity in the lower 
part of the stem for the larger yields (above about 200 kt). 


Cumulative Activity with Height in the Nuclear Cloud for Air Bursts 


With the use of the stem-activity curve in Fig. 3, an estimate of 
the cumulative activity with height in the nuclear cloud was derived. 
The solid portion of the curve in Fig. 5 was obtained from the stem- 
activity curve by using an average stem height of 40,000 ft and by as- 
suming the stem-top height (or visual cloud base) to be 63% of the 
cloud~top height (the average for the Dominic I series), Since the en- 
tire stem appears to contain less than 1% of the total activity, it is 
obvious that the activity must increase rapidly with height at or above 
the base of the cloud, The dashed portion of the curve represents a 
subjective estimate (based, in part, on the Redwing data for surface 
detonations) of the distribution of activity in the mushroom portion of 
the cloud. The activity in the mushroom is assumed to be distributed 
as follows: 


Layer, % of cloud-top height Fraction of total activity, % 
65 to 70 0.6. 
70° to 75 14 
75 to 80 25 
80 to 85 25 
85 to 90 15 
90 to 95 15 
95 to 100 3 
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For air bursts it appears reasonably certain that less than 1% of 
the total activity is present in the stem and that less than 0.1% sta- 
bilizes between the surface of the earth and one-half of the cloud-top 
altitude. The fraction of activity per unit altitude increases with height 
throughout the stem, and the region of maximum vertical gradient, 
which might be termed “the radiological base of the cloud,” probably 
occurs somewhat above the visual cloud base. The peak activity per 
unit altitude is assumed to occur between 75 and 85% of the distance 
from the surface to the cloud top, It is also assumed that for air bursts 
the distribution of activity relative to the cloud-top height does not 
vary with nuclear yield, burst height, or atmospheric conditions, 

Actually, the interaction of these factors must exert some in-— 
fluence on the activity distribution, The estimated top and base heights 
of the Dominic I clouds indicate that the ratio of base height to top 
height has a tendency to decrease with increasing yield. However, it is 
uncertain to what extent these indications are valid since the variation 
among detonations of about the same yield is almost as great as that 
for the range of yields from 40 kt to several megatons. The mean 
ratio is 63%, with individual clouds varying from 53 to 73%. Some of 
the variation may be attributed to errors in the estimates of the base 
and top heights, but part of it is undoubtedly real. There is a similar 
uncertainty concerning the height of the radiological base. 


Partition of Activity Between Stratosphere and Troposphere 


The height of the tropopause varies with latitude, season, and 
daily atmospheric changes. The daily and seasonal variations are less 
in tropical latitudes than elsewhere. The tropopause heights averaged 
about 54,000 ft above sea-level for the Dominic I tests, varying between 
50,000 and 58,000 ft-on individual shot days. This variation is repre- 
sentative of the tropical tropopause. 

With the use of the activity distribution in Fig. 5, a mean tropo- 
pause height of 54,000 ft, and the mean cloud-height curve in Fig. 1, a 
typical curve of the percent of totaldebris in the troposphere as a func- 
tion of yield has been calculated. The curve, shown in Fig. 6, is intended 
to be valid at the time of cloud stabilization for air bursts in a tropical 
atmosphere. Another curve has been drawn to indicate the likely maxi- 
mum tropospheric fraction, based on a high tropopause (58,000 ft) and 
low cloud heights (from the lower curve in Fig. 1). These extremes do 
not represent absolute maximums Since higher tropopauses and lower 
clouds may occur occasionally. In addition, the uncertainties in the 
activity—height curve (Fig. 5) make it impossible to define a mean- 
ingful and useful absolute maximum tropospheric fraction. No attempt 
has been made to estimate the minimum tropospheric fraction, but, in 
the megaton-yield range, it could be several orders of magnitude below 
the typical fraction. 
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Fig. 6—Percent of total activity initially injected in the troposphere 
as a function of total yield for air bursts in a tropical atmosphere. 
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The most critical uncertainty in the estimates is for the range of 
yields from about 700 kt to about 5 Mt, where the radiological cloud 
base may lie in the vicinity of the tropopause. For yields less than 700 
kt, the tropospheric fraction (at cloud stabilization) can be estimated to 
within a factor of 2 or less. For yields above about 5 Mt, the fraction 
in the troposphere becomes very small, although precisely how small 
it may be has not yet been determined. 

An estimate of the kiloton equivalent of fission products stabilized 
_ in the troposphere as a function of total yield for air bursts is shown in 
_ Fig. 7. The typical and maximum curves are derived from the curves 
in Fig. 6, based on the assumption that the yield is entirely due to fis- 
sion. Several interesting features may be noted. The maximum tropo- 
spheric contamination is produced by bursts in the low megaton range 
(assuming 100% fission yield). With typical cloud heights and an aver- 
age tropopause height of 54,000 ft, the maximum tropospheric con- 
tamination is about 500 kt for yields between about 800 kt and 2 Mt. As 
the yield increases, the tropospheric debris decreases rapidly and then 
levels off at about 5kt of fission equivalent for yields from 10 to 100 Mt. 
The maximum curve, based on a high tropopause and low cloud heights, 
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Fig. 7—Amount of debris (kilotons equivalent) initially injected in the 
troposphere as a function of total yield for air bursts in a tropical 
atmosphere. 


is quite similar with a maximum, tropospheric contamination of about 
1.5 Mt for yields between 2 and 3 Mt, all fission. This curve also de- 
creases rapidly and then levels off at about 12 kt of tropospheric debris 
for yields between 15 and 100 Mt, It should be recalled here that these 
curves are based on the activity—height curve given in Fig. 5 and are 
subject to the same uncertainties. For thermonuclear devices the 
amount in the troposphere should be multiplied by the fission fraction 
of the device. a ye 


Dominic | Debris in the Troposphere 


With the use of the estimated tropospheric fraction and the fission | 
yield for the individual Christmas Island detonations, it is estimated 
that about one-third of the total radioactivity initially stabilized in the 
troposphere. Since a half-residence time of one month is generally 
accepted for tropospheric debris,’° about one-sixth of the total activity 
might be expected to be deposited at the surface, mostly in tropical 
latitudes, within a month after the conclusion of this test series. A rough 
integration!! of the activity found in the AEC Health and Safety Labora- 
tory monthly fallout-deposition collections indicated that only about half 
this- amount was deposited in the latitude band from 30°N to 30°S 
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through August 1962. The difference between the expected and the mea- 
sured amounts may be due to the uncertainties in the estimates of the 
tropospheric fraction and the deposition. However, there are several 
reasons for believing that the amount deposited in the latitude band was 
actually less than that.estimated from the tropospheric fraction at time 
of cloud stabilization, 

First, some of the debris which initially stabilized below the tropo- 
pause may have ascended into the stratosphere in convective cells or as 
a result of thermally induced direct circulation. Second, some debris 
was transported to mid-latitudes at altitudes below the tropical tropo- 
pause. Since there is a polar tropopause in mid-latitudes, generally 
between 30,000 and 40,000 ft, the debris which was transported away 
from the equatorial region at altitudes from about 40,000 to 55,000 ft 
would have become incorporated into the mid-latitude stratosphere. An 
‘interesting example of such transport was provided by the interception 
of one of the Dominic I clouds by sampling aircraft over the western 
United States.” In addition, sampling of the lower stratosphere over 
the United States indicates that it contained fresh debris from the 
Dominic I tests during most of the month of May 1962. Finally, the evi- 
dence for a half-residence time of one month for tropospheric debris 
may actually apply only to debris below the polar tropopause. The 
residence time for debris inthe troposphere, above 40,000 ft, in tropical 
latitudes has not been established. Only a very small fraction of the 
debris from the Dominic I tests stabilized below 40,000 ft. The fraction 
was much smaller than that for previous Pacific test series, which 
consisted primarily of surface bursts, 

In any case, it has become increasingly evident that the potential 
hazard due to short-lived fission products is not attributable solely to 
the portion initially injected in the troposphere since there is an ex- 
change of air between the stratosphere and troposphere. Therefore the 
three-dimensional trajectory of the debris-laden air would have to be 
considered in determining the fate of a particular debris cloud. 

It has also been shown” that severe thunderstorms that penetrate 
the lower stratosphere provide an effective mechanism for bringing 
stratospheric debris directly to the ground. It appears that thunder- 
storm scavenging of stratospheric debris from the Dominic I tests 
accounted for most of the ‘I found in milk in the midwestern United 
States in May 1962, 


RECOMMENDATIONS FOR FUTURE WORK 


Project Stemwinder has shown in-cloud dose-rate monitoring by 
aircraft to be a relatively simple and economical way to obtain infor- 
mation on the distribution of radioactive debris in nuclear clouds. 


644 FERBER 


Through the use of this technique, it may be possible to greatly reduce 
the amount of radiochemical analyses which would be required to ob- 
tain a reliable inventory of radioactivity in a nuclear cloud. This type 
of monitoring might be particularly well suited to clouds from large 
nuclear cratering detonations such as those contemplated for the Plow- 
share program. Such clouds should be visible for at least several hours 
and would probably be confined to the troposphere. 

Tentative answers have been found for the questions that led to 
Project Stemwinder, but large uncertainties still exist. The experience 
gained in the project indicates that for air bursts the lower stem should 
be monitored soon after cloud stabilization while it is still visible and 
that several penetrations should be made at each altitude to ensure that 
representative readings are obtained. Additional data are particularly 
needed for yields in the megaton range. 

An obvious limitation of Project Stemwinder was the aircraft ceil- 
ing of 50,000 ft. A determination of,the amount of debris initially sta- 
bilized in the troposphere requires sampling to an altitude of 60,000 ft. 
Aircraft with the capability of operating at this altitude have been used 
for cloud sampling, but they were not available for Project Stemwinder. 

The following recommendations are offered for any future opera- 
tions: 

1. A continuous-recording gamma-intensity instrument package 
with a range from 1 mr/hr to 1000 r/hr should be used for aircraft 
cloud penetrations, 

2. Experimental determination of the dose-rate reduction due to the 
aircraft should be attempted. 

3. An attempt should be made to monitor the entire cloud from a 
low- or an intermediate-yield detonation. Thus the distribution of ac- 
tivity throughout the cloud can be ascertained, and the total computed 
activity can be compared with the fission yield of the device as a 
check on the method. 
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ABSTRACT 


Marine and plant foods used by natives living on Rongelap in the 
Marshall Islands contain higher levels of long-lived fission-product 
radionuclides than do diets of people living in the United States due to 
residual contamination from fallout in 1954. During the 1963 medical 
survey of the Rongelap population, three food items indigenous to the 
Rongelap diet were brought back to the United States and consumed 
over a Seven-day period by a member of the medical team. The inges- 
tion of these foods introduced levels of Sr and '’Cs which were 20 and 
60 times higher, respectively, than in the normaldiet and was therefore’ 
considered in terms of an acute intake of two fission-product radio- 
nuclides that are important from a radiological standpoint. Urinary 
and fecal collections were analyzed separately, and whole-body '*"Cs 
measurements were made with a whole-body counter. The urine was 
the principal excretory route for the ICs. whereas the feces was the 
main removal means for the Sr. The retention of Sr could be 
represented by a Series of exponentials, whereas the retention of 
cesium as determined by whole-body counting indicated that a single 
long-term component with a biological half-life of 74 days describes 
the removal process. Reasonably good agreement was obtained be- 
tween retention as determined by whole-body counting and by excre- 
tion measurements. It is estimated that about 25% of the Sr from the 
Rongelap food was retained by the body at the end of 190 days. The 
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Fig. 5— Percent retention of 137Cs from Rongelap food. 
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Fig. 6—Percent retention of 90Sy from Rongelap food as measured 
from excretion. 


EARLY FOOD-CHAIN. KINETICS 
OF RADIONUCLIDES FOLLOWING CLOSE-IN FALLOUT 
~ FROM A SINGLE NUCLEAR DETONATION 


WILLIAM E. MARTIN ; 
University of California at Los Angeles, Los Angeles, California 


_ ABSTRACT - 


Plant samples and rabbits were collected from representative locations 
in the Sedan fallout field before and at various times after the detona- 
tion. Radiochemical and statistical analyses indicated highly significant 
correlations between estimates of gamma dose rates and maximum 
concentrations of ™Sr or '] in plant samples and in the stomach con- 
tents, bone ash, or thyroids of rabbits collected between 15 and 110 
miles from ground zero. 

The effective half-lives of ®Sr and I on fallout-contaminated 

plants were approximately 18 and 5.0 to 5.5 days, respectively. Maxi- 
mum concentrations of ®Sr in rabbit bone ash occurred about 30 days 
after the detonation and remained high for at least 60 days; but maxi- 
mum concentrations of ‘I in rabbit thyroids occurred by or before 
five days and then declined to pre-Sedan levels in less than 60 days 
after the detonation. 

Deterministic exponential models were formulated and found to 
function satisfactorily, with parameter values derived from the data, in 
providing a partial explanation of the quantitative kinetic relations be- 
tween initial concentrations of ™Sr and ''I-on plants and subsequent 
concentrations in the bone ash or thyroids of rabbits collected in the 
Sedan fallout field. Major sources of error in the estimation of input- 
parameter values and in the use of such models to make predictions 

. are described and discussed. 

Similar models were proposed for the study of radionuclide kinet- 

ics in human food chains (j.e., pasture plants, cow milk, and human 
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tissues or organs) following environmental contamination by a single 
fallout event. The results of hypothetical calculations were compared 
with the Radiation Protection Guides recommended by the Federal Ra- 
diation Council. 


INTRODUCTION 


To explain the kinetics of radionuclide transfers in a given food 
chain, it is first necessary to describe the routes and rates of transfer 
affecting specific food-chain compartments. The major route and some 
of the minor routes of radionuclide transfer and exchange in a terres- 
trial ecosystem are shown in Fig. 1. : 

Some of the fallout particles initially deposited on soil surfaces 
may be redeposited by wind or rain on plants,‘ but most of them are 
mechanically trapped and not susceptible to redeposition.” Material 
deposited on soil may enter food chains by stem-base absorption,’ by 
downward leaching in the soil profile and subsequent root uptake,*»> or 
by accidental ingestion, especially by burrowing animals.® Most of the 
fallout deposited directly on soil can be regarded as unavailable for 
rapid entry into major food chains. ; 

Some of the radioactivity contained in the fallout particles, mostly 
<44 p in diameter,’! retained on foliage may be assimilated by foliar 
absorption, *° but a much larger fraction is subject to fairly rapid re- 
moval by wind or rain.!°!1 While it remains on plant surfaces, this 
material may be ingested by herbivores. 


- ATMOSPHERE 


PLANTS ANIMALS 


INPUT = OUTPUT 


MAJOR C——> 
— 


MINOR -—-> 


Fig. 1—Major and minor routes of radionuciide transfer and exchange 
in a terrestrial ecosystem contaminated by radioactive fallout. 
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HerbivorouSs mammals may be externally contaminated by direct 
exposure to fallout or by contact with contaminated plants or soils. 
Radionuclides may accumulate in animal tissues via inhalation, which, 
in many cases, is relatively unimportant or via ingestion of con- 
taminated materials. Although inhalation, ingestion of contaminated 
soil or water, and ingestion of fallout particles while the animal is 
preening cannot be dismissed entirely, it is probably reasonable to 
assume that externally contaminated plants are the major sources of 
radionuclides for assimilation by herbivorous mammals (e.g., jack- 
rabbits or dairy cattle) feeding in a fallout-contaminated environment 
during the first 30 to 90 days after close-in fallout. 

The summer of 1962 provided an excellent opportunity to study the 
food-chain kinetics of ™Sr and ‘I in relation to desert shrubs and 
rabbits in the Sedan fallout ‘field. These studies included the formula- 
tion and testing of mathematical models that can be shown to provide 
at least a partial explanation of the kinetic relations between initial 
concentrations of ®Sr and '84I on fallout-contaminated plants and sub- 
sequent concentrations of ®Sr in the bone ash or of !*'I in the thyroids 
of rabbits collected in the Sedan fallout field. Although the results of 
these studies!!!’ are not conclusive, partly because the causes of vari- 
ation are imperfectly understood, they do provide evidence that math- 
ematical models similar to those used by radiochemists to explain 
decay-chain kinetics,'® by physiologists to explain tracer kinetics,’ 
and by the International Commission on Radiological Protection (ICRP) 
to establish maximum permissible concentrations”? can also be used to 
study food-chain kinetics under field conditions. 

The objectives of this paper are (1) to summarize some of the data 
related to ®Sr and I on desert shrubs and in rabbit tissues following 
fallout from Project Sedan, (2) to present the deterministic models that 
provide a partial explanation of these results, and (3) to illustrate the 
potential value of similar models in studying the food-chain kinetics of 
radionuclides on pasture plants, in cow milk, and in human tissues 
following a single fallout event. , 


METHODS 


Project Sedan, a peaceful nuclear-explosives test, involved the 
detonation of a 100 + 15 kt thermonuclear device at a depth of 635 ft in 
alluvium and tuff at the north end of the Nevada Test Site on July 6, 
1962. As predicted, the early fallout was relatively light and occurred 
primarily within a 150-mile sector, N60°W to N60°E, from ground zero 
in Yucca Flat. ; aor 

Before and at various times after the detonation, plant samples 
and rabbits were collected from representative. locations in the Sedan 
fallout field (see Fig. 2). Each sampling station was marked bya metal 
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Fig. 2—Approximate locations of study areas and sampling stations in 
the Sedan fallout field. Dose-rate contours (mr/hy at H+ 24) ave based 
on aerial surveys”? , 
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post to facilitate its relocation and positive identification. Plant sam- 
ples were collected by clipping twigs and foliage from the crowns of 
desert shrubs growing within a radius of 150 ft from each station 
marker. Artemisia tridentata (sagebrush) was the species collected at 
all stations in Groom Valley and the Currant Area, and Alriplex con- 
fertifolia (shadscale) was collected at all but three of the stations in 
Penoyer and Railroad valleys. Most of the rabbits representing a given 
location were shot within a few hundred yards of a station marker, but 
rabbits taken within a 1-mile radius were accepted as representatives 


of the general location. Most of the rabbits collected were blacktailed . 


jackrabbits (Lepus californicus), but a fewcottontail rabbits (Sv/vilagus 
audubont) were taken in mountainous areas where jackrabbits were not 
always available. ce 

Radiochemical analyses were made to determine the “Sr and '!I 
contents of plant samples, the "Sr content of rabbit bone ash, and the 
1317 content of rabbit thyroids. The stomach contents of rabbits were 
analyzed for I, but the samples remaining after aliquots were re- 
moved for this purpose were too small to be used in making analyses 
for ®sr. Following each determination, a decay correction was made to 
indicate the Sr or ''I concentration at the time the sample was col- 
lected. The procedures followed in making these analyses are de- 
scribed elsewhere ‘’’'®:”! in considerable detail. 

For simplicity, the only samples considered in tidepaper are those 
collected at the 20 representative locations indicated in Fig. 2. Addi- 
tional data are given in other publications.’*-'® Estimates of the initial 
gamma dose rates (Ry = mr/hr at 3 ft aboveground at H + 24) were ob- 
tained from a tracing of Guillou’s”* original large-scale map prepared 
from aerial surveys made before and after the detonation. 

Correlation and regression analyses were made to examine the 
relation between initial gamma dose rates, Ro, and initial concentra- 
ations of Sr or ‘I on plants, Py; between P, and maximum concentra- 
tions of ®sr in rabbit bone ash, By; and between Py and maximum 
concentrations of '*"I in rabbit thyroids, A;. Similar analyses were also 
made to determine the relations between S) and Ro, As, and So, where S, 
represents the initial concentration of '!I in rabbit-stomach contents. 

Estimates of *sr and ‘I retention by plants in relation to theo- 
retical deposition rates, different parts of the fallout field, and dif- 
ferent plant species were based on Eq. 1: 


a ; 
= <2 
f,=3 (i) 
where f, = retention index of N, deposited on plants, sq ft/g [N.B. This 
index multiplied by G = g (dry weight) of plants per square 
foot would give the fraction of N, deposited on plants | 
a, = pe N;/g (dry)/Ry deposited on plants 


o 
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a, = pe N,/sq ft/Rp deposited on uniform plane 
N, =a given radionuclide (e.g., ®Sr or '°'I) 


Estimates of a, were based on the results of radiochemical anal- 
yses, P,, and estimates of a, (theoretical values) were based on Eq. 2: 


9 = 0:093R.¥ (2) 
: T D,CR, 


where F, = 1.43 x 10”* fissions (?°U)/kt (Ref. 23) 

Y =percent fission yield (atoms/fission) of N, (4.79% ®Sr, 3.1% 
'31T) (Ref. 24) 

T, = radioactive half-life of N; in days (50.5-day Sr, 8.04~day 
1317) 

D, = 3.2 x 10° dis/pe/day 

C = 2.79 x 10’ sq ft/sq mile 

R,=4.5 x 104 mr/hr/kt/sq mile at 3 ft aboveground at H + 24 
(Ref. 25) 


The relations between initial gamma dose rates, Ry, and initial 
concentrations of ®Sr or 1! on fallout-contaminated plants, Pp, were 
tentatively defined by linear regression formulas. For example, where 
P, and R, are the means of Py and Ry and b,, is the coefficient of re- 
gression of Py on Ro, 


I 


Py = Py Be by, (Ro a Ry) (3) 


The deterministic exponential models used to describe and explain 
the time-specific relation between ®sr or !8!I concentrations on plants, 
P,, and of *Sr in bone ash, B,, or of '*'I in thyroids, A,, are given (if at 
t=0, P = Py, B = 0, and A = 0) by Eqs. 4 to 6: 


P, = Pye? (4) 
W. = t —~Alt 
Be PD, rfp [ eo aie »y (5) 
We \Ap— Ap Ap— Ab 
-—A,t ~Agt 
- A,= PW,f, (= Peo (6) 
* NA. Ap Apo Aa 


where t = days after fallout deposition 
P, = pe ®Sr/g (dry) or pe 'I/g (dry) on plants at t > 0 
B: = pe ®Sr/g of rabbit bone ash at t > 0 
A;,=pe !I (total) in rabbit thyroid at t > 0 
W,=g (dry) of contaminated plant material consumed per rabbit 
per day 
W, = g (dry) of bone ash per rabbit 
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f, = fraction of ingested **Sr deposited in the rabbit’s skeleton 

f, = fraction of ingested ''I deposited in the rabbit’s thyroid 

Ap = 0.693/T,, where T, = effective half-life of Sr or ‘I on plants 

A, = 0.693/T, where T, = effective half-life of *Sr in rabbit’s 
skeleton 

rd, = 0.693/T, where T, = effective half-life of ‘I in rabbit’s 
thyroid 


RESULTS 


Estimates of average gamma dose rates, of ®*Sr concentrations in 
plant samples and in the bone ash of rabbits, and of I concentrations 
in plant samples, in the stomach contents of rabbits, and in rabbit thy- 
roids are given in Tables 1 and 2. The approximate locations at which 
these samples were collected are shown in Fig. 2. 

The wide range of values and the relatively large standard errors 
of the méans given in Tables 1 and 2 indicate a high degree of appar- 
ently inherent variability.* The wide range of ™Sr and ''I concentra- 
tions in plant samples and in rabbit tissues or stomach contents is 
partly a reflection of the wide range of initial contamination levels as 
indicated by estimated gamma dose rates; but much of the inherent 
variability is undoubtedly related to the probabilistic nature of the en- 
vironmental and biological processes that influence the external con- 
tamination of plants and the accumulation of radionuclides in animal 
tissues. The probabilistic approach to the study of food-chain kinetics 
and the use of stochastic models to simulate these processes are con- 
sidered in this symposium by Turner.f 

The results of correlation and regression analyses given in 
Table 3 show that the quantitative interrelations of initial dose-rate 
estimates and estimates of maximum ®Sr or ‘11 concentrations in 
plants and animal tissues are highly significant in spite of their in- 
herent. variability. As might be expected, the correlation between 
initial gamma dose rates, Ry, and initial concentrations of ®Sr or ‘11 
on plants, Py, are highly significant. The correlations and regressions 
of plant contamination on maximum concentrations of **Sr in bone ash, 
By, and of I in thyroids, A;, are also highly significant. The some- 
what higher correlation between '*/] concentrations in stomach con- 
tents, Sp, and in thyroids, A;, suggests some difference between the 
plants collected for radiochemical analyses and those actually eaten 
by the rabbits in the Sedan fallout field. 


*Since the standard deviations (= sx ¥n) are large in relation to the means, 
the frequency distributions of these variates are sharply skewed. Preliminary 
investigations indicate that they may be lognormal rather than normal. 

~See paper by Frederick B. Turner, this volume.'5 
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Table 1— AVERAGE GAMMA DOSE RATES, Rp, AND AVERAGE CONCENTRA- 
TIONS OF Sr IN PLANT SAMPLES AND IN THE BONE ASH OF RABBITS 
COLLECTED FROM DIFFERENT PARTS OF THE SEDAN FALLOUT FIELD 
AT VARIOUS TIMES AFTER THE DETONATION* 


Initia] gamma 


Plant samples, Rabbit bone ash, 


wakes Days Cp : , BG. 
Study J eeber alee y after pe ®Sr/g (dry) pe ®Sr/g (dry) 
areas x SX n detonation x SX n x Sx n 
Groom 45.0 +£30% 5 5 4059 +30% 5 1459 + 42% 5 
Valley 15 2716 +38% 5 3648 +62% 5 
30 1544 +50% 5 3667 +41% 5 
60 788 +38% 5 3581 +74% 5 
Penoyer 16.8 +45% 5 5 948 +37% 5 1363 £51% 5 
Valley 15 470 30% 5 2024 +387% 5 
30 332 +39% 5 2283 +61% 5 
60 283 +44% 5 552 +34% 5 
Railroad 6.8 +32% 5 5 397 +14% ) 334 +30% 5 
Valley 15 269 +13% 5 620 +68% 5 
30 164 +19% 5 783 +33% 5 
60 127 +30% 5 462 +£33% 5 
Currant 1.5 +33% 5 5 318 +20% 5 295 +54% 5 
Area 15 183 + 238% 5 323 +31% 5 
30 136 +£18% 5 427 +26% 5 
60 55 +15% 5 280 +31% 5 
All areas 17.5 +£30% 20 5 14386 432% 20 863 +29% 20 
15 909 +37% 20 1680 +38% 20 
30 544 +40% 20 2097 +30% 20 
60 313 +32% 20 1389 +34% 20 
Pre-Sedan concentrations 30 +£17% 14 56 +30% 16 


*xX = mean, SX = Standard error expressed as a percentage of the mean, and 


n= number of samples. (Based on data given by Martin and Turner.'8) 
+Ro = mr/hr at 3 ft aboveground at H + 24. 


When Py or Sp, Br, and A; aregiven, the parameter values required 
to solve Eqs. 4 to 6 are Wy, Wo, Tp, Ty, Ta, fo, and f,. Some of these 
parameter values can be obtained from experimental results reported 
in the literature, and others can be obtained from the data given in 
Tables 1 and 2. For our purposes those obtained from the field data may 
be more accurate, . 

Various studies?®.?? have indicated that adult jackrabbits consume 
approximately 100 g of dry plant material per day; therefore W, = 100 g. 
Unpublished data collected during this and previous studies!:?»*8 indicate 
that the average adult jackrabbit in Nevada weighs about 2000 g and 
has a skeleton weighing 200 g (fresh). Since the ratio of fresh bone 
weight to bone-ash weight is approximately 4 tol, Wp = 50 g. The ratio, 
W,/W, = 100/50, indicates a feeding rate of approximately 2.0 g of dry 
plant material per gram of bone ash per rabbit per day. The feeding 


Table 2— AVERAGE CONCENTRATIONS OF 1311] IN PLANT SAMPLES AND IN THE 
STOMACH CONTENTS AND THE THYROIDS OF RABBITS COLLECTED FROM 
DIFFERENT PARTS OF THE SEDAN FALLOUT FIELD AT VARIOUS TIMES 

. AFTER THE DETONATION* 


Rabbit stomach 
Plant samples, contents, Rabbit thyroids, 


Sails Poe __pe !I/g (dry) pe 1I/g (dry) ne }31[ per thyroid 
areas detonation xX sk on 4 Sx no’ xX Sx n 
Groom 5 11966 +33% 5 5194 +29% 5 467 +34% 5 
Valley 10 6671 +39% 5 1283 +32% 4 386 +438% 4 
tes 3192 +34% 5 1018 +47% 5 139 452% 5 
20 1224 +49% 5 424 +28% 5 83 435% 5 
25 526 +=36% 5 161 +50% 5 61 +31% 5 
30 273 +30% 5 149 +32% 5 16 +31% 5 
Penoyer 5 1244 +338% 5 1853 +31% 5 257 +34% 5 
Valley 10 659 +30% 5 1212 +20% 4 141 414% 4 
15 412 +36% 5 638 +17% 4 88 +34% 5 
20 227 +46% 5 311 +40% 5 40 +48% 5 
25 192 +64% 5 115 +44% 5 29 452% -5 
30 81 +41% 5 81 +52% 4 1.5 +70% 5 
Railroad 5 713 411% 5 693 +41% 2 111 +65% 4 
Valley 10 376 +10% 5 604 +34% 4 85 +36% 5 
15 211 411% 5 348 +39% 3 61 +34% 5 
20 119 +17% 5 148 +31% 5 22 +45% 5 
25 96 414% 5 75 +19% 5 29 452% 5 
30 48 +25% 5 33 +24% 4 ll +£27% 5 
Currant 5 501 +15% 5 322 +37% 3 26 +27% 5 
Area 10 200 +20% 5 57 +12% 3 7 +23% 5 
15 121 457% 5 46 +28% 4 6 +50% 5 
20 56 427% 5 42 +36% 5 6 +30% 5 
25 35 +26% 5 28 +50% 4 1.6 +10% 5 
: 30 26 415% 5 7 +43% 5 ll 434% 5 
All areas 5 3606 +40% 20 2506 +29% 15 221 +28% 19 
10 1977 +44% 20 838 +21% 15 143 +34% 18 
15 984 +40% 20 554 431% 16 74 +36% 20 
20 406 +43% 20 231 +23% 20 38 +26% 20 
25 “212 +32% 20 96 +27% 19 30 427% 20 
30 113 +27% 20 69 +29% 18 12 +50% 20 
Pre-Sedan 312 457% 14 207 +51% 17 5.4 +61% 15 

concentrations 


*X = mean, sx = standard error expressed as a percentage of the mean, and 
n= number of samples. (Based on data given by Turner and Martin.'8) 


rate applicable to the accumulation of ''I in the thyroid was taken as 
100 g of dry plant material per whole thyroid per rabbit per day. 

As shown in Table 1 and Fig. 3, the apparent rate of 89Sr loss from 
fallout-contaminated plants tended to decrease with increasing time 
after the detonation, but the average effective half-life, T,, from D+ 5 
to D+ 30 (5 to 30 days after the detonation) was approximately 18 days. 
The average concentration of ®°Sr in bone ash, B,, increased rapidly 
from D + 5 to D + 15 and then more slowly to a maximum (apparent 
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Table 3—-RESULTS OF CORRELATION AND REGRESSION ANALYSES BASED 
ON ESTIMATES OF INITIAL‘-GAMMA DOSE RATES AT 20 REPRESENTA- 
TIVE LOCATIONS IN THE SEDAN FALLOUT FIELD AND ON THE CONCEN- 
TRATIONS OF 8Sr AND !3![ IN PLANT SAMPLES OR IN THE STOMACH 
CONTENTS, BONE ASH, OR THYROIDS OF RABBITS COLLECTED 
AT THE SAME LOCATIONS 


Variables* Correlatién Regression Regression formulas 
[x(y)] coefficients coefficients [k=K+byfy-y)] — 
Bsr 
Po(Ro) 0.755t 83.75 + 16.657 Py = 83.75 Rg + 335 
Bgo(Po) 0.698} 0.742 0.2507 Bg = 0.742 Py + 761 
131) 
Pg (Rg) 0.7387 ° 0.405 + 0.1877 Py = 0.405 Rg + 0.147 
Sp (Ro) 0.8357 0.184 + 0.0297 + Sp = 0.184 Rg + 0.817 : 
As(Pp) 0.626} 12.76 + 3.72t As = 12.76 Py + 122 
Ag (Sp) 0.7694 38.82 + 7.597 Ag = 38.82 Sp + 57 


*Ry = mr/hr at 3 ft aboveground at H+ 24; Pg(®Sr) = 1.26 P; (pc/g on plants 
at t= 0); Pp(!11) = 2.00 P; (nc/g on plants at t= 0); Sp(!3!I) = 2.00 S, (nc/g in 
stomachs at t =.0); Byg(®Sr) = observed values (pc/g in bone ash at t = 30 days); 
and A; (1811) = observed values [nc(total) in thyroids at t = 5 days]. 

+Statistically significant at the 1% level of probability. 
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P, = PC %5¢/G (DRY) IN PLANT SAMPLES 
By = PC 8Sr/G (DRY) IN RABBIT BONE ASH 


Fig. 3—Average concen- 
trations of Sr in plant 
samples, P,,andinthe bone 
ash of rabbits, Br, collected 
from representative loca- 
tions (see Fig. 2) in the 
Sedan fallout field at vari- 
ous times after the detona- 
tion. 
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equilibrium) about D + 30. The rate ofdecline from D+ 30 to D + 60 ap- 
proximated the radioactive half-life of **Sr. Estimates of the effective 
half-life of ®°Sr in rabbit bone (T, = 20 days) were based on the average 
rates of loss from plants and the average rates of accumulation in bone 
ash, After estimates of the other required parameter values were ob- 
tained, these and the observed values for B, were used to solve Eq. 5 
for f,. The average value thus obtained was fp = 0.0575. 

As shown in Table 2, the average concentrations of ‘I in plant 
samples (sagebrush) from Groom Valley and the Currant Area were 
higher than the. average concentrations in the stomach contents of 
rabbits from the same locations. In the Penoyer and Railroad valleys, 
1311 concentrations in the stomach contents of rabbits were somewhat 
higher than those in plant samples (mostly shadscale) from the same 
locations. Apparently the rabbits in these areas were feeding on plants 
other than sagebrush and shadscale, 

This supposition is Supported by the estimates of fp, given in Table 
4. The average f, values based on the plant samples from all locations 
are approximately the same for **Sr and *1I, However, the average 
f, value for 1317 based on concentrations in plant samples is higher than 
the average based on concentrations in the stomach contents of rabbits. 
For both **Sr and "I, the f,, values for sagebrush are higher than those 
for shadscale. The differences between f, values for **Sr and ‘I in the 
different study areas are probably due "to disparities between actual 
and theoretical, a,, deposition rates, i.e., to fractionation or to errors 
in the estimation of average gamma dose rates. In general, these re- 
sults indicate-that the plants (probably grasses and broad-leaved-herbs) 
representing the average rabbit’s diet were less efficient than sagebrush 
but somewhat more efficient than shadscale in regard to the intercep- 
tion of fallout particles. 

As shown in Table 2 and Fig. 4, the average concentrations of 1317 
in plant samples from all stations were higher than the average el 
concentrations in the stomach contents of rabbits from the same areas; 


Table 4— ESTIMATES OF f) BASED ON EQS. 1 AND 2 AND OF INITIAL 
CONCENTRATIONS OF Bgr AND !311 IN PLANT SAMPLES, Py, AND OF 
131] IN THE STOMACH CONTENTS OF RABBITS, S,* 


Study areast Plant species ®™Sr Po 1317 Py 1317 Sy 
Groom Valley Sagebrush 4.86 x<10% 5.61 «10% 2.42 x103 
Penoyer Valley Shadscalet 3.03 x10 1.56 107% 2.32 x1073 
Railroad Valley Shadscale 3.14 x103 2.21 x10 2.24 x103 
Currant Area Sagebrush 11.9 x10% 7.02 x103 4.51 x103 


Average for all locations 4.40 x10% 4,34 x10 3.02 x1073 


*Py(®Sr) = 1.26 Ps, Po(!34I) = 2.00 Ps, and S)(!341) = 2.00 Ss. 
jLocations of study areas are shown in Fig. 2. 
tSome of the plant samples from this area were Grayia spinosa. 
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Pe=PC 7G (DRY) IN 
PLANT SAMPLES 

S,=PC B/G (DRY) IN 
STOMACH CONTENTS 


Ap=NC 8 (TOTAL) IN 
RABBIT THYROIDS 
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Fig. 4—Average concentrations of 'I in plant samples, P,, in the 
stomach contents, S,, and in thyroids,A,, of rabbits collected from re p- 
resentative locations (see Fig. 2) in the Sedan fallout field at variuus 
times after the detonation. 


put the apparent rates of loss, indicating an effective half-life, T,, of 
5.0 to 5.5 days, were about the same. The average concentrations of . 
317 in rabbit thyroids attained a maximum by or before D + 5 and de- 
clined at a half-time rate of 5 to 6 days. 

French?” and French and Van Middlesworth’? did some experimental 
work with jackrabbits near the National Reactor Testing Station in 
Idaho. Their reports give estimates of T, ranging from 1.5 to 2.5 days 
and of f, ranging from 0.29 + 0.09 to 0.158 + 0.086. A more recent 
report by McBride gives estimates of T, = 2.0 days and f, = 0.29. By 
trial and error we have found that estimates of T, = 5.0 days and T, = 
2.0 days are reasonable approximations in relation to the data given in 
Table 2 and Fig. 4. With the use of these parameter values and the ob- 
served average values of A,, Eq. 6 was solved for f,. The results based 
on plant data, Pp, indicate f, = 0.18, and the results based on stomach- 
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contents data, S,, indicate f, = 0.26. (N.B. If stomach-contents data 
were available for *Sr, one might expect to find a similar difference in 
estimates of f, required to fit the observed data points.) 

Figures 5 and 6 show the relatively close agreement between ob- 
served average concentrations of ®°Sr in rabbit bone ash and of '*4I in 
rabbit thyroids and the hypothetical values obtained through solution of 
Eqs. 5 and 6. Because of variations within and between the different 
study areas (Groom, Penoyer, Railroad, and Currant), the parameter 
values that apply to the Sedan fallout field as a whole may or may not 
apply with equal accuracy to different areas withinthe fallout field. For 
example, if they were based on initial '*4I concentrations in plant sam- 
ples, estimates of ‘42 concentrations in thyroids, A,, would be higher 
than those observed in the Groom and Currant study areas but lower 
than those observed in the Penoyer and Railroad valleys, However, the 
disparities between hypothetical and observed ®*Sr concentrations in 
the bone ash of rabbits from different parts of the Sedan fallout field, 
based on Ry rather than physiography, are generally less than the 
standard errors of the observed means.'? 
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O—O OBSERVED MEANS (n= 20) 
4-—- HYPOTHETICAL MEANS BASED ON 
? 
P, =1800 PC/G Tp= 18 DAYS 
Wp= 100 G/DAY Th= 20 DAYS 
Wb= 50 G PER RABBIT (EQUATION 5) 


Fig. 5—-Observed and hy- 
pothetical average concen- 
trations of ®Sr in the bone 
ash of rabbits collected 
“from representative loca- 
tions (see Fig. 2) in the 
Sedan fallout field at various 
times after the detonation. 


PC %s-/G BONE ASH 


DAYS AFTER DETONATION 


o 
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O-——O OBSERVED MEANS (n=18 TO 20) 
4-~—- HYPOTHETICAL MEANS BASED ON 


P,=7200 PC/G (DRY) T,=5.0 DAYS 
Wp:100 G/DAY T,=2.0 DAYS 
f,=0.18 { EQUATION 6) 


TOTAL NC ! PER THYROID 


0 5 10 15 20 25 30 35 
DAYS AFTER DETONATION ; 


Fig. 6— Observed and hypothetical average concentrations of *3"I in the 
thyroids of rabbits collected from representative locations (see Fig. 2) 
in the Sedan fallout field at various times after the detonation. 


DISCUSSION 


These results seem to indicate that Eqs. 4 to 6 function satisfac- 
torily, with estimated parameter values based on observed means, in 
explaining the early food-chain kinetics of ®°Sr and ‘41 on plants and in 
the skeletons or thyroids of rabbits collected from 20 representative 
locations ‘in the Sedan fallout field at between 5 and 60 days after the 
detonation, The input parameters for Eqs. 4 to 6 (W,, W, f,, f2, Tp, Th, 
and T,) are known to vary; but the means of these values can be used, 
with fair success, as though they were constants. This deterministic 
approach to the problem permits us to describe and explain certain 
aspects of food-chain kinetics in relatively simple mathematicalterms, 
It may also permit us to make errors in estimating parameter values; 
these errors are difficult to detect because the several parameters are 
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‘interdependent. The variations within and between different study areas 
and the twofold difference between estimates of f, based on !!I con- 
centrations in plant samples and in the stomach contents of rabbits are 
good examples of the kinds of disparities to be anticipated. 

Although it is quite tempting to use such models to predict the 
biological consequences of close-in fallout, there are many good 
reasons for proceeding with caution. Some of the uncertainties, possi- 
ble sources of error, and other obstacles to the achievement of this 
goal are described and discussed below. 

From our analyses, the relation between initial gamma dose rates, 
Ro, and initial concentrations of radionuclides on fallout-contaminated 
plants, Py, appears to take the form of'a linear regression formula. 
There are several reasons for viewing this apparently simple relation 
with suspicion. For example, Guillou”? has noted that the probable 
accuracy of Ry estimates based on aerial-survey data is no greater 
than +50%. Consequently the regression coefficients given in Table 3 
could be in error by +50%, and errors associated with extrapolations 
beyond the standard regression formula limits could be even greater. 

Because of nuclide fractionation and downwind variations in the 
particle-size composition of fallout, one might not expect to find sig- 
nificant correlations between gross gamma dose rates and concentra- 
tions of specific radionuclides on plants or soils in fallout fields re- 
sulting from surface or low-altitude nuclear detonations. Previous 
studies '*"»78 along the hot lines of fallout from balloon- and tower- 
supported detonations have shown correlations between plant and 
animal contamination and the distribution of fallout particles <44 u in 
diameter.* In several cases, the highest percentages of fallout <44 yin 
diameter and the highest levels of plant and animal contamination were 
found at intermediate distances from ground zero and hence at inter- 
mediate levels of gross gamma dose rate. Under these conditions, it 
would not be possible to predict plant contamination levels by means of 
linear regression formulas involving Ry as the independent variable. 

As shown in Table 4, it is also necessary to consider those factors 
which influence a plant’s ability to intercept and retain the fallout par- 
ticles deposited on it. Other things being equal, it may be assumed that 
plants with dense foliage composed of oily, resinous, or pubescent 
leaves should have higher f, values than plants with sparse foliage 
composed of smooth or waxy leaves; but more data are needed to de- 
termine f, values for a variety of wild and cultivated plant species in 
relation to specific morphological features, the particle-size com- 
position of fallout, and various deposition rates. 


*Autoradiographs and microscopic examinations have shown that virtually 
all the radioactive particles on plant foliage in the Sedan fallout field were <50 u 
in diameter and that about half of those measured were < 20 yu in diameter. 


FOOD-CHAIN KINETICS OF RADIONUCLIDES 773 


Also needed are data to indicate the effects on f, of plant density 
(e.g.,.¢/sq ft) and other phytosociological characteristics of vegetation 
such as species composition, community structure, and phenology. In 
desert-shrub communities, for example, shrubs are widely spaced, and 
it is probably safe to assume that the f, value for plants in a given area 
is independent of plant density. In pasture or forest, where plants are 


’ more crowded and where some may grow in the shade of others, f, may 


decrease as density increases. Some evidence of this relation in regard 
to “I has been cited recently by Straub and Fooks.*! 

_ Many of these factors (fractionation, particle-size composition and 
distribution, and retention factors, a, = {,) are considered in the fallout 
model proposed by Miller. If suitable plant data were available for a 
given situation, Miller’s model would surely provide a theoretically 
sounder basis than unqualified linear regression formulas for the 
prediction of the initial concentrations of radionuclides on plants at 
different locations in a close-in fallout field. If one wanted to use the 
food-chain model described earlier (Eqs. 3 to 6) to make predictions of 
Py, Pr, Br, and Ac, a fallout model such as Miller’s should be substi- 
tuted for Eq. 3. ; 

After the deposition of fallout, the concentration of a given radio- 
nuclide on fallout-contaminated plants can be expected to decline at a 
rate significantly faster than would be predicted on the basis of its 
radioactive half-life. Our estimates of effective half-lives on plants in 
the Sedan fallout field) 18 days for “Sr and 5.0 to 5.5 days for “41, 
indicate environmental half-lives (i.e., half-time rates of loss due to 
all causes other than radioactive decay) of approximately 28 days for 
88Sr and 15 days for "I. Since there was little or no rain in the area of 
the Sedan fallout field during the period ofthis study, the environmental 
half-life of ®Sr on plants can be attributed primarily to wind action that 
removed particles from foliage or foliage from plants. The shorter 
environmental half-life of '‘I on plants may reflect the combined ef- 
fects of wind action and sublimation.''” es 

Other studies have indicated that our estimates of environmental 


‘and effective half-lives may not be applicable to other situations. For 


example, Bartlett et al.'° sprayed solutions of different fission products 
on grass that was then exposed to both wind and rain for periods up to 
60 days. Their results indicated an average environmental half-life of 
about 14 days for each of the radionuclides used. The difference be- 
tween their results and ours is probably due to the effects of rain in ~ 
removing soluble materials, but rain should have similar effects on 
particulate materials. 

Considering the large number of variables involved (e.g., plant 
species, local weather conditions, the particle-size composition of 
material deposited on plants, the developmental stage of foliage, etc.), 
one should expect to find a wide range of apparent effective half-lives 
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of different radionuclides on fallout-contaminated plants. In the absence 
of pertinent measurements, it is reasonable to assume an environ- 
mental half-life of 14 days or less for ‘I and for other radionuclides 
on ‘fallout-contaminated plants in humid regions. Except for ca an 
environmental half-life of 28 days should be more accurate for arid 
regions. In general, the effective half-life of a given radionuclide on 
plants, Ty, could be estimated, where T, is the environmental half-life 
and T, is the radioactive half-life, by 


PT. (7) 


Our method’® of estimating the effective half-life of Sr in rabbit 
bone (T, = 20 days) was made necessary by the absence of pertinent 
experimental data, and our only confidence in the accuracy of our es- 
timate is based on the results obtained (Fig. 5) when we used it to 
solve Eq. 5. Our estimate of the effective half-life of I in rabbit 
thyroids (T.2 = 2.0 to 2.5 days) was well within the range of experi- 
mental results reported by French’’ (T, = 1.5 to 2.5 days): Whenever 
possible, the effective half-lives of radionuclides in animal tissues or 
organs should be determined by experimental as well as by empirical 
methods. One can then judge which of several possible values should 
be applied to a given set of circumstances... : 

Animal retention factors, e.g., f, and f,, are especially difficult to 
evaluate; and, as shown in Eqs. 5 and 6, errors in the estimation of 
these parameter values would result in proportional errors in the 
prediction of tissue burdens. Our estimates of retention factors for 
Sr (£, = 5.75%) and for "I @, = 18% if based on Py or f, = 26% if based 
on Sj) represent mathematically arbitrary numbers calculated to ob- 
tain reasonably good fits between hypothetical and observed average 
tissue burdens. Because of the methods used in the estimation of these 
values, their physiological significance is doubtful; but these values 
may be just as useful as those obtained from feeding experiments. For 
example, French’s experimental results indicated f, values ranging 
from <10 to >30% for jackrabbits, whereas the averages obtained for 
Dutch rabbits under laboratory conditions were only half as high. We 
collected samples of sagebrush and shadscale from the Sedan fallout 
field and fed them to Dutch rabbits in the laboratory. The results 
reported by Turner" indicated f, values ranging from 2.0 + 0.6% (based 
on samples from the Currant Area) to 12.0 + 8.0% (based on samples 
from Groom Valley). Because of the difference in animal species 
(Dutch rabbits vs. jackrabbits) and because our field data indicate that 
jackrabbits in the Sedan fallout field were feeding on plants other than 
sagebrush and shadscale, these results are of dubious value in relation 
to the food-chain model. 


wo 
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DETERMINISTIC MODELS FOR HUMAN-FOOD-CHAIN KINETICS 


The data summarized and discussed in the preceding pages have 
served to illustrate the usefulness of deterministic exponential models 
in providing at least a partial explanation of the early food-chain ki- 
netics of radionuclides following a single fallout event. Although the 
results of these studies are not conclusive, they are promising and 
provide a basis for the supposition that similar models should be use- 
ful in the study of radionuclide transfers in food chains leading to man. 

For example, suppose that an ordinary pasture is contaminated by 
fallout from a single nuclear detonation or from a reactor accident and 
that the milk produced by dairy cattle grazing in the contaminated 
pasture is consumed by people living on the farm or in a nearby vil- 
lage. If an estimate of the initial concentration of a given radionuclide, 
N;, on pasture plants, P),-can be obtained from direct measurement or 
predicted by a suitable fallout model, the subsequent concentration of 


‘N, on pasture plants, P,, in cow milk, M,, and in the human tissue or 


organ of reference, H;, can be estimated, with the assumptions that. 
P = Py, M = 0, and H = 0 at t = 0, by Eqs. 8 to 10: 


P, = Pye” pt 4 
—Apt dant 
My = PoKinfm (4 : = =| i 
i emt eT ‘mt 
ee epee acres (creer 
4 e (10) 


GQ, aaow) 


where t = days after fallout 

N, = a given radionuclide (e.g., Sr or ‘**I) 

P, = pe N,/g (dry) on pasture plants at t = 0 

P,=pe N,/g (dry) on pasture plants at t > 0 

M,= pe N,/ml (fresh) in milk produced at t > 0 

i = pe Ni/g (fresh)in the human tissue or organ of reference att>0 

K,, =dry weight (g) of pani consumed/volume (m1) of milk pro- 
duced per day 

K,,= volume (ml) of milk consumed per day/weight (g) of human 
organ or tissue 

= fraction of ingested N,; secreted in cow milk 
: =fraction of ingested N;, deposited in human tissue or organ of 


reference 
A = 0.693/T 
.T, = effective half-life of N, on pasture plants 
T, = effective half-life of N; in cow milk production 


-T,, = effective half-life of N, in human tissue or organof reference 


776 MARTIN 


To illustrate the potential value of the model formulated above, we 
have adopted the parameter values listed in Table 5. Most of these 
values are based on experimental data reported in the literature; but, 
when two or more estimates of a given parameter value were found, we 
usually selected the one resulting in-higher concentrations in milk or 


3 


Table 5—-PARAMETER VALUES USED WITH EQS. 8 TO 10 TO 
CALCULATE HYPOTHETICAL CONCENTRATIONS OF 89Sr 
AND %!] ON PASTURE PLANTS, IN COW MILK, AND IN HU- 
MAN TISSUES FOLLOWING A SINGLE FALLOUT EVENT* 


Parametert ®Sr 131] References 


Py 100 pe/g 400 pe/g Arbitrary values 
Km 1.4 g/ml/day 1.4 g/ml/day 33 

fe 0.02 0.06 34, 35 

K,, (adult) 1/7 ml/day/g 50 ml/day/g 20 

K ;, (infant) 10/7 ml/day/g 500 ml/day/g 20 

fy G:21* 0.30 20 

Krt 1.43 x 107+ 1,17 x 107 20 

Ty 18.0 days 5.5 days 13, 16 

hes 2.5 days 2.0 days 34 

Ty 50.4 days 7.5 days 20 


*Results are shown in Figs. 7 and 8. 
+Parameters are defined in the text. 
tUsed with Eq. 11 to calculate tissue doses. 


human tissues. The results of calculations based on Eqs. 8 to 10 and 
the parameter values listed in Table 5 are shown in Figs. 7 and 8. Ob- 
viously, for different values of Py, the corresponding values of M, and 
H, would be proportionally larger or smaller than shown in Figs. 7 and 
8. Using the same effective half-lives but different estimates of K,,, 
f.,, Ky, or f, would also result in proportionally higher or lower esti- 
mates of M, and H,. 

The total dose, D,, delivered by N; to the human tissue or organ of 
reference is given, fort =0 tot = %, by 


1—e7*r 1—e sn 
D+ = Py9Kinfm Kuf.K, |_—————- + — 
PE [Nevin Ne) Mh a Ap). Men Op = Nin) n= Naa) 
1 — eT Mat 


u AWA = Ay) Og, — Ay) a 
where D; = total dose (rem) delivered by Ni, from t = 0 to t > 0, to the 
human organ or tissue of reference ‘ 
K, = [(3.20 x 10° dis/pe/day) (E)/(6.24 x 10" Mev/100 ergs/g)| 
E = effective absorbed energy (Mev) per disintegration of N, 
[Mev .absorbed/disintegrations of N; x relative biological 
effectiveness (RBE)| 


1, =PC ®S:/KG ON PASTURE PLANTS 
M,= PC 51/10 LITERS OF COW MILK 
Hy=PC Sr (TOTAL) IN HUMAN SKELETON 
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PC 


10° 


107 
0 20 40 60 80 100 120 140 
DAYS AFTER FALLOUT 


Fig. 7—Hypothetical concentrations of ®Sr on pas- 
ture plants, in cow milk, and in human skeletons fol- 
lowing environmental contamination by a single fallout 
event. 


ie" | T 


P, =PC 8) /KG ON PASTURE PLANTS 
My= PC UL / LITER IN COW MILK 
Hy = PC 311 (TOTAL) IN HUMAN THYROID 


| 
0 a 10 1S 20 25 30 35 
DAYS AFTER FALLOUT 


Fig. 8—Hypothelical concentrations of *'I on pasture 
plants, in cow milk, and in human thyroids following 
environmental contamination by a single fallout event. 
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. If the concentration of N; on pasture plants is not known, estimates 


of doses to human tissue can be based on concentrations in milk (M, = 


pe N,/ml). If the time of fallout is known and M, is given for a spec- 
ified time of milking, Eq. 9 can be simplified, as shown in Eq. 12, to 


obtain an integration constant, I), that can then be substituted for the © 


expression P)K,f,, in Eqs. 10 and 11 to obtain estimates of concentra- 
tions in, H,, and doses to, D,, the human tissue or organ of reference: 


er t — ert 
M, = 2—2= (12) 
m Pp 


James*® recently used a similar method to estimate human thyroid 
doses resulting from a Single fallout event, but many of his parameter 
values (i-e., f, = 0.25, K,= 1.0 x 107°, T, = 5.0 days, and T,, = 1.9 days) 
were not as pessimistic as those given in Table 5. James concluded 
from his investigation that a maximum concentration of 2800 pc of 13/1 
per liter of milk resulting from a single fallout event would, if a milk 
consumption of 1 liter/day were assumed, result in a total dose of 0.5 
rad to a 2,0-g thyroid. Our results (Table 6) indicate that a maximum 
concentration of only 1850 pc of ‘“!I per liter of milk would be required 


for a potential dose of 0.5 rem to a 2.0-g human thyroid. (N.B. In this ° 


case 1 rem = 1 rad x RBE, and RBE = 1.0; therefore 1 rem = 1 rad.) 

Pendelton et al.* based their estimates of thyroid doses on esti- 
mates of total ‘I intake, I,, by people consuming 1 liter of milk per 
day in various parts of Utah during the summer of 1962. According to 
our hypothesis, the value of I, would be given by 


I, = IV; (13) 


ApQAm — Ap) m\“p m) 


where I) is an integration constant obtained from .Eq. 12 and V,, is the 
volume of milk in milliliters consumed per person per day, 

With the use of I, from Eq. 13, the total thyroid dose would be 
given, where W, is the thyroid weight (fresh) in grams, by 


os I,K, fp 


Ber AnWh 


(14) 


The Radiation Protection Guides (RPG’s) recommended by the 
Federal Radiation Council (FRC)*® are said to represent “...a rea- 
sonable balance between biological risk and benefit to be derived from 
useful applications of radiation and atomic energy.” The RPG’s for 
human bone and thyroid are based on average rates of radiostrontium 
. and radioiodine intakes which, in the opinion of the FRC, should result 
in doses no greater than 1.5 rem/year to “... individuals in the gen- 


; 
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Table 6—-SUMMARY OF HYPOTHETICAL VALUES THAT, IF INDICATED 
BY MEASUREMENTS MADE AFTER ENVIRONMENTAL CONTAMINA- 
TION BY A SINGLE FALLOUT EVENT, WOULD IMPLY TOTAL DOSES 

OF 0.5 REM TO THE SKELETONS OR THYROIDS OF INFANTS 
CONSUMING 1 LITER OF MILK PER DAY 


Hypothetical values* 89Sr 1317 
Initial concentrations on pasture plants, Pg 61.0 pe/g ' 13.7 pe/g 
Maximum concentrations in milk, M, 4500 pe/liter 1850 pe/liter 
Time after fallout, t wax 8 days 4 days 
Total intake (to t = ~) : 1.60 x 108 pe 2.63 x 104 pe 
Maximum concentration in human tissue, H; 27 pe/gt 1580 pe/gt 
Time after fallout, t pax 50 days 15 days 
Total dose (at t = ©) . 0.5 remj 0.5 remt 


*Calculations were based on Eqs. 8 to 14 and the parameter values listed 
in Table 5. 
+Based on a 700-g skeleton or a 2.0-g thyroid. 


eral population...” or an average of 0.5 rem/year “... to be applied 
to suitable samples of an exposed population group.” Average annual 
intakes of >2000 pe *Sr/day or >100 pe ''I/day “... would be presumed 
to result in exposures exceeding the RPG... .” 

With the assumptions of a skeleton weight of 700 g and a thyroid 
weight of 2.0 g, Eq. 14 and the parameter values listed in Table 5 indi- 
cate that total annual intakes of 160,000 pc of **Sr or of 26,300 pc of 
817 could result in doses of 0.5 rem/year to the bones and the thyroids 
of infants, These values are lower than thoSe indicated by the FRC 
(i.e., 720,000 pc **Sr/year and 36,500 pe '4I/year) because the param- 
eter values in Table 5 are more pessimistic than those adopted by the 
FRC, — 

For further comparison with the FRC’s recommendations, we have 
estimated the hypothetical values of various measurements which, if 
obtained following a single fallout event, could result in doses of 0.5 
rem to the bones and thyroids of infants. These hypothetical values are 
given in Table 6. 

Tables 1 and 2 show that initial concentrations of ‘I on plants in 
various parts of the Sedan fallout field were about four times higher 
than the initial concentrations of ®*Sr. The estimates given in Table 6 
indicate that initial concentrations of ®°Sr about 4.5 times higher than 
those of ‘*'I are required to deliver comparable doses to human tissues. 
We might therefore suppose that ‘I is considerably more hazardous 
than **sr in an area contaminated by close-in fallout. 

The data given in Table 2 indicate that the initial concentration of 
317 on plants in the vicinity of Currant, Nev., was approximately 1000 + 
150 pc/g. Similar plant samples collected in a hot spot (centered near 
Fruitland, Utah), which was discovered after the Small Boy detonation 
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on July 14, 1962, indicated an initial ‘I concentration (on July 15) of 
approximately 800 + 200 pc (n-= 6). If we assume an f, value of about 
3.0 x 107? for pasture plants compared to aopresiniabely 7.5 x 107 for 
desert shrubs, the initial concentrations of '*'I on pasture plants in the 
vicinity of Currant, Nev., and Fruitland, Utah, could have been 400 and 
320 pe/g, respectively. 

If our hypothetical pastures, dairy herds, and human populations 
had been located in these areas, the total '°4l intakes by people drinking 
1 liter of milk per day could have ranged from 670,000 to 825,000 pc, 
and the doses to the thyroids of children could have been 11.4 to 14.6 
rem. With the use of the more optimistic parameters adopted by the 
FRC, the estimated doses based on these intakes would be 9.2 and 11.3 
rem. (N.B. These estimates are comparable to the highest estimates 
reported by Pendelton et al.*’ for stations located in other parts of 
Utah), 

Actually, there are very few milk cattle in the area of the Sedan 
fallout field, and, since we collected no milk samples, our treatment 
of the problem is strictly hypothetical. However, one sample collected 
near Ely, Nev., and analyzed by the U. S. Public Health Service®’ con- 
tained 2800 pc '*4I/liter on July 24, 18 days after the detonation. If this 
value represents the concentration at the time of milking, it could 
indicate a possible thyroid dose of 2.6 rems. 

Solutions to Eq. 13 indicate that ‘I intakes during periods of 7, 
14, or 21 days following environmental contamination by a single fall- 
out event should account for 42, 85, or 97%, respectively, of the total 
potential intakes from t = 0 tot =~. Therefore the simplest counter- 
measures to avoid 85% or more of the potential biological hazard 
related to ‘I would be (1) to feed cattle on stored feed for a period of 
two or more weeks after a detonation or (2) to use the milk produced 
during that period for making cheese or other dairy products that 
would not be consumed for a period of three or more weeks after milk 
production. 
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